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CHAPTER  1 


INTRODUCTION 


1.1.  Dynamic  Data  Structures 

In  recent  years,  the  design  of  dynamic  data  structures  has  received  increasing  attention,  and 
is  one  of  the  most  challenging  areas  of  theoretical  computer  science. 

As  a  typical  example,  consider  search  problems  where  a  set  5  of  objects  is  given.  We  want 
to  answer  queries  about  these  objects  in  an  on-line  fashion,  i.e.,  the  queries  are  not  known  in 
advance,  and  the  answer  to  a  query  must  be  given  before  the  next  query  can  be  processed.  Typi¬ 
cally,  the  set  S  is  represented  by  means  of  a  data  structure  D  on  which  the  query  operations  can 
be  efficiently  performed.  The  complexity  of  such  a  data  structure  is  usually  measured  in  terms 
of  the  space  required  for  storing  D,  the  preprocessing  time  needed  to  construct  D,  and  the  query 
tune.  In  several  applications  the  set  5  is  not  static,  but  dynamically  evolves  in  time.  This  can  be 
modeled  by  means  of  update  operations  that  modify  S  and  can  be  arbitrarily  interleaved  with  the 
queries.  Therefore,  the  problem  arises  of  designing  a  dynamic  data  structure  that  supports  on¬ 
line  both  queries  and  updates  on  5.  In  this  case,  the  previous  complexity  measures  are  supple¬ 
mented  by  the  time  bounds  for  the  update  operations. 

The  classical  dictionary  problem  is  the  simplest  example  of  a  dynamic  search  problem. 
Here,  S  is  a  subset  of  an  ordered  set  U,  the  queries  ask  whether  a  given  element  of  U  is  in  S,  and 
updates  are  insertions  and  deletions  of  elements.  For  this  problem,  the  earliest  example  of  an 
efficient  dynamic  data  structure  is  the  AVL  tree  [1],  which  requires  linear  space  and  supports 
queries,  insertions,  and  deletions  in  logarithmic  time. 

Our  investigations  focus  on  the  problem  of  searching  in  a  subdivision  of  the  plane,  a  funda¬ 
mental  issue  in  computational  geometry,  for  it  represents  the  2-dimensional  extension  of  the  dic¬ 
tionary  problem.  Three  specific  problems  have  been  addressed  in  this  area.  The  first  problem, 
dynamic  point  location,  considers  a  geometric  subdivision  of  the  plane  into  polygonal  regions. 
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and  asks  for  the  region  that  contains  a  given  query  point  The  second  problem,  dynamic  planar 
embedding,  considers  a  topological  subdivision  of  the  plane  induced  by  a  planar  embedding  of  a 
graph  G,  and  asks  for  a  region  that  contains  two  given  query  vertices  of  G  on  its  boundary  (if 
one  exists).  The  third  problem,  dynamic  transitive  closure  considers  a  planar  acyclic  digraph  G 
embedded  in  the  plane,  and  asks  for  testing  the  existence  of  or  reporting  a  directed  path  between 
two  query  vertices  of  G.  In  all  of  the  three  problems  the  update  operations  consist  of  inserting  or 
deleting  vertices  and  edges. 

L2.  Overview 

This  thesis  is  conceptually  subdivided  into  two  parts.  In  the  first  part,  consisting  of 
Chapters  2,  3,  4,  and  3,  we  consider  searching  in  a  geometric  structure,  and  present  three 
dynamic  techniques  for  planar  point  location.  In  the  second  part,  consisting  of  Chapters  6  and 
7,  we  consider  searching  in  a  topological  structure,  and  present  dynamic  techniques  for  main¬ 
taining  the  planar  embedding  and  the  transitive  closure  of  a  planar  graph. 

In  Chapter  2,  we  provide  basic  geometric  definitions  and  summarize  previous  research  on 
planar  point  location,  a  fundamental  geometric  search  problem  which  is  used  as  a  subroutine  in  a 
variety  of  algorithms.  Most  of  the  past  research  on  planar  point  location  has  been  targeted  to  the 
static  case,  where  the  subdivision  is  fixed,  and  point  location  queries  have  to  be  answered  on¬ 
line.  For  this  instance  of  the  problem  several  efficient  methods  are  available.  The  analogy  with 
one-dimensional  search,  for  which  both  static  and  dynamic  optimal  techniques  have  long  been 
known,  naturally  motivates  the  desire  to  develop  techniques  for  dynamic  planar  point  location, 
where  the  planar  subdivision  can  be  modified  by  insertions  and  deletions  of  points  and  segments. 
Work  on  dynamic  point  location  is  a  rather  new  undertaking,  and  in  a  recent  paper  [50]  Samak 
and  Taijan  indicate  as  one  of  the  most  challenging  problems  in  computational  geometry  the 
development  of  a  fully  dynamic  point  location  data  structure  whose  space  and  query  time  perfor¬ 
mance  are  of  the  same  order  as  that  of  the  best  known  static  techniques  for  this  problem.  As 
shown  in  the  following,  we  come  very  close  to  this  goal. 

In  Chapter  3,  we  present  a  dynamic  point  location  technique  for  monotone  subdivisions. 
The  update  operations  are  the  insertion  of  a  vertex  on  an  edge  and  of  a  chain  of  edges  between 
two  vertices,  and  their  reverse  operations.  Let  n  be  the  number  of  vertices  of  the  subdivision. 
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The  data  structure  uses  space  O  ( n ).  The  query  time  is  O  (log^n),  the  time  for  insertion/deletion 
of  a  vertex  is  0(logn),  and  the  time  for  insertion/deletion  of  a  chain  with  k  edges  is 
OOog2/!  +k),  all  worst  case.  This  technique  is  based  on  the  chain  method  of  Lee-Preparata  [35]. 
The  emergence  of  full  dynamic  capabilities  is  afforded  by  a  subtle  choice  of  the  chain  set 

(separators),  which  induces  a  total  order  on  the  set  of  regions  of  the  planar  subdivision. 

r. 

In  Chapter  4,  we  present  a  dynamic  point  location  technique  for  a  class  of  monotone  subdi¬ 
visions  (which  includes  convex  subdivisions),  whose  vertices  lie  on  a  fixed  set  of  N  horizontal 
lines.  The  supported  update  operations  are  insertion/deletion  of  vertices  and  edges,  and  (hor¬ 
izontal)  translation  of  vertices.  This  technique  achieves  query  time  0(logn  +  log  AO,  space 
0(n  log /V),  and  insertion/deletion  time  0  (log  n  log  N),  where  n  is  the  number  of  vertices  of  the 
subdivision.  Hence,  whenever  N  =  0(n),  the  query  time  is  0(logn),  which  is  optimal.  This 
r  *  technique  is  based  on  the  trapezoid  method  of  Preparata  [46],  which  has  been  experimentally 

shown  to  be  the  fastest  point  location  method  among  those  with  asymptotically  optimal  query 
time  [12].  It  is  easily  realized  that  in  many  significant  applications  the  most  frequent  operation 
is  point  location  query,  while  updates  are  more  rarely  executed.  Hence,  this  technique  provides 
®  the  most  efficient  solution  for  such  applications. 

• 

In  Chapter  5,  we  present  a  dynamic  point  location  technique  for  triangulations,  which  is 
based  on  the  incremental  reconstruction  method  of  Bentley-Saxe  [3]  and  Overmars  [42]  for 
I  decomposable  searching  problems.  It  achieves  O  (log2 n /log  log  n )  query  time,  O(n)  space 

requirement,  and  O  (log3n/log  logn)  update  time,  where  n  is  the  number  of  vertices.  This  tech¬ 
nique  can  be  used  in  conjunction  with  any  static  point  location  algorithm. 

In  Chapter  6,  we  present  a  dynamic  technique  for  maintaining  on-line  the  transitive  closure 

*  _ 

of  a  planar  sr-graph,  i.e.,  a  planar  acyclic  digraph  embedded  in  the  plane  with  exactly  one  source, 
s,  and  one  sink,  t,  both  on  the  external  face.  This  class  of  graphs  was  first  introduced  in  the 
planarity  testing  algorithm  of  Lempel  et  al.  [36],  and  was  fruitfully  used  in  a  number  of  applica¬ 
tions,  which  include  planar  graph  embedding,  graph  planarization,  graph  drawing,  floor  plan¬ 
ning,  planar  point  location,  visibility  representations,  motion  planning,  and  VLSI  layout  compac¬ 
tion.  Also,  planar  rr- graphs  are  important  in  the  theory  of  partially  ordered  sets,  since  they  are 
I  associated  with  planar  lattices  [32].  We  show  that  a  planar  sr-graph  G  admits  two  total  orders 
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(called  leftist  and  rightist,  respectively)  on  the  set  V\jE  uF,  where  V ,  E,  and  F  are  respectively 
the  set  of  vertices,  edges,  and  faces  of  G,  with  I V  I  =  n.  Assuming  that  G  is  to  be  dynamically 
modified  by  means  of  insertions  of  edges  and  expansions  of  vertices  (and  their  inverses),  we 
exhibit  a  O  (n)-space  dynamic  data  structure  for  the  maintenance  of  these  orders,  such  that  an 
update  can  be  performed  in  time  O  (log  n).  From  this  result,  we  derive  a  dynamic  data  structure 
for  the  transitive-closure  query  problem,  which  consists  of  testing  for  the  existence  of  (or  report¬ 
ing)  a  directed  path  between  two  vertices  u  and  v  of  G.  The  space  is  O(n),  and  the  time  for 
queries  and  updates  is  O  (log  n)  (worst-case).  Also,  we  show  that  the  discovered  dynamic  pro¬ 
perties  of  planar  sr-graphs  provide  the  topological  underpinning  of  the  planar  point  location  tech¬ 
nique  for  monotone  subdivisions  presented  in  Chapter  3,  and  can  be  applied  to  the  problem  of 
dynamic  contact  chain  queries. 

In  Chapter  7,  we  present  a  dynamic  technique  that  allows  for  incrementally  constructing  a 
planar  embedding  of  a  planar  graph  with  n  vertices.  This  problem,  referred  to  as  dynamic 
embedding  problem ,  naturally  arises  in  interactive  CAD  layout  environments.  A  query  operation 
asks  whether  a  new  edge  (u,v)  can  be  added  to  the  embedding  without  introducing  crossings, 
i.e.,  whether  there  is  a  face  of  the  embedding  whose  boundary  contains  both  u  and  v.  The  update 
operations  are  the  insertion  and  deletion  of  vertices  and  edges.  We  present  a  data  structure  for 
the  dynamic  embedding  problem  that  uses  0(n)  space  and  supports  queries  and  updates  in  time 
G(Iogn)  (worst-case).  These  results  are  obtained  by  maintaining  on-line  an  orientation  of  the 
graph,  called  spherical  st-orientation,  a  notion  that  generalizes  planar  sr-graphs.  This  technique 
has  applications  to  circuit  layout,  graphics,  motion  planning,  and  computer-aided  design.  It  also 
constitutes  a  first  step  toward  the  development  of  an  efficient  algorithm  for  testing  on-line  the 
planarity  of  a  graph,  under  insertions  and  deletions  of  vertices  and  edges. 

As  a  general  remark,  we  would  like  to  underscore  that  in  addition  to  their  good  theoretical 
space/time  performance,  all  the  data  structures  and  algorithms  presented  in  this  thesis  are  also 
practical  and  easy  to  implement,  and  therefore  suited  for  real-world  applications. 

The  material  in  the  thesis  is  presented  so  that  each  chapter  can  be  read  independently  from 
the  others.  Chapter  2  provides  background  material  for  Chapters  3,  4,  and  5.  For  the  reader’s 
convenience,  some  definitions  are  repeated  wherever  appropriate. 
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CHAPTER  2 


PLANAR  POINT  LOCATION 

2.1.  Planar  Subdivisions 

The  motivation  for  the  geometric  definitions  given  below  is  readily  obtained  if  we  view  a 
point  of  the  plane  as  the  central  projection  of  a  point  of  a  hemisphere  to  whose  pole  this  plane  is 
tangent. 

A  vertex  v  in  the  plane  is  either  a  finite  point  or  a  point  at  infinity  (the  latter  is  the  projection 
of  a  point  on  the  hemisphere  equator).  An  edge  e  =(u,v)  is  the  portion  of  the  straight  line 
between  u  and  v,  with  the  only  restriction  that  u  and  v  are  not  points  at  infinity  associated  with 
the  same  direction.  Thus,  e  is  either  a  segment  or  a  straight-line  ray,  but  not  a  whole  straight 
line.  When  both  u  and  v  are  at  infinity,  then  e  is  an  edge  at  infinity,  i.e.,  a  portion  of  the  line  at 
infinity  (the  projection  of  an  arc  of  the  equator).  A  {polygonal)  chain  y  is  a  sequence 
(«,•  :«i  =  (Vj,vi  +  1),i  =  1,  •  •  •  ,p  - 1)  of  edges;  it  is  simple  if  nonintersecting;  it  is  monotone  if  any 
line  parallel  to  the  x-axis  intersects  y  in  either  a  point  or  a  segment  In  the  following,  the  notion 
of  “left”  and  “right”  refer,  when  not  specified  otherwise,  to  a  bottom-up  orientation  of  the 
entity  being  considered  (a  chain,  or,  later  on,  a  separator,  an  edge,  etc.). 

A  simple  polygon  r  is  a  region  of  the  plane  delimited  by  a  simple  chain  with  vp  =  v  i ,  called 
the  boundary  of  r.  Note  that  r  could  be  unbounded;  in  this  case  the  boundary  of  r  contains  one 
or  more  edges  belonging  to  the  line  at  infinity.  A  polygon  r  is  monotone  if  its  boundary  is  parti- 
tionable  into  two  monotone  chains  yi  and  Y2,  called  the  left  chain  and  right  chain  of  r,  respec¬ 
tively  (see  Fig.  2.1).  Chains  yi  and  %  share  two  vertices,  referred  to  as  HIGH(r)  and  LOW(r), 
respectively,  with  y{HIGH{r))>y{LOW{r)).  In  other  words,  HIGH(r)  and  LOW(r)  are, 
respectively,  points  of  maximum  and  minimum  ordinates  in  polygon  r,  each  of  them  is  unique 
unless  it  is  the  extreme  of  either  a  horizontal  edge  or  an  edge  at  infinity,  in  which  case  the  selec¬ 
tion  between  the  two  edge  extremes  is  arbitrary.  A  convex  polygon  is  a  simple  polygon  r  such 
that  for  any  two  vertices  u  and  v  of  r  the  segment  with  endpoints  u  and  v  lies  entirely  inside  r.  A 
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HIGH  (r) 


Figure  2.1  Nomenclature  for  a  monotone  polygon, 
convex  polygon  is  also  a  monotone  polygon. 

A  planar  subdivision  9?  is  a  partition  of  the  plane  into  simple  polygons,  called  the  regions 
of  91.  It  is  easily  realized  that  91  is  determined  by  a  planar  graph  G  embedded  in  the  plane 
whose  edges  are  either  segments  or  rays  of  straight  lines  (referred  to  as  a  “planar  straight-line 
graph”  in  [35]):  edges,  vertices,  and  chains  of  G  are  referred  to  as  edges,  vertices,  and  chains  of 
91.  The  vertices  of  91  are  both  the  finite  ones  and  those  at  infinity.  This  ensures  the  validity  of 
the  well-known  Euler’s  formula  and  its  corollaries: 

I VI  +  I/?  I  =  l£l  +2;  IE  I  £3  I VI  -6;  IE  I  £2  I VI -4, 
where  V,  E,  and  R  are  respectively  the  set  of  vertices,  edges,  and  regions  of  91. 

A  monotone  subdivision  91  is  a  planar  subdivision  whose  regions  are  monotone  polygons 
(see  Fig.  2.2(a)).  Monotone  subdivisions  form  an  important  class  of  planar  subdivisions,  since 
they  include  convex  subdivisions ,  whose  regions  are  convex  polygons  (see  Fig.  2(b)),  and  tri¬ 
angulations,  whose  regions  are  triangles  (see  Fig.  2(c)). 


2.2.  Point  Location  Techniques 

The  point  location  problem  is  formulated  as  follows:  Given  a  planar  subdivision  SR  with  n 
vertices,  determine  to  which  region  of  SR  a  query  point  q  belongs.  If  the  query  point  q  is  on 
some  vertex/edge  of  SR,  then  that  vertex/edge  is  reported.  The  repetitive  use  of  SR  and  the  on¬ 
line  requirement  on  the  answers  call  for  a  preprocessing  of  51  that  may  ease  the  query  operation, 
analogous  to  sorting  and  binary  search  in  one-dimensional  search.  The  history  of  planar  point 
location  research  spans  more  than  a  decade  and  is  dense  in  results:  the  reader  is  referred  to  the 
extensive  literature  on  this  subject  [1 1-13, 34, 35, 37, 46, 47, 50]. 

Most  of  the  past  research  on  the  topic  has  focused  on  the  static  case  of  planar  point  loca¬ 
tion,  where  the  subdivision  SR  is  fixed.  For  this  instance  of  the  problem,  several  practical  tech¬ 
niques  are  available  today  (e.g.,  [12, 13,35,46,50]),  some  of  which  have  O(n)  space  require¬ 
ment  and  O  (log n)  query  time  [13, 50],  and  are  thus  provably  optimal  in  the  asymptotic  sense. 

Work  on  dynamic  point  location  is  a  rather  recent  undertaking.  Overmars  [43]  proposed  a 
technique  for  the  case  where  the  n  vertices  of  51  arc  given,  the  boundary  of  each  region  has  a 
bounded  number  of  edges,  and  only  edges  can  be  easily  inserted  or  deleted.  The  basic  entities 
used  in  Overmars’  method  are  the  edges  themselves;  each  edge  currently  in  51  is  stored  in  a  seg¬ 
ment  tree  defined  on  the  fixed  set  of  vertex  abscissae,  and  the  edge  fragments  assigned  to  a  given 
node  of  the  segment  tree  form  a  totally  ordered  set  and  arc  therefore  efficiently  searchable.  This 
approach  yields  O  (/ilog/i)  space  requirement,  and  O  (log2/!)  query  and  update  times.  Note  that 
these  measures  are  unrelated  to  the  current  number  of  edges  in  51.  This  technique  can  be 
extended  to  support  insertions  and  deletions  of  vertices  in  O  Gog2/!)  amortized  time,  at  the 
expense  of  deploying  a  rather  complicated  and  not  very  practical  data  structure. 

Another  interesting  dynamic  point  location  technique,  allowing  both  edge  and  vertex 
updates,  is  presented  by  Fries  and  Mehlhom  [18, 19,39],  Their  approach  achieves  O(n)  space 
requirement,  Oflog2/!)  query  time,  and  0(log4/t)  amortized  update  time.  If  only  insertions  are 
considered,  the  update  time  is  reduced  to  0  (log2/!)  ( amortized)  [39,  pp.  135-143]. 

Tables  2.1  and  2.2  summarize  the  performance  of  the  aforementioned  static  and  dynamic 
point  location  techniques.  The  new  results  on  dynamic  point  location  presented  in  Chapters  3, 4, 
and  5  of  this  thesis  are  included  in  Table  2.2. 
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Table  2.1  Static  point  location  techniques 


E£9 

Query  Time 

Source 

Notes 

«2 

log  ft 

Dobkin-Lipton  76 

slab  method,  not  practical 

n 

log2  ft 

Lee-Preparata  76 

chain  method,  practical 

« 

log  ft 

Lipton-Tarjan  77 

not  practical 

nlogn 

log  ft 

Preparata  81 

trapezoid  method, 
practical 

n 

log  ft 

Kirkpatrick  83 

not  practical 

#i 

log  ft 

Edelsbrunner  et  aL  85 

chain  method  with 
fractional  cascading, 
practical 

n 

log  ft 

Samak-Tarjan  85 

slab  method  with 
persistent  search  trees, 
practical 

Table  23.  Dynamic  point  location  techniques.  Notation: 
S(n):  space;  Q(n):  query  time  Ie(n):  time  for  edge 
insertion;  Iv(n):  time  for  vertex  insertion;  De(n):  time  for 
edge  deletion;  Dv(n):  time  for  vertex  deletion. 


ESI 

Q(n) 

I.M 

warn 

DM) 

Source 

Notes 

nlogn 

lOg2/! 

lOg2#! 

lOg2#! 

lOg2#! 

lOg2/! 

Overmars  85 

fixed  vertices,  regions 
with  0(1)  edges,  not 
practical  vertex 
updates 

■ 

lOg2/! 

lOg4#! 

k>g4ft 

k)g4#! 

lOg4#! 

Fries-Mehlhom  85 

amortized  update 
times 

■ 

lOg2#! 

lOg2#! 

lOg2#! 

— 

— 

Fries-Mehlhom  84 

only  insertions, 
amortized  update  time 

■ 

lOg2/! 

log 2n+* 

lOgft 

lOg2#!** 

lOgft 

This  Thesis 
Chapter  3 

monotone 
subdivision, 
insertions/deletions  of 
k-edge  chains 

ftlogN 

logn+logN 

log#!  log  A/ 

logn  logN 

log/!  log// 

_ 

log#!  log jV 

This  Thesis 
Chapter  4 

convex  subdivision, 
vertices  on  N  fixed 
horizontal  lines 

n 

Knai 

BToi 

wmm 

warm 

msrjm 

This  Thesis 
Chapter  5 

triangulations 

log  log  ft 

log  log  ft 

log  logn 

log  log  ft 

log  log  ft 

■ 

1  . 

7'°i" 

— /iclog/i 
t 

—  ft'logft 
e 

This  Thesis 
Chapter  5 

triangulations 

CHAPTER  3 


DYNAMIC  POINT  LOCATION:  CHAIN  METHOD 


3.1.  Introduction 

In  this  chapter  we  present  a  fully  dynamic  point  location  technique  for  monotone  subdivi¬ 
sions.  The  central  result  is  expressed  by  the  following  theorem: 

Theorem  3.1  Let  SI  be  a  monotone  planar  subdivision  with  n  vertices.  There  exists  a  dynamic 
point  location  data  structure  with  0(n)  space  requirement  and  OQog*n)  query  time,  which 
allows  for  insertion/deletion  of  a  vertex  in  time  O(logn)  and  insertion/deletion  of  a  chain  of  k 
edges  in  time  O  (log 2n  +k),  all  time  bounds  being  worst-case. 

It  must  be  underscored  that  our  method  allows  for  arbitrary  insertions  and  deletions  of  ver¬ 
tices  and  edges,  the  only  condition  being  that  monotonicity  of  the  subdivision  be  preserved.  It  is 
also  interesting  to  observe  that  the  dynamic  technique  is  based  on  the  same  geometric  objects, 
the  separating  chains,  which  yielded  the  first  practical,  albeit  suboptimal,  point  location  tech¬ 
nique  of  Lee-Preparata  [35],  and  later  on  the  practical  and  optimal  algorithm  of  Edelsbrunner- 
Guibas-Stolfi  [13] 

Our  technique  represents  a  reasonably  efficient  solution  of  the  dynamic  point  location  prob¬ 
lem  in  monotone  subdivisions.  It  requires  no  new  sophisticated  or  bizarre  data  structures,  and  it 
appears  eminently  practical. 

It  remains  an  open  problem  whether  O(logn)  optimal  performance  is  achievable  for 
query/update  times;  in  particular,  whether  the  technique  of  fractional  cascading,  which  achieved 
optimality  for  its  suboptimal  static  predecessor  [13],  can  also  be  successfully  applied  to  the  tech¬ 
nique  discussed  in  this  chapter. 

The  rest  of  the  chapter  is  organized  as  follows.  In  Section  3.2  we  review  the  technical 
background  and  formulate  the  problem.  In  Section  3.3  we  introduce  a  total  ordering  of  the 


regions  of  SR,  which  plays  a  crucial  role  in  the  dynamic  point  location  algorithm  illustrated  in 
Section  3.4. 

3.2.  Preliminaries 

We  recall  from  the  preceding  chapter  that  a  monotone  subdivision  is  a  partition  of  the  plane 
into  monotone  polygons  (see  Fig.  3.1(a)).  Given  a  monotone  subdivision  SR,  a  separator  a  of  SR 
is  a  monotone  chain  (vj,  •  •  ■  ,vp)  of  91  with  the  property  that  Vj  and  vp  are  points  at  infinity 
(hence,  each  horizontal  line  intersects  a  separator  either  in  a  point  or  in  a  segment).  A  separator 
of  SR  is  illustrated  with  bold  line  segments  in  Fig.  3.1(a).  Given  separators  Oi  and  02  of  91,  we 
say  that  is  to  the  left  of  a 2,  denoted  Oj  <02,  if,  for  any  horizontal  line  l  intersecting  both  0j 
and  02  ^  a  single  point,  the  abscissa  of  the  intersection  of  /  with  02  is  no  smaller  than  the 
corresponding  one  with  0^  A  partial  subdivision  is  the  portion  of  a  monotone  subdivision  con¬ 
tained  between  two  distinct  separators  0i  and  02,  with  0i  <  02.  A  complete  family  of  separators 


Figure  3.1  (a))  A  monotone  subdivision  SR.  The  dashed 
circle  represents  the  line  at  infinity.  A  separator  is  shown 
with  bold  line  segments,  (b)  The  regular  subdivision  SR* 
obtained  from  SR  by  forming  all  maximal  clusters. 


I  for  91  is  a  sequence  (01,02,  •  •  •  ,ot)  with  Oi  <  02  <  •  •  •  <  o„  such  that  every  edge  of  91  is 
contained  in  at  least  one  separator  of  £.  Notice  that  from  Euler’s  formula  t-0 (n).  As  shown  in 
[35],  every  monotone  subdivision  admits  a  complete  family  of  separators. 

Given  a  complete  family  of  separators  £  for  SR,  it  is  well-known  [35]  how  to  use  £  to  per¬ 
form  planar  point  location  in  91.  If  n  is  the  number  of  vertices  of  91,  then  in  time  O  (logn)  we 
can  decide  on  which  side  of  a  separator  the  query  point  q  lies;  applying  this  operation  as  a  primi¬ 
tive,  a  bisection  search  on  £  determines  in  time  0  (log t-  logn)  two  consecutive  separators 
between  which  q  lies.  This  process  can  be  adapted  or  supplemented  to  determine  the  actual 
region  r  to  which  q  belongs. 

Since  £  is  used  in  a  binary  search  fashion,  each  separator  is  assigned  to  a  node  of  a  binary 
search  tree,  called  the  separator-tree  3 .  With  a  minor  abuse  of  language,  we  call  “node  0”  the 
node  of  3  to  which  a  has  been  assigned.  An  edge  e  of  91  belongs,  in  general,  to  a  nonempty 
interval  (0,-,0,  +  i,  •  •  •  ,0y)  of  separators.  Let  node  ck,  i<,k<>  j,  be  the  common  ancestor  of 
nodes  0,-,0j+i,  •  •  •  ,ay,  then  e  is  called  a  proper  edge  of  a *  and  is  stored  only  once  at  node  0*. 
We  denote  by  proper(ak )  the  set  of  proper  edges  of  0*.  i.e.,  the  edges  of  0*  stored  at  node  0*. 
This  yields  0(n)  storage  space  while  guaranteeing  the  correctness  of  the  technique  (see 
[13, 35]).  Note  that  edges  whose  extremes  are  both  at  infinity  need  not  be  stored. 

We  now  illustrate  that  a  planar  subdivision  91  can  be  constructed  by  an  appropriate 
sequence  of  the  following  operations: 

INSERTPOINT (v,e;ei,e2):  Split  the  edge  e  =(u,w)  into  two  edges  e  1  =(u,v)  and 
e2  =  (v,w),  by  inserting  vertex  v. 

REMOVEPOJKT  (v;e):  Let  v  be  a  vertex  of  degree  2  whose  incident  edges,  e\  =(u,v)  and 
ej  =  (v,w),  are  on  the  same  straight  line.  Remove  v  and  replace  e  \  and  e^  with  edge  e  -  (u,w). 

INSERTCHAIN (y,V\,V2,r;r i,r2):  Add  the  monotone  chain  •  •  •  ,wk_1,v2), 

withy(v1)Sy(v2),  to  91  inside  region  r  of  91,  which  is  decomposed  into  regions  r  1  and  r2,  with 
r  1  and  r2,  respectively,  to  the  left  and  to  the  right  of  y,  directed  from  vx  tov2. 

REMOVECHAIN  (y ,r):  Let  y  be  a  monotone  chain  whose  nonextreme  vertices  have 
degree  2.  Remove  y  and  merge  the  regions  rx  and  r2  formerly  on  the  two  sides  of  y  into  region 
r,  [  The  operation  is  allowed  only  if  the  subdivision  91'  so  obtained  is  monotone.  ] 
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With  the  above  repertory  of  operations,  we  claim  that  a  monotone  subdivision  91  can 
always  be  transformed  into  a  monotone  subdivision  91'  having  either  fewer  vertices  or  fewer 
edges.  Then  by  0(n)  such  transformations  we  obtain  the  trivial  subdivision,  whose  only  region 
is  the  entire  plane  (bounded  by  the  line  at  infinity). 

Indeed,  let  a  be  a  separator  of  91,  and  imagine  traversing  a  from  -«©  to  +<».  Let  deg“(v) 
and  deg+(v),  respectively,  denote  the  numbers  of  edges  incident  on  v  and  lying  in  the  halfplanes 
y  < y(y)  and  y  >  y(v)  (assume  for  simplicity  that  no  edge  is  horizontal).  Let  (v  i ,  •  •  •  ,vp)  be  the 
sequence  of  finite  vertices  of  o  with  deg+(v,)  +  deg"(v,)  ;>  3,  as  encountered  in  the  traversal  of  o. 
If  this  sequence  is  empty,  the  entire  chain  can  be  trivially  removed.  Therefore,  assume  p  £  1.  If 
there  are  two  consecutive  vertices  v,-  and  v, + 1  such  that  deg+  (v,)  £  2  and  deg*  (v,  + 1 )  £  2,  then  the 
chain  y  (of  degree-2  vertices)  between  v(-  and  vi+1  can  be  deleted  by  REMOVECHAIN  while 
preserving  monotonicity.  Suppose  that  there  are  no  two  such  vertices.  If  deg“(v1)S2,  then  the 
portion  of  o  from  -«» to  v i  can  be  deleted.  Otherwise,  deg~(v  j)  =  1,  and  the  preceding  condi¬ 
tions  give  raise  to  the  following  chain  of  implications: 

[degTyx^lJ  s»  ^deg+(vj)£2j  =>  ^deg+(v2)^^  =>  •••  =>  ^deg+(v,)£^  ; 

the  latter  shows  that  the  portion  of  a  from  vp  to  +  «» can  be  deleted.  This  establishes  our  claim. 

When  all  finite  vertices  have  disappeared,  the  resulting  subdivision  consists  of  a  closed 
chain  of  arcs  whose  union  is  the  line  at  infinity.  Removal  of  the  vertices  at  infinity  completes  the 
transformation.  Since  all  of  the  above  operations  are  reversible,  this  shows  that  any  monotone 
subdivision  91  with  n  vertices  can  be  constructed  by  O(n)  operations  of  the  above  repertory: 

Theorem  3.2  An  arbitrary  planar  subdivision  91  with  n  vertices  can  be  assembled  starting  from 
the  empty  subdivision  by  a  sequence  of  0(n)  INSERTPOfNT  and  INSERTCHAIN  operations, 
and  can  be  disassembled  to  obtain  the  empty  subdivision  by  a  sequence  of  O  (n)  REMO- 
VEPOINT  and  REMOVECHAIN  operations. 

Although  the  above  operations  are  sufficient  to  assemble  and  disassemble  any  monotone 
subdivision,  the  following  operation  is  also  profitably  included  in  the  repertory: 


MOVEPOINT  (yjc,y):  Translate  a  degree-2  vertex  v  from  its  present  location  to  point 
(x ,y).  [  The  operation  is  allowed  only  if  the  subdivision  SR'  so  obtained  is  monotone  and  topo¬ 
logically  equivalent  to  SR.  ] 

33.  Ordering  the  Regions  of  a  Monotone  Subdivision 

Let  SR  be  a  monotone  subdivision,  and  assume  for  simplicity  that  none  of  its  finite  edges  is 
horizontal.  Given  two  regions  rx  and  r2  of  SR,  we  say  that  rx  is  left-adjacent  to  r2,  denoted 
r  i<^2,  if  r  i  and  r2  share  an  edge  e,  and  any  separator  of  SR  containing  e  leaves  rx  to  its  left  and 
r2  to  its  right.  Notice  that  relation  «r  is  trivially  antisymmetric.  But  we  can  also  show 

Lemma  3.1  Relation  c  on  the  regions  of  SR  is  acyclic. 

Proof:  Assume  r\<r2<  •  •  •  Let  I  be  a  complete  family  of  separators  and  let 

{<Ji,  •  •  •  cE  be  such  that  O/  separates  rx  and  rl  +  1,  so  that  ax  <a2  <  •  •  •  <o*_i.  If 

oe  Z  leaves  rk  to  its  left  and  rx  to  its  right,  we  have  o*_  j  <  a  and  a<  Oj,  a  contradiction  since 
the  separators  are  ordered.  □ 

Thus,  the  transitive  closure  of  c  is  a  partial  order,  referred  to  as  to  the  left  of,  and  denoted 
— Specifically,  rx-*r2  if  there  is  a  path  from  rx  to  r2  in  the  directed  graph  of  the  relation  « 
on  the  set  of  regions.  Correspondingly,  given  two  regions  rx  and  r2  of  SR,  we  say  that  rx  is 
below  r2,  denoted  r1fr2,  if  there  is  a  monotone  chain  from  HIGH  (r  x )  to  LOW (r2).  Obviously, 
T  is  a  partial  order  on  the  set  of  regions.  The  following  lemma  shows  that  these  two  partial  ord¬ 
ers  are  complementary. 

Lemma  3.2  Let  rx  and  r2  be  two  regions  of  SR.  Then  one  and  only  one  of  the  following  holds: 

ri-*r2,  r2-wi,  r2Trt. 

Proof:  Let  be  the  leftmost  separator  that  contains  the  left  chain  of  the  boundary  of  rx  and, 
analogously,  let  a*  be  the  rightmost  separator  containing  the  right  chain  of  the  boundary  of  rx. 
These  separators  partition  SR  into  five  portions,  each  a  partial  subdivision:  one  of  them  is  rx 
itself,  and  the  others  are  denoted  L,  R,  B,  and  T  (see  Fig.  3.2).  Now,  we  consider  four  mutually 
exclusive  cases  for  r2,  one  of  which  must  occur 


Figure  3.2  For  the  proof  of  Lemma  3.2. 

(1)  r2  e  L.  Consider  any  sequence  of  regions  (r^,  •  •  •  ,/s)  such  that  r2-*r'i,  »r',  +  1  for 

» =  1,  •  •  •  ,j  - 1,  and  the  right  chain  of  has  a  nonempty  intersection  with  aL.  If  the  right 
chain  of  r's  has  an  edge  that  is  also  on  the  left  boundary  of  rx,  then  r2->r  x .  Otherwise,  by 
the  definition  of  a*,,  there  is  a  sequence  /,+  i,  "  Sp  of  regions  such  that  /j-t/j+i  for 
j=s+ 1,  •  •  •  ,p-l,  and /p-w-!.  Thus,  in  aD  cases, r2-wi. 

(2)  r2  €  /?:  Arguing  as  in  (1),  we  establish  rx-¥r2. 

(3)  r2  €  B.  Since  LOW  (rj  is  the  highest  ordinate  vertex  in  B,  there  is  a  monotone  chain  from 
any  vertex  in  B  to  LOW (r^,  and,  in  particular,  from  HIGH  (r2)  to  LOW (r  1).  Thus  r2tr  i . 

(4)  r 2  e  T.  Arguing  as  in  (3),  we  establish  r1tr2.  □ 

We  say  that  rx  precedes  r2,  denoted  rx<r2,  if  either  rx  -*r2  or  rx  T r2. 

Theorem  3.3  The  relation  <  on  the  regions  of  51  is  a  total  order. 

As  an  example,  the  region  subscripts  in  Fig.  3.1(a)  reflect  the  order  < . 
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A  regular  subdivision  is  a  monotone  subdivision  having  no  pair  (rltr2)  of  regions  such  that 

riT r2.  For  example,  in  Fig.  3.1(a)  we  have  r9tr  10,  which  shows  that  the  illustrated  monotone 

\ 

subdivision  is  not  regular.  An  example  of  regular  subdivision  is  given  in  Fig.  3.3(a). 

The  significance  of  regular  subdivisions  is  expressed  by  Theorem  3.4  below.  It  is  easily 
realized  that  there  is  a  unique,  complete  family  Z  =  (Ci,  -  *  *  ,at)  of  separators  for  a  regular  subdi¬ 
vision  SR.  By  the  definition  of  separator,  all  regions  to  the  left  of  Oj  precede  all  those  to  its  right 
in  the  order  < .  Let  3  be  the  separator-tree  for  the  above  family  Z.  Recalling  the  rule  for  stor¬ 
ing  the  edges  of  separator  a  in  3 ,  as  reviewed  in  Section  3.2,  we  have 

Theorem  3.4  In  a  regular  subdivision  SR,  the  edges  of  proper  (a)  in  3  form  a  single  chain  (see 
Fig.  3.3(b)). 

Proof:  Assume  for  a  contradiction  that  a  contains  a  chain  y  which  is  the  bottom-to-top  concate¬ 
nation  of  three  nonempty  chains  yi ,  ft,  and  ft ,  where  Yi  and  ft  consist  of  proper  edges  of  o,  and 
ft  contains  no  such  edge.  Let  vj  and  v2  be  the  bottom  and  top  vertices,  respectively,  of  ft  and 
let  e'e  yx  and  e“  e  y2  be  the  edges  of  C  incident  on  Vi.  Since  e"  e  proper(o),  we  must  have 


(*)  0>) 


Figure  3.3  (a)  A  regular  subdivision,  and  (b)  its  chains  [proper (<r) :  as  Z). 


e"  €  proper  (o'),  where  node  o'  is  an  ancestor  of  node  o  in  3 .  We  claim  there  is  a  region  r  \ 
such  that  vx^HlGH(r\).  Otherwise,  each  separator  containing  e"  -  and  so  o'  -  also  contains 
«',  contrary  to  the  hypothesis  that  node  o  is  closest  to  the  root  among  the  nodes  whose  separator 
contains  e'.  Analogously,  we  show  that  there  is  a  region  r 2  for  which  v2  =LOW(r2).  Since  ft 
is  a  monotone  chain  from  HIGH  (r  x )  to  LOW  (r2),  then  whence  a  contradiction.  □ 

This  theorem  shows  that  a  regular  subdivision  has  a  particularly  simple  separator-tree.  In 
the  next  section  we  shall  show  that  the  property  expressed  by  Theorem  3.4  is  crucial  for  the 
efficient  dynamization  of  the  chain  method  for  point  location.  We  now  slightly  generalize  the 
notion  of  region  in  a  way  that  will  enable  us  to  show  that  any  monotone  subdivision  embeds  a 
unique  regular  subdivision.  We  say  that  two  regions,  rx  and  r 2,  consecutive  in  < ,  with  ritr2, 
are  vertically  consecutive.  We  then  have 

Lemma  3.3  If  rx  and  r2  are  two  vertically  consecutive  regions  of  a  monotone  subdivision  SR, 
then  the  chain  from  HIGH  (r  1 )  to  LOW (r2)  is  unique. 

Proof:  The  lemma  holds  trivially  if  HIGH  (r  £ )  =LOW (r2).  Thus,  assume  the  contrary.  Since 
r  1  Tr2,  there  is  a  monotone  chain  y  from  HIGH  (rj)  to  LOW  (r  2).  Suppose  now,  for  a  contradic¬ 
tion,  that  there  is  a  monotone  chain  /  from  HIGH  (r  1 )  to  LOW  (r2)  distinct  from  y.  Then  y u1/ 
defines  the  boundary  of  a  partial  subdivision  that  contains  at  least  one  region  of  SR.  For  any 
region  r  inside  this  partial  subdivision,  there  are  (possibly  empty)  chains  from  HIGH(rx)  to 
LOW(r)  and  from  HIGH(r)  to  LOW(r2),  so  that  n  trtr2,  contrary  to  the  hypothesis  that  rx 
and  r2  are  consecutive  in  < .  □ 

Given  two  vertically  consecutive  regions,  rx  and  r2,  in  SR,  with  rx^r2,  the  unique  chain 
from  HIGH  (r ! )  to  LOW  ( r2 )  is  called  a  channel. 

Lemma  3.4  All  channels  are  pairwise  vertex-disjoint 

Proof:  Assume,  for  a  contradiction,  that  there  are  two  channels  yj  and  ft  that  are  not  vertex- 
disjoint,  where  yi  connects  regions  rx  and  r2,  ft  connects  regions  r3  and  r4,  and 
rt  <r2<r3  <r4  (see  Fig.  3.4).  Since  yi  and  ft  share  a  vertex  x,  there  is  a  chain  from 
HIGH(rj)  to  LOW(r2),  which  consists  of  the  portion  of  ft  from  HIGH  (r 3)  tox,  and  the  portion 
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of  Yi  from  x  to  LOW (ri).  Hence,  we  have  r3  <  r 2,  a  contradiction. 


□ 


Given  two  vertically  consecutive  regions,  r\  and  with  riTrj,  we  imagine  duplicating 
the  channel  from  rj  to  r  2  and  view  the  measure- zero  region  delimited  by  the  two  replicas  as  a 
degenerate  polygon  joining  r  1  and  and  merging  them  into  a  new  region  ri  ur2  (see  Fig.  3.5). 
Gearly,  we  can  merge  in  this  fashion  any  sequence  of  vertically  consecutive  pairs.  This  is  for¬ 
mulated  in  the  following  definition: 

Clusters  are  recursively  defined  as  follows: 

(1)  An  individual  region  r  is  a  cluster, 

(2)  Given  two  vertically  consecutive  clusters  Xi  and  X2>  with  Xi  ?X2*  their  union  is  a  cluster  x, 
denoted  Xi-%2  (the  horizontal  bar  denotes  the  channel). 

A  maximal  cluster  x  is  one  which  is  not  properly  contained  in  any  other  cluster. 


Figure  3.4  Example  for  the  proof  of  Lemma  3.4. 


r 


Figure  3.5  Creation  of  a  channel  between  two  vertically  consecutive  regions. 

^  The  unique  subdivision  resulting  by  forming  all  maximal  clusters  of  91  is  denoted  91* .  Fig¬ 

ure  3.1(b)  illustrates  the  regular  subdivision  91*  corresponding  to  the  subdivision  91  of  Fig. 
3.1(a).  Notice  the  clusters  r2-r3,  r 6-r 7,  and  r^-r  10-r  n . 

®  The  above  definition  leads  us  to  a  convenient  string  notation  for  the  order  < ,  as  well  as  for 

the  cluster  structure  where  appropriate.  Normally,  we  shall  use  lower-case  roman  letters  for 
individual  regions,  lower-case  greek  letters  for  clusters,  and  upper-case  roman  letters  otherwise 
(i.e.,  for  collections  of  consecutive  regions  not  forming  a  single  cluster).  Specifically,  we  have 

(1)  A  cluster  (possibly,  a  region)  is  a  string. 

(2)  Given  two  strings  A  and  B,  such  that  the  rightmost  cluster  of  A  and  the  leftmost  cluster  of  B 

'  are  consecutive  (contiguity),  then  AS  is  a  string. 

A  subdivision  may  be  represented  by  means  of  its  structural  decomposition.  For  example,  the 
subdivision  of  Fig.  3.1(b)  is  described  by  the  string 


I 
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where  %\  =r2~ri,  Xi  =r6-r7,  and'x3=r9-rlQ-rn. 

Later  on,  we  will  find  it  convenient  to  explicitly  indicate  that  two  consecutive  regions  r  x 
and  r2  may  or  may  not  form  a  cluster.  We  shall  denote  this  with  the  string  notation  ri--r2, 
where  means  “potential  channel.” 

We  conclude  this  section  with  the  following  straightforward  observation: 

Theorem  3.5  The  subdivision  91*  obtained  by  forming  all  maximal  clusters  of  a  monotone  sub¬ 
division  91  is  regular. 

Note  that  in  the  transformation  of  91  to  91*  only  the  edges  of  channels  arc  duplicated.  By 
Lemma  3.4,  each  edge  is  duplicated  at  most  once,  thereby  ensuring  that  the  number  of  edges 
remains  O  (rt). 

3.4.  Dynamic  Point  Location 

3.4.1.  Data  structure  In  the  following  description,  we  assume  that  all  sorted  lists  are  stored  as 
red-black  trees  [20, 55].  We  recall  the  following  properties  of  red-black  trees,  which  are  impor¬ 
tant  in  the  subsequent  time  complexity  analyses. 

( 1 )  Only  0(1)  rotations  are  needed  to  rebalance  the  tree  after  an  insertion/deletion. 

(2)  The  data  structure  can  be  used  to  implement  concatenable  queues.  Operations  SPUCE  and 
SPLIT  of  concatenable  queues  take  O(logn)  time  and  need  O(logn)  rotations  each  for 
rebalancing. 

The  search  data  structure  consists  of  a  main  component,  called  the  augmented  separator 
tree,  and  an  auxiliary  component,  called  the  dictionary.  The  augmented  separator-tree  3  has  a 
primary  and  secondary  structure.  The  primary  structure  is  a  separator-tree  for  91* ,  i.e.,  each  of 
its  leaves  is  associated  with  a  region  of  91*  (a  maximal  cluster  of  91),  and  each  of  its  internal 
nodes  is  associated  with  a  separator  of  91* .  (The  left-to-right  order  of  the  leaves  of  the  primary 
structure  of  3  corresponds  to  the  order  <  on  the  regions  of  91*.)  The  secondary  structure  is  a 
collection  of  lists,  each  realized  as  a  search  tree.  Specifically,  node  a  points  to  the  list  proper  (a) 
sorted  from  bottom  to  top,  and  the  leaf  associated  with  cluster  %  (briefly  called  “leaf  x”)  points 
to  the  list  regions  (x)  of  the  regions  that  form  cluster  x,  also  sorted  from  bottom  to  top. 
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Given  two  regions  r\  and  consecutive  in  < ,  the  separator  o  between  r  i  and  is  asso¬ 
ciated  with  the  least  common  ancestor  of  the  leaves  associated  with  the  respective  clusters  of  rj 
and  ri  in  3 .  By  the  definition  of  separator-tree,  the  edges  of  c  are  stored  in  the  path  from  node 
a  to  the  root  of  3 ;  by  Theorem  3.4,  in  a  regular  subdivision,  each  extreme  vertex  of  proper{o) 
splits  proper(&),  for  some  ancestor  node  o'  of  node  o,  into  two  chains.  More  precisely,  the  fol¬ 
lowing  simple  lemma,  stated  without  ;,roof,  makes  explicit  the  allocation  of  the  edges  of  c  to  the 
nodes  of  3 . 

Lemma  3.5  Let  a  be  a  separator  of  91*,  and  Oi,  *  *  *  ,<Jk  be  the  sequence  of  nodes  of  3  on  the 
path  from  the  root  (=Oj)  to  o.  Then 

o=(alta2,  •  *  •  ,ak,proper(o),ph,pk-lt  •  •  •  .fo), 
where  a,  and  (5X-  are  (possibly  empty)  subchains  of  proper(a ;),  i  =  1,  •  •  •  ,h. 

In  order  to  dynamically  maintain  the  channels,  it  is  convenient  to  keep  two  representatives 
e'  and  e"  of  each  edge  e,  which  are  created  when  e  is  inserted  into  91.  If  e  does  not  belong  to  a 
channel,  e'  and  e"  are  joined  into  a  double  edge  and  belong  to  the  same  proper  (a).  If  instead  e 
is  part  of  a  channel,  then  e’  and  e"  are  single  edges  and  belong  to  distinct  proper  (o')  and 
proper  (a").  In  the  latter  case  e'  and  e"  are  on  the  boundary  of  the  same  cluster  x,  so  that  nodes 
o'  and  o"  are  on  the  path  from  leaf  x  to  the  root  of  3 .  Therefore,  we  represent  proper(o)  by 
means  of  two  lists,  denoted  strandl  (o)  and  strand2  (o).  List  strandl  (s),  called  primary  strand, 
stores  a  representative  for  each  edge  of  proper{o),  in  bottom-to-top  order.  List  strand2  (o), 
called  secondary  strand ,  stores  a  representative  for  each  double  edge  of  proper  (a),  in  bottom- 
to-top  order. 

Moreover,  associated  with  each  chain  proper  (a)  there  are  two  boolean  indicators  t(o)  and 
b(o),  corresponding  respectively  to  the  topmost  and  bottommost  vertices  of  proper(o). 
Specifically,  let  o'  be  the  ancestor  of  o  such  that  the  topmost  vertex  of  proper  (o)  is  an  internal 
vertex  of  the  chain  proper  (o')  (for  the  special  case  where  the  topmost  vertex  of  o  is  at  infinity, 
we  let  o'  be  the  father  of  o).  We  define  t(o)~left  if  o  is  to  the  left  of  o',  and  t(a)= right  if  o  is 
to  the  right  of  o'.  Parameter  b(a)  is  analogously  defined. 


The  dictionary  contains  the  sorted  lists  of  the  vertices,  edges,  and  regions  of  SR,  each  sorted 
according  to  the  alphabetic  order  of  their  names.  With  each  vertex  v  we  store  a  pointer  to  the 
representative  of  v  in  the  (at  most  two)  chains  proper  (a)  and  proper  (a)  of  which  v  is  a  nonex- 
treme  vertex.  With  each  edge  e  we  store  pointers  to  the  two  representatives  of  e  in  the  data 
structure.  Finally,  with  each  region  r  we  store  the  vertices  HIGH(r)  and  LOW  (r),  and  a  pointer 
to  the  representative  of  r  in  the  list  regions  (x)  such  that  %  is  die  maximal  cluster  containing  r. 
The  dynamic  maintenance  of  the  dictionary  in  the  various  operations  can  be  trivially  performed 
in  O  Gog  n)  time,  and  will  not  be  explicitly  mentioned  in  the  following. 

To  analyze  the  storage  used  by  the  data  structure,  we  note:  the  primary  structure  of  3  has 
0(n )  nodes,  since  there  are  O  (n)  regions  (by  Euler’s  formula)  and  therefore  O(n)  separators; 
the  secondary  structure  of  3  also  has  size  0(n),  since  there  are  O(n)  edges  in 
lproper(o) : a e  3 }  and  O(n)  regions  in  [regions  (y) :  Xe  3  }  (again,  by  Euler’s  formula);  the 
auxiliary  component  has  one  record  of  bounded  size  per  vertex,  edge,  and  region.  Therefore,  we 
conclude 

Theorem  3.6  The  data  structure  for  dynamic  point  location  has  storage  space  O  ( n ). 

Note  that  the  above  data  structure  is  essentially  identical  with  the  one  originally  proposed 
for  the  static  version  of  the  technique  [35].  What  is  remarkable  is  that  the  single-chain  structure 
of  the  proper  edges  of  any  given  separator,  due  to  our  specific  choice  of  the  separator  family,  is 
the  key  for  the  emergence  of  full  dynamic  capabilities. 

We  now  show  that  the  property  expressed  by  Theorem  3.4  allows  us  to  establish  an  impor¬ 
tant  dynamic  feature  of  the  data  structure.  According  to  standard  terminology,  a  rotation  at  node 
v  of  a  binary  search  tree  is  the  restructuring  of  the  subtree  rooted  at  v  so  that  one  of  the  children 
of  v  becomes  the  root  thereof.  A  rotation  is  either  left  or  right  depending  upon  whether  the  right 
or  left  child  becomes  the  new  root,  respectively.  We  then  have 

Lemma  3.6  A  rotation  at  a  node  of  3  can  be  performed  in  O  (log  n)  time. 
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Proof:  Without  loss  of  generality,  we  consider  a  left  rotation  as  illustrated  in  Fig.  3.6.  Clearly, 
the  separators  stored  at  nodes  outside  the  subtree  rooted  at  02  in  Fig.  3.6(a)  are  not  affected  by 
the  rotation,  nor  are  those  in  the  subtrees  rooted  at  <5\,  03,  and  03.  Thus  the  only  alterations 
involve  separators  02  and  04.  If  04  nproper(a 2)-0,  then  the  modification  is  trivial.  Suppose 
then  that  04  n proper (Oi)#0.  In  this  case  the  set  04  n  proper  {o 2)  forms  either  the  initial  or  the 
final  segment  of  the  chain  proper(a 2),  or  both,  for,  otherwise,  its  removal  from  proper(c2) 
would  yield  an  updated  (after  the  rotation)  proper (02)  violating  Theorem  3.4.  Thus  the  update 
is  accomplished  by  (1)  splitting  proper  (a2)  into  “fc  =  04  ^ proper  (a2)  and  its  relative  comple¬ 
ment  Yi ;  (2)  splicing  ft  with  proper (04)  to  form  the  updated  proper (04);  and  (3)  setting  the 
updated  proper (a2)  equal  to  Yi-  Note  that  the  extreme  vertices  of  proper (04)  are  obtained  in 
time  O  (1),  and  the  splitting  vertices  of  proper  (02)  are  determined  in  time  O  Gog  n).  Since  data 
structures  for  proper(a2)  and  proper (04),  (Le.,  the  red-black  trees  associated  with  lists 
strandl  (0)  and  strand 2  (0))  are  also  concatenate  queues,  the  splitting  and  splicing  operations 
are  executed  in  time  0(logn)  as  well.  The  parameters  r(0)  and  6(0)  for  the  resulting  separators 
are  updated  in  O  (1)  time  by  means  of  straightforward  rules.  □ 


Hereafter,  the  red-black  tree  3  is  assumed  to  be  balanced.  The  rest  of  this  section  is 
devoted  to  the  discussion  of  the  algorithms  to  perform  searches,  insertions,  and  deletions.  We 
have  shown  in  Section  3.2  that  the  four  operations  INSERTCHAIN,  REMOVECHAIN,  INSERT- 
POINT,  and  REMOVEPOINT  are  sufficient  to  generate  any  planar  subdivision.  The  measures  of 
time  complexity  ought  to  be  expressed  as  functions  of  the  form  f(n,k),  where  n  is  the  current 
size  of  31,  and  k  is  the  size  of  the  chain  y  to  be  inserted  or  deleted. 

3.4.2.  Query 

To  perform  a  point  location  search  for  a  query  point  q,  we  use  essentially  the  same  method 
as  [35].  The  search  consists  of  tracing  a  path  from  the  root  to  a  leaf  %  of  3 .  At  each  internal 
node  0  we  discriminate  q  against  separator  0.  Three  cases  may  occur 

(1)  q  €  0:  we  return  the  edge  of  0  that  contains  q  and  stop; 

(2)  q  is  to  the  left  of  0:  we  proceed  to  the  left  child  of  0; 


(3)  q  is  to  the  right  of  c:  we  proceed  to  the  right  child  of  o. 

Once  we  reach  a  leaf  x,  we  know  that  q  belongs  to  a  region  of  X-  Since  the  regions  of  x  are 
sorted  from  bottom  to  top,  such  region  is  determined  by  searching  in  the  list  regions  (x).  The 
above  technique  can  be  viewed  as  a  “horizontal”  binary  search  in  the  set  of  separators  of  91*, 
followed  by  a  “vertical”  binary  search  in  the  set  of  regions  of  the  leaf  %■ 

Let  e  be  the  edge  of  a  whose  vertical  span  contains  y(q).  When  e  e  proper  (a),  the 
discrimination  of  q  against  o  is  a  conventional  search  in  strandl  (a).  When  e  €  proper  (a)  then 
we  use  the  pair  (r(c),b(o)):  for  example,  when  e  is  above  proper  (a),  if  t(o)=left,  then  q  is 
discriminated  to  the  right  of  o,  and  to  its  left  otherwise.  (This  is  a  minor  variant  of  the  criterion 
adopted  in  [35]).  The  case  when  e  is  below  properic)  is  treated  analogously.  This  simple 
analysis  confirms  that  the  time  spent  at  each  node  is  O(logn).  We  have  [35] 

Theorem  3.7  The  time  complexity  of  the  query  operation  is  O  (log2 a). 

3.4.3.  Insertion 

We  shall  first  show  that  the  effect  of  operation  1NSERTCHAJN (y,vx,v2,r;rx,r2)  on  the 
order  <  of  the  regions  of  91  can  be  expressed  as  a  syntactical  transformation  between  the  strings 
expressing  the  order  before  and  after  the  update.  The  situation  is  illustrated  in  Fig.  3.7. 

On  the  boundary  of  r  there  are  two  distinguished  vertices:  HIGH  (r  x )  and  LOW (r2).  Note 
that  HIGH (r ! ) = HIGH (r)  if  v2  is  on  the  right  chain  of  the  boundary  of  r  (and  similarly 
LOW (r2) -LOW (r)  if  vj  is  on  the  left  chain).  Thus,  in  general,  HIGH(rx)  is  on  the  left  chain 
of  the  boundary  of  r,  and  LOW  (r2)  is  on  the  right  chain.  Using  the  string  notation  introduced  in 
Section  3.3,  let  L  and  R  be  the  strings  corresponding  to  the  regions  that  respectively  precede  and 
follow  r  in  < .  Thus,  the  subdivision  91*  is  described  by  the  string  LrR. 

Let  e  \  be  the  edge  of  91*  on  the  left  boundary  of  r  incident  on  HIGH(rx )  from  below,  and 
let  x  be  the  maximal  cluster  on  the  left  of  ex.  In  general,  this  cluster  consists  of  two  portions,  Xi 
and  Xz  (such  that  Xi“%2)>  where  X2  consists  exactly  of  the  regions  q'  of  x  for  which 
y(LOW (q'))'Zy(HIGH (jr x)).  Thus,  we  have  L  =L'%L”.  We  now  distinguish  three  cases  and 
define  substrings  Xi ,  X2,  L  x ,  and  L2  as  follows. 


Figure  3.7  (a)  Canonical  partition  of  subdivision  9t*  with 
reference  to  region  r  and  vertices  vj  and  vj,  and  (b)  the 
restructured  subdivision  after  the  insertion  of  chain  y 
between  and  V2- 


(1)  X2*&-  Lct‘Ki=X\,‘^2=X2fLi^L\L2=L”,so\hatL=Li\i-'K2L2. 

(2)  X2  =  0-  Let  q  be  the  region  preceding  r  (note  q  could  form  a  cluster  with  r).  We  further 
distinguish: 

(2.1) y  (LOW  (q))~Zy  (HIGH  (r \)).  In  this  case  we  let  L  =L  1X2^2.  where  X2  is  die  maxi¬ 
mal  cluster  immediately  following  %. 

(2.2)  y  (LOW  (q))  <y  (HIGH  (r  1)).  In  this  case  we  let  L=LiX\--,  where  X!  is  the  right¬ 
most  maximal  cluster  of  L  (but  not  necessarily  a  maximal  cluster  in  91*). 

The  three  cases  are  conveniently  summarized  by  the  notation 

L  —L\X\  -  -K2L2. 

Analogously,  string  R  can  be  reformulated  as 

R  =R  iPi-*P2#2 

with  straightforward  meanings  of  the  symbols. 
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Thus,  in  general,  for  any  given  region  r  and  choice  of  vj  and  v2  on  its  boundary,  we  have 
the  following  canonical  string  decomposition  of  9?*: 

LjXj  -  lPl  ‘  "P2-^2- 

The  corresponding  partition  of  the  subdivision  is  illustrated  in  Fig.  3.7(a).  Examples  of 
configurations  corresponding  to  cases  (2.1)  and  (2.2)  are  shown  in  Fig.  3.8.  Namely,  part  (a) 
shows  case  (2.1)  for  L  and  case  (1)  for  R,  part  (b)  shows  case  (2.2)  for  both  L  and  R,  and  pan  (c) 
shows  case  (2.2)  for  L  and  case  (1)  for  R. 

We  now  investigate  the  rearrangement  of  this  order  caused  by  the  insertion  of  chain  y  into 
r.  Referring  to  Fig.  3.7(b),  it  is  immediately  observed  that  the  order  after  the  update  is  as  fol¬ 
lows: 

L1<X1<TJ<RJ<PJ  <X2<L2<T2<p2</?2- 

To  obtain  the  string  description  of  the  updated  subdivision  we  must  determine  whether  any  new 
channel  has  been  created.  Any  such  channel  can  only  arise  in  correspondence  with  a  new  adja¬ 
cency  caused  by  the  update,  specifically  for  the  following  pairs:  (Xj.rj),  (pi,X2),  and  (r2,P2)- 
The  channel  from  Xj  to  r  j  exists  only  if  y  (HIGH  (Xj ))  £y(LOW  {r  j )),  and  analogously  for  the 
channel  from  r2  to  p2.  Instead,  since  y(HIGH (pj)) <y(LOW (X2)),  the  cluster  pj-X2  always 
exists.  Therefore,  the  order  caused  by  the  insertion  of  yis  represented  by  the  string 

LjXj-  -TjR  jPj— X2L2T2*  -p2R2- 

(In  purely  syntactic  terms,  this  transformation  corresponds  to  rewriting  r  as  r2-r  j  and  then 
exchanging  substrings  rj/?  j  pj  and  X2L2r2.)  This  is  summarized  as  follows: 

Theorem  3.8  Let  L  j  Xj  -  -X2 L  2r/?  j  pj  -  -p2R  2  be  the  string  description  of  the  order  of  91* ,  where 
Lj,  r,  and  R2  are  nonempty.  After  operation  INSERTCHAIN  (y,vi,v2,r;ri,r2)  the  new  order  is 
described  by 

L  j  Xj  -  -r  j/?ipi-X2L2r2*-p2^2- 

The  algorithm  for  the  INSERTCHAIN  (y,vi,v2,T;ri,/,2)  operation  implements  the  syntacti¬ 
cal  transformation  of  the  string  description  by  decomposing  the  subdivision  91  into  its  com¬ 
ponents  L j ,  Xj ,  X2,  L2,  t,  R  j ,  pj ,  P2,  and  R2>  which  are  subsequently  reassembled  according  to 
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the  new  order  given  by  Theorem  3.8. 

To  formally  describe  the  algorithm,  we  denote  by  91(5)  the  partial  subdivision  associated 
with  a  string  5  of  consecutive  regions  of  SR.  We  can  represent  91(5)  with  the  same  data  structure 
described  in  Section  3.4.1,  and  we  denote  with  3  (5)  the  augmented  separator-tree  for  5.  Partial 
subdivisions  can  be  cut  and  merged  with  the  same  rules  as  for  the  decomposition  and  concatena¬ 
tion  of  the  corresponding  strings. 

Let  91(5 1),  91(52),  and  91(5)  be  partial  subdivisions  such  that  5  =5j52-  We  show  in  the 
following  how  to  merge  3  (5 1 )  and  3  (52)  into  3  (5),  and  how  to  cur  3  (5)  to  produce  3  (5 1 )  and 
3(52).  The  merge  operation  needs  also  the  separator  a  at  the  boundary  between  91(5 1)  and 
91(52),  represented  by  its  primary  and  secondary  strands.  The  cut  operation  returns  the  separator 
a.  These  operations  can  be  implemented  by  means  of  the  following  six  primitives: 

Procedure  MERGE1  (5 1 ,0,52:5) 

[  It  merges  partial  subdivisions  91(5 1)  and  91(52),  with  5i— >52',  o  is  the  separator  between 
91(50  and  9?(52).  ] 

(1)  Construct  a  separator-tree  3  (5)  for  91(5),  by  placing  a  at  the  root,  and  making  3  (5 1 )  and 
3  (5  2)  the  left  and  right  subtrees  of  O,  respectively. 

[  3  (5)  is  a  legal  separator-tree  for  91(5),  but  might  be  unbalanced.  ] 

(2)  Rebalance  3  (5)  by  means  of  rotations. 

Procedure  MERGE2(Xi>ct,x2,X ) 

[  It  merges  partial  subdivisions  91(Xi)  and  9?(x2)  such  that  Xi^X2  into  9t(x),  where  X=Xi~X2, 
and  a  is  the  channel  between  Yi  and  X2-  ] 

(1)  Separate  the  two  strands  of  a,  and  make  the  secondary  strand  become  a  new  primary 
strand. 

(2)  Splice  regions (Xi)  and  regions  (xi)  to  form  regions  (x)- 


Lemma 3.7  Operations  MERGE1(S  i,a,S2‘,S)  and  Af£/?G£2(Xi  ,a,X2'0C)  have  time  complexity 
0  (log2;*)  and  0  (log  n),  respectively. 

Proof:  The  time  bound  for  operation  MERGE2  follows  immediately  from  the  properties  of  con- 
catenable  queues.  With  regard  to  MERGE1,  Step  1  consists  of  forming  3  (S)  by  joining  the  pri¬ 
mary  structures  of  3  (S  i)  and  3  (S 2)  through  node  o,  which  takes  0(1)  time.  Since  we  use  red- 
black  trees,  we  can  rebalance  3(5)  with  O(logn)  rotations  [55,  pp.  52-53].  By  Lemma  3.6, 
each  such  rotation  takes  O  (log  n )  time,  so  that  the  total  time  complexity  is  O  (log2n).  □ 


Procedure  CUT1  (S,Xi,X2',Si,g,S2) 

[  It  cuts  partial  subdivision  9t(S)  into  9t(S  1)  with  rightmost  cluster  Xi  and  9t(S 2)  with  leftmost 

cluster  %2>  such  that  Xi~*X2»  and  also  returns  the  separator  o  between  SR(Si)  and  91(52).  ] 

(1)  Find  the  node  o  of  3  (5)  that  is  the  lowest  common  ancestor  of  leaves  Xi  and  %2- 

(2)  Perform  a  sequence  of  rotations  to  bring  O  to  the  root  of  3  (S),  where  after  each  rotation  we 
rebalance  the  subtree  of  a  involved  in  the  rotation,  namely,  the  left  subtree  for  a  left  rota¬ 
tion  and  the  right  subtree  for  a  right  rotation  (see  Fig.  3.9). 

(3)  Set  3  (5 1)  as  the  left  subtree  of  a  and  3  (52)  as  the  right  subtree  of  a.  Return  the  chain 
proper  (a). 

Procedure  CUT2(x>X\&Xi) 

[  It  cuts  partial  subdivision  9l(x)  into  9l(Xi)  and  9l(X2)*  where  X=Xi~%2  and  a  is  the  channel 

between  Xi  and  X2- 1 

(1)  Join  the  two  previously  separated  strands  of  a,  so  that  the  rightmost  one  becomes  the 
secondary  strand  of  the  other. 

(2)  Split  regions (x)  into  regions(%  1 )  and  regionsix 2). 

Lemma 3.8  Operations  CUT1  (5,Xi,X2;*Si,c,52)  a™*  CUT2 (x;Xi>Y»X2)  have  time  complexity 

0(log2n)  and  O  (log  n),  respectively. 
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Proof:  The  time  bound  for  operation  CUT2  is  immediate.  With  regard  to  operation  CUT1,  Step 
1  takes  O  (height (Si  (S)))  =  0 (log  n )  time.  In  Step  2,  we  perform  no  more  than  height (3  (S)) 
rotations  to  bring  a  to  the  root  After  each  such  rotation,  we  have  to  rebalance  a  subtree  3 ' 
whose  left  and  right  subtrees  are  already  balanced,  so  that  the  number  of  rotations  required  for 
rebalancing  is  proportional  to  the  difference  of  height  of  the  subtrees  of  3 '.  Such  differences 
form  a  sequence  whose  sum  is  proportional  to  height (3  (S))  [55,  p.  53].  To  see  this,  consider  the 
path  of  3  (S)  from  node  o  to  the  root  and  let  X  i ,  •  •  •  ,Xp  be  the  sequence  of  left  subtrees  of  the 
nodes  of  this  path  whose  right  child  also  belongs  to  the  path,  in  bottom-to-top  order.  Also, 
define  X\  =Xi.  The  t-th  left-rotation  performed  in  bringing  o  to  the  root  splices  trees  Xi+\  and 
X'i  to  form  a  new  tree,  which  after  rebalancing  will  be  denoted  X',  + 1 .  We  recall  that  each  node 
p  of  a  red-black  tree  is  associated  with  an  integer  rank(ii)  such  that 

rank(n)  <,  height  ([i)  £  2  rank(\i). 

Also,  the  rank  of  a  subtree  is  equal  to  the  rank  of  its  tool  From  the  theory  of  red-black  trees  [55, 
p.  53]  we  have 

(1)  rank (X, )  £  rank  (X,-  + 1 ); 

(2)  rank(Xi)  £  rank(X'i)  £  rank(Xt )  + 1; 

(3)  the  number  r,  of  rotations  needed  to  obtain  the  balanced  tree  X',  +  1  is  at  most 
r,  I  rank(Xi + 1 )  ~  rank  (X'i)  I . 

The  total  number  t  of  rebalancing  rotations  is 

t  =P£  n  =P£  I rank(Xi +l)- rank(X'i)  I . 

i  =  1  i  =  1 

Now,  if  rank(X'i)  =  rank(Xi)  or  rankQii)  <  rank(Xi+l),  we  have 

I  rank(Xi +i)~  rank(X'i)  I  =  rank(X{  +1)- rank  (X'i)  £  rank(Xt +1)- rank(X{). 

Otherwise  (i.e.,  rank  (X'i)  =  rank(Xi)  + 1  and  rank(Xi)  =  rankQii  + 1 )),  we  have 
I rank(Xi+i)~rank(X'i)\  =  \rank(Xi+i)-rank(Xi)-l  I  =1. 

Hence,  in  both  cases  the  following  inequality  holds: 


32 


D 


r. 


r 


■ 


I 


\rank(Xi+i)-rank(X'i)\  £rank(Xi+l)- rank(Xi)+ 1, 

from  which  we  obtain 

(rank(Xi+1)-rank(Xi)+ 1  )Zrank(Xp)+(p  -  l)Z2height&  (S))  =  0(log/i). 
i=  1 

By  Lemma  3.6,  each  rotation  takes  0  (log n)  time,  so  that  the  total  time  complexity  is  O  (log2/!). 

□ 

Procedure  FINDLEFT  (c;%) 

[  It  finds  the  cluster  x  to  the  left  of  edge  e.  If  e  is  part  of  a  channel,  then  %  is  the  cluster  that  con¬ 
tains  such  a  channel.  ] 

(1)  Perform  a  point  location  search  for  (any  point  of)  edge  e.  The  search  will  stop  at  a  node  a 
of  3  that  stores  (a  representative  of)  e. 

(2)  If  e  is  a  double  edge  of  c  [  i.e.,  e  does  not  belong  to  a  channel  ],  resume  the  point  location 
search  in  the  left  subtree  of  o  and  return  the  leaf  x  where  the  search  terminates.  [  This 
corresponds  to  searching  for  a  point  p'  immediately  to  the  left  of  edge  e.  ] 

(3)  Otherwise  [  i.e.,  e  is  a  single  edge  of  o  and  belongs  to  a  channel  ]  resume  the  point  location 
search  in  both  subtrees  of  o.  [  This  corresponds  to  searching  for  points  p~  and p+  immedi¬ 
ately  to  the  left  and  right  of  edge  e,  respectively.  ]  One  of  the  two  searches,  say  the  left  one, 
will  terminate  in  a  leaf,  while  the  other  search,  say  the  right  one,  will  stop  at  a  node  o'  that 
stores  the  other  representative  of  e.  [  Recall  that  the  two  nodes  storing  e  are  on  the  path 
from  leaf  x  to  the  root  of  3 .  ]  We  resume  the  search  in  the  left  subtree  of  o'  and  return  the 
leaf  x  where  the  search  terminates.  [  The  case  where  the  right  search  out  of  a  terminates  in 
a  leaf  is  analogous.  ] 

Procedure  FINDRIGHT  (e  ;x) 

[  It  finds  the  cluster  x  to  the  right  of  edge  e.  If  e  is  part  of  a  channel,  then  x  is  the  cluster  that 

contains  such  a  channel  ] 

( Analogous  to  FINDLEFT  ] 
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Lemma  3.9  Operations  FINDLEFT  (e;%)  and  FIND  RIGHT  (e ;%)  have  each  time  complexity 
0  (log^n). 

Proof:  Since  each  edge  has  two  representatives,  there  are  at  most  two  nodes  of  3  where  we 
proceed  to  both  children.  Hence,  we  visit  a  total  of  0(log/i)  nodes,  spending  OQogn)  time  at 
each  node.  C3 

The  complete  algorithm  for  operation  INSERTCHAIN  (y,v  i  ,v2,r;r  t  ,r2)  is  as  follows: 
Algorithm  INSERTCHAIN  (Y,vi,v2,r;r1,r2) 

(1)  Find  regions  q  and  s  immediately  preceding  and  following  r,  respectively;  also,  find  clus¬ 
ters  Xl  and  Xr  by  means  of  FINDLEFT (e  \  and  FINDRJGHT (c2;X/?)-  From  these 
obtain  ^ ,  X2,  pi ,  and  p2. 

(2)  Perform  a  sequence  of  CUT1  and  CUT2  operations  to  decompose  SR=SR(LiX1-- 

L2L2r/?1p1--p2fl2)  into  9i(Z,i),  $0*),  9J(L2),  9?(r),  WO,  9l(pi).  9l(p2),  and 

9t(/?2).  The  primary  and  secondary  strands  returned  by  each  such  operation,  which  form 
the  boundaries  of  the  above  partial  subdivisions,  are  collected  into  a  list  p. 

(3)  Construct  the  primary  and  secondary  strands  of  chain  y  and  add  them  to  p. 

(4)  Destroy  9?(r)  and  create  5K(r  j )  and  SR(r2). 

(5)  Test  for  channels  Xi-fT  and  r2-p2,  and  perform  a  sequence  of  MERGE1  and  MERGE2 
operations  to  construct  the  updated  subdivision  9?(LjX.i--rj/?iPi-X.2L2r2--p2/?2).  The 
separators  and  channels  needed  to  perform  each  such  merge  are  obtained  by  splitting  and 
splicing  the  appropriate  strands  of  p. 

Theorem  3.9  The  time  complexity  of  operation  INSERTCHAIN  (Y.vi  ,v2,r;r!,r2),  where  Ycon- 
sists  of  k  edges,  is  0(log2n  +k). 

Proof:  In  Step  1,  finding  q  and  s  takes  0  (log  n)  time.  In  fact,  q  is  either  in  the  cluster  of  r  or  in 
the  cluster  immediately  preceding  the  one  of  r,  and  analogously  for  s.  By  Lemma  3.9,  finding 
Xl  and  xr  takes  O (log2/:)  time.  The  remaining  computation  of  X2,  pi,  and  p2  can  be  done 
in  0(logn)  time.  By  Lemma  3.8,  Step  2  takes  0(log2n)  time.  Notice  that  the  list  p  has  0(1) 
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elements.  Step  3  can  be  clearly  performed  in  time  0(k).  Step  4  takes  O  (1)  time  since  r,  r  ^  and 
r2  are  single-region  structures.  In  Step  5,  testing  for  channels  and  r2-  -p2  takes  0(1) 

time.  Since  the  list  p  has  O  (1)  elements,  we  can  construct  in  O  (logn)  time  the  separators  and 
channels  needed  for  each  merge  operation  of  Step  5.  By  Lemma  3.7  the  total  time  for  such 
merges  is  O  (log2/!).  □ 

With  regard  to  the  INSERTPOINT  operation,  we  locate  the  edge  e  in  the  dictionary,  and 
replace  each  of  the  two  representatives  of  e  in  the  data  structure  with  the  chain  (e  j,v,e2).  This 
corresponds  to  performing  two  insertions  into  sorted  lists,  so  that  we  have 

Theorem  3.10  The  time  complexity  of  operation  INSERTPOINT  (v,e;e  j  ,e2)  is  O(logn). 

A  similar  argument  shows  that 

Theorem  3.11  The  time  complexity  of  operation  MOVEPOINT  iy\x,y)  is  <9(log/i). 

3.4.4.  Deletion 

The  transformations  involved  in  a  REMOVECHAIN  operation  are  exactly  the  reverse  of  the 
ones  for  the  INSERTCHAIN  operation.  We  observe  that  all  the  updates  performed  in  the  latter 
case  are  totally  reversible,  which  establishes 

Theorem  3.12  The  time  complexity  of  operation  REMOVECHAIN  (y ,r),  where  y  consists  of  k 
edges,  is  O  (log2  n+k). 

The  same  situation  arises  with  respect  to  the  INSERTPOINT  and  REMOVEPOINT  opera¬ 
tions,  so  that  we  have 

Theorem  3.13  The  time  complexity  of  operation  REMOVEPOINT  (v;e)  is  <9 (log  n). 

Theorem  3.1  stated  in  Section  3.1  results  from  the  combination  of  the  above  Theorems  3.6, 


CHAPTER  4 


DYNAMIC  POINT  LOCATION:  TRAPEZOID  METHOD 


4.1.  Introduction 

In  this  chapter  we  present  a  dynamic  point  location  technique  with  optimal  query  time. 
This  technique  considers  convex  subdivisions  whose  vertices  lie  on  a  fixed  set  of  N  horizontal 
lines.  The  supported  update  operations  are  insertion/deletion  of  vertices  and  edges,  and  (hor¬ 
izontal)  translation  of  vertices.  The  method  achieves  query  time  0(\ogn+\ogN),  space 
0(n  log N),  and  insertion/deletion  time  O  (log  a  logN).  Hence,  whenever  N  =  0(n),  the  query 
time  is  0(logn),  which  is  optimal.  In  addition  to  its  good  theoretical  performance,  this  tech¬ 
nique  is  also  easy  to  implement  and  very  efficient  in  practice.  It  is  based  on  the  trapezoid 
method  of  Prepara ta  [46],  which  has  been  experimentally  shown  to  be  the  fastest  point  location 
method  among  those  with  asymptotically  optimal  query  time  [12].  It  is  easily  realized  that  in 
many  significant  applications  the  most  frequent  operation  is  the  point  location  query,  while 
updates  are  more  rarely  executed.  Hence,  this  technique  provides  the  most  efficient  solution  for 
such  applications. 

The  rest  of  the  chapter  is  subdivided  into  seven  sections.  In  Section  4.2  we  formally  state 
the  problem.  Section  4.3  describes  the  geometric  foundations  of  the  technique.  Sections  4.4  and 
4.5  present  the  dynamic  data  structure  and  the  query  algorithm.  The  update  algorithms  are  dis¬ 
cussed  in  Sections  4.6  and  4.7.  Finally,  Section  4.8  discusses  extensions  of  the  technique  and 
directions  of  further  research. 

4.2.  Preliminaries 

Let  L  =  {/0,/i,  ■  •  •  ,//y)  be  a  set  of  horizontal  lines,  in  this  order  from  bottom  to  top.  The 
lines  of  L  partition  the  plane  into  horizontal  strips,  called  elementary  slabs,  two  of  which  are  (the 
bottom  and  top)  half-planes.  A  slab  is  either  an  elementary  slab  or  the  union  of  two  contiguous 
slabs.  We  consider  a  convex  subdivision  91  whose  n  vertices  are  on  the  lines  of  L  (see  Fig.  4.1). 
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This  implies  that  all  finite  edges  of  SR  occur  in  the  plane  strip  between  / o  and  fa,  and  that  all 
infinite  edges  are  horizontal  rays. 

The  following  update  operations  on  SR  are  defined: 

INSERTPOINT  (v,e;e !  ,e2): 

Split  the  edge  e  =  (u,w)  into  two  edges  ei=(u,v)  and  e2  =  (v,w),  by  inserting  vertex  v, 
which  lies  on  some  line  of  L. 

REMOVEPOINT(v,eue2;e ): 

Let  v  be  a  vertex  of  degree  2  whose  incident  edges,  e  \  =  (u,v)  and  e2=(v,w),  are  on  the 
same  straight  line.  Remove  v  and  replace  e  t  and  e  2  with  edge  e  =  (u,w). 

INSERTSEGMENT  (e,v  x,v2trtrx  ,r  2): 

Add  the  edge  e  =(v!,v2),  with  y(v!)£y(v2),  to  SR  inside  region  r  of  31,  which  is  decom¬ 
posed  into  regions  rx  and  r2,  with  rx  and  r2,  respectively,  to  the  left  and  to  the  right  of  e, 
directed  from  v x  to  v2. 


Figure  4.1  Example  of  convex  subdivision  with  vertices 
on  the  set  L  of  horizontal  lines. 


REMOVESEGMENT  (e,vx,v2,r  i  ,r2  ;r): 

Remove  edge  e  and  merge  the  regions  rx  and  r2  formerly  on  the  two  sides  of  e  into  region 
r.  [  The  operation  is  allowed  only  if  the  subdivision  51'  so  obtained  is  convex.  ] 

MOVEPOINT  (yjt): 

Horizontally  translate  a  vertex  v  from  its  present  location  to  abscissa  x.  [  The  operation  is 
allowed  only  if  the  subdivision  9?'  so  obtained  is  convex  and  topologically  equivalent  to 
*.] 

The  following  theorem  demonstrates  the  adequacy  of  the  above  repertory: 

Theorem  4.1  An  arbitrary  convex  subdivision  91  with  n  vertices  can  be  assembled  starting  from 
the  empty  subdivision,  and  disassembled  to  obtain  the  empty  subdivision,  by  a  sequence  of  O  (n) 
operations  from  the  above  repertory. 

Proof:  It  suffices  to  constructively  demonstrate  the  assembly  of  91,  since  its  disassembly  is 
obtainable  by  replacing  each  operation  with  its  inverse  and  by  reversing  their  order.  First,  we 
recall  the  definition  of  the  acyclic  relation  between  the  regions  of  91  (see  Section  3.3  of 
Chapter  3).  Given  rx  and  r2  in  91,  r  j  is  left-adjacent  to  r2  (rj*cr2)  if  their  boundaries  share  an 
edge  e,  and  rx  and  r2  are  respectively  to  the  left  and  to  the  right  of  e.  (If  rx  and  r2  share  a  hor¬ 
izontal  edge,  then  r!<r2  if  rx  is  below  e .)  Consider  an  order  (r1,r2,  •  •  •  ,/y)  of  the  regions  of 
91  which  is  consistent  with  We  build  91  by  adding  its  regions  one  by  one,  according  to  the 
chosen  order.  At  the  generic  step  of  the  construction,  after  adding  region  r,  we  call  91,  the 
resulting  subdivision.  (Obviously,  9ty=91.)  In  91,  there  is  a  monotone  polygonal  line  a„  which 
forms  the  boundary  of  the  set  of  regions  (r1,r2,  •  •  •  ,r, }  (refer  to  Fig.  4.2).  Consider  the  (exter¬ 
nal)  angles  to  the  right  of  o,-  at  each  nonextreme  vertex  of  a,-:  each  angle  >  n  is  bisected  by  a 
horizontal  ray  to  ensure  the  property  of  convexity  for  91;.  Consider  now  region  r,+1:  obviously, 
the  left  chain  of  r,+i  is  already  present  in  9t„  and  let  the  right  chain  of  r,  +  1  be  v'x,  •  •  •  ,v',.  By 
means  of  INSERTSEGMENT  we  insert  the  edge  e  =  ( LOW  (r;  +  j  ),HIGH  (r,- + x )).  (This  operation 
is  performed  as  many  times  as  there  are  regions.)  Next,  with  the  sequence  INSERTPOINT, 
INSERTSEGMENT,  MOVEPOINT  we  insert  vertex  v'2  and  edge  (v'j,v  2)  into  their  final  posi¬ 
tion,  and  partition  the  region  to  the  right  by  a  horizontal  ray  to  preserve  convexity.  This  process 


is  performed  for  v'3,  •  •  •  ,v's~\.  (This  sequence  of  operations  is  performed  at  most  once  for 
each  vertex  of  91.)  Finally,  by  REMOVESEGMENT  we  remove  any  horizontal  ray  originating  at 
v'i  and  v's  if  the  corresponding  external  angles  have  become  Ire.  (This  operation  is  performed 
at  most  twice  per  region.)  The  resulting  subdivision  is  +  y  and  exhibits  the  same  invariant  as 
regards  the  region  to  the  right  of  o,+ 1.  The  comments  about  the  frequency  of  each  executed 
operation  complete  the  proof.  □ 

Let  r  be  a  convex  polygon.  Polygon  r  is  said  to  be  1 -supported  if  it  has  an  edge  e  such  that 
any  horizontal  line  intersecting  r  intersects  also  e,  i.e.,  r  contains  the  edge  (LOW  ( r),HlGH  (r)) 
(see  Fig.  4.3(a)).  Polygon  r  is  said  2 -supported  if  it  has  two  edges,  e\  and  e2.  such  that  any  hor¬ 
izontal  line  intersecting  r  intersects  at  least  one  of  e  1  and  e2  (see  Fig.  4.3(b)).  (Notice  that  a  1- 
supported  polygon  is  also  2-supported.)  /»  convex  subdivision  is  said  to  be  1 -supported  (respec¬ 
tively  2-supported)  if  all  its  regions  are  1-supported  (respectively  2-supported). 

An  important  property  of  i-supportedness  (i  =  1,2)  is  expressed  by  the  following  straightfor¬ 
ward  lemma: 

Lemma  4.1  Let  convex  polygon  r  be  cut  by  a  horizontal  line  into  two  convex  polygons  r  1  and 
r2.  We  have: 

(i)  If  r  is  1-supported,  then  both  r  1  and  r2  are  1-supported. 

(ii)  If  r  is  2-supported,  then  at  least  one  of  r  1  and  r2  is  1-supported. 

In  the  description  of  the  data  structure  for  dynamic  point  location  we  will  assume  that  the 
subdivision  91  is  2-supported.  As  shown  in  the  following,  this  is  not  restrictive.  Let  91*  be  the 
subdivision  obtained  from  91  by  adding  to  each  region  r  which  is  not  1 -supported  the  edge  from 
LOW  (r)  to  HIGH  (r).  (Notice  that  91  may  have  degenerate  regions  of  measure  zero.)  Clearly 
91*  is  1-supported.  The  answer  to  a  point  location  query  in  91*  can  be  immediately  converted  to 
an  answer  in  91.  Regarding  updates,  we  have 
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Lemma  4.2  For  each  update  operation  performed  on  SR  there  exists  a  sequence  of  at  most  four 
operations  which  update  SR*  so  that  after  each  intermediate  operation  the  resulting  subdivision  is 
2-supported  (possibly  1-supported). 

Proof:  Consider  operation  INSERTSEGMENT  (e,vi,v2,r;ri,r2)  performed  on  91.  The 
sequence  of  operations  in  91*  is  as  follows  (see  Fig.  4.4): 

( 1)  insert  edge  s  t  =  ( LOW  (r  ^  ),HIGH  (r  i )); 

(2)  insert  edge  s2  =  (LOW  (r2)MIGH  (r 2)); 

(3)  remove  edge  s  =  ( LOW  (r),H!GH  (r)),  unless  s=s\  or  r  =  s2; 

(4)  insert  edge  e. 


HIGH  (r) 


(a) 


HIGH  (r) 


(b) 


Figure  4.3  (a)  A  1-supported  polygon,  (b)  A  2-supported 
polygon. 
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(We  assume  that  an  insertion  is  not  performed  if  the  edge  to  be  inserted  already  exists.)  Similar 
arguments  hold  for  operations  REMOVESEGMENT. 

Regarding  operation  INSERTPOINT  (v.e.e\,ei\.  let  r  j  and  r2  be  the  regions  to  the  left  and 
right  of  r,  respectively.  The  sequence  of  operations  in  SR*  is  as  follows: 

(1)  insert  vertex  v; 

(2)  insert  edge  s  \  =  ( LOW  (r  j  ),HIGH  (r  \ )); 

(3)  insert  edge  $2  =  (LOW  {r2),HIGH  (r^)). 

Again,  the  case  of  operation  REMOVEPOINT  is  similar. 

Finally,  operation  MOVEPOINT  does  not  affect  the  supportedness  of  any  region  of  91,  so 
that  no  further  operation  is  required  in  91*.  □ 


4.3.  Trapezoids 

Throughout  this  section  we  denote  by  91  a  2-supported  convex  subdivision  whose  n  vertices 
are  on  the  lines  of  L.  We  show  that  the  trapezoid  method  for  point  location  [46]  can  be  modified 
to  perform  efficiently  updates  on  91. 

The  trapezoid  method  is  based  on  a  recursive  decomposition  of  the  subdivision  91  into  com¬ 
ponents,  called  trapezoids.  A  trapezoid  x  is  a  convex  subdivision  whose  external  boundary  is  a 
quadrilateral  with  two  (not  necessarily  bounded)  horizontal  sides,  which  belong  to  the  lines  of  L. 
We  denote  the  four  sides  of  x  by  LEFT  (x),  RIGHT  (x),  BOT  (x),  and  TOP  (x).  Let  BOT  (x)  be  on 
line  /,-,  and  TOP  (x)  on  lj.  The  median  line  of  X,  denoted  MEDIAN  (x),  is  the  line  /*  with 
k  =  [(/  +y)/2j  .  The  subdivision  91  is  the  trapezoid  with  BOT  (91)=/ o,  TOP  (91  )=/y,  and  having 
LEFT  (x)  and  RIGHT  (x)  at  infinity. 

A  spanning  edge  of  trapezoid  x  is  an  edge  of  x  with  endpoints  on  the  top  and  bottom  sides 
of  x.  A  spanning  edge  s  partitions  trapezoid  x  into  two  trapezoids,  denoted  Xi  and  Xg,  with 
s  =  RIGHT  (Xl)=LEFT  (Xg).  This  is  called  a  vertical  cut  (see  Fig.  4.5(a)).  If  trapezoid  x  does 
not  have  spanning  edges,  we  decompose  it  by  means  of  a  horizontal  line  through  its  median  line, 
and  obtain  trapezoids  xB  and  Xj,  with  MEDIAN  (x)=TOP  (xb)=BOT  (X7).  This  is  called  a  hor¬ 
izontal  cut  (see  Fig.  4.5(b)).  A  trapezoid  can  always  be  decomposed  unless  it  is  empty.  We 
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assume  that  a  vertical  cut  takes  precedence  over  a  horizontal  cut,  so  that  the  decomposition  pro¬ 
cess  is  uniquely  determined  to  within  the  specification  of  the  sequence  of  consecutive  vertical 
cuts  of  the  same  trapezoid. 

The  above  decomposition  process  can  be  conveniently  represented  by  means  of  a  string 
notation.  Namely,  a  vertical  cut  can  be  represented  by  the  transformation  x-»X/,  s  x#,  where  s  is 
the  spanning  edge;  and  a  horizontal  cut  by  the  transformation  x— Kx^)/  (X7),  where  /  is  the 
median  line  of  t.  Finally,  an  empty  trapezoid  x  is  eliminated  using  the  transformation  T— »£, 
where  e  is  the  null  string.  Hence,  the  recursive  decomposition  of  a  trapezoid  x  can  be  viewed  as 
a  sequence  of  rewritings  using  one  of  the  above  rules.  In  this  framework,  trapezoids  correspond 
to  nonterminal  symbols  (variables),  while  spanning  edges  and  lines  correspond  to  terminal  sym¬ 
bols.  The  initial  string  is  the  nonterminal  symbol  x.  The  final  string  is  a  sequence  of  spanning 


Figure  4  J  (a)  Example  of  vertical  cut.  (b)  Example  of 
horizontal  cut. 
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edges,  lines,  and  parentheses.  For  example,  the  recursive  decomposition  of  the  trapezoid  shown 
in  Fig.  4.6  is  described  by  the  string 

(  (« 1  «3)  / 1  (C2«4«7)£8g9  )  1 2  («2(«4)  *3  («5  «6)«10ell  )• 

Notice  that  the  parentheses  have  been  introduced  in  the  horizontal  decomposition  to  eliminate 
ambiguities  in  the  parsing  of  such  strings. 

The  canonical  decomposition  of  a  trapezoid  x  with  spanning  edges  is  defined  by  (see  Fig. 


X  -4  ToOjTiC^  •  •  •  Okxk 


where 


(i)  each  x,-,  i  *  1,  •  •  ♦  ,k  - 1,  is  a  nonempty  trapezoid  without  spanning  edges; 

(i ii )  each  of  to  and  xk  either  is  empty  or  is  a  nonempty  trapezoid  without  spanning  edges;  and 
(w)  each  a7,  j - 1,  •  •  •  ,k,  is  a  maximal  sequence  of  spanning  edges  that  delimit  empty  tra- 


«2\  e5\ 


h  /— 


Figure  4.7  Canonical  decomposition  of  a  trapezoid. 


4.4.  Data  Structure 

The  data  structure  for  point  location  makes  use  of  the  so-called  biased  binary  tree  [2], 

which  is  a  dynamic  binary  search  tree  that  represents  a  sorted  list  of  weighted  items. 

Specifically,  the  i-th  item  in  the  sorted  list,  with  weight  w,-,  is  represented  by  the  i-th  leaf  of  the 
k 

tree.  Let  w  =  £ w;  be  the  total  weight  of  the  items  in  the  tree.  Some  important  properties  of 
i=l 

biased  binary  trees,  which  are  relevant  to  this  work,  are  listed  below: 

( 1 )  The  depth  of  the  i-th  leaf  is  at  most  log  (w/w; ) + 2. 

(2)  Dynamic  operations  such  as  insertion/deletion  of  an  item,  change  of  weight  of  an  item,  and 
split/splice  of  biased  trees,  can  be  performed  in  time  O  Gog  w).  (For  the  splice  operation,  w 
is  the  total  weight  of  the  trees  to  be  spliced). 

The  point  location  data  structure  consists  of  a  main  component  and  of  an  auxiliary  com¬ 
ponent.  The  main  component  is  a  binary  tree  3  (SR),  called  trapezoid  tree,  which  represents  the 
recursive  decomposition  of  the  subdivision  SR. 

The  nodes  of  the  trapezoid  tree  3  (SR)  are  of  three  types  (see  Fig.  4.8): 
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(C) 


Figure  4.8  Nodes  of  a  trapezoid  tree:  (a)  V-node;  (b)  O- 
node;  (c)  e-node. 

(1)  V-node:  each  such  node  is  associated  with  a  nonempty  trapezoid  x  which  is  partitioned  by 
a  horizontal  cut;  it  stores  the  line  MEDIAN  (x). 

(2)  O-node :  each  such  node  is  associated  with  a  nonempty  trapezoid  x  which  is  partitioned  by 
vertical  cuts;  it  stores  a  maximal  sequence  o  of  spanning  edges  in  the  canonical  decomposi¬ 
tion  of  x,  represented  by  a  balanced  tree. 

(3)  e-node:  each  such  node  is  associated  with  an  empty  trapezoid,  and  is  a  leaf  of  the  tra¬ 
pezoid  tree. 

In  addition,  each  node  ji  stores  an  integer  weighty),  which  denotes  the  number  of  vertices  in  the 
interior  of  the  trapezoid  associated  with  |i.. 

Each  node  of  the  trapezoid  tree  is  uniquely  associated  with  a  trapezoid,  so  that  we  will 
sometimes  use  the  same  name  for  a  node  of  the  tree  and  for  the  associated  trapezoid. 

The  trapezoid  tree  3  (x)  for  a  nonempty  trapezoid  x  is  recursively  defined  as  follows  (see 
Fig.  4.9): 

(1)  If  x  has  no  spanning  edges,  the  root  of  3  (x)  is  a  V-node,  and  the  left  and  right  subtrees  of 
3  (x)  are  trapezoid  trees  for  Xb  and  Xj,  respectively. 
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Figure  4.9  Trapezoid  tree  for  the  trapezoid  of  Fig.  4.7. 

(2)  Otherwise  (x  has  spanning  edges),  we  consider  the  canonical  decomposition  of  t, 
XoOlXi02  •  •  •  a*x*,  and  define  T(x)  as  a  biased  binary  tree  for  the  items  Xo,  •  •  •  .x*.  Tree 
T(x)  is  called  the  decomposition  tree  of  x.  The  i-th  leaf  of  7*(x)  is  a  V-node  for  trapezoid  X;. 
For  the  special  case  when  i  =0  or  /  =k  and  x,-  is  empty,  leaf  /  is  a  e-node,  and  has  nominal 
unit  weight.  The  i-th  internal  node  of  T,  denoted  p.„  stores  the  sequence  o,  of  spanning 
edges,  and  is  associated  with  the  trapezoid  formed  by  the  union  of  the  trapezoids  associated 
with  the  leaves  of  the  subtree  of  T(x)  rooted  at  |i;.  Each  node  of  T(x)  is  also  equipped  with 
pointers  pred  and  succ,  which  form  a  doubly-linked  list  of  the  nodes  of  T(x)  in  their  left-to- 
right  order.  Finally,  the  trapezoid  tree  3  (x)  is  obtained  by  replacing  each  leaf  x,  of  T(x) 
with  the  trapezoid  tree  3  (x,). 
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It  is  interesting  to  observe  that  each  edge  e  of  the  subdivision  9?  may  be  stored  in  the 
sequences  of  spanning  edges  of  several  O-nodes  of  tree  3  (91).  Each  such  representative  is 
called  a  fragment  of  edge  e.  The  decomposition  of  e  into  fragments  is  performed  according  to 
the  well-known  segment-tree  scheme  [47,  pp.  13-15],  so  that  the  number  of  fragments  of  each 
edge  is  at  most  2  log  N  - 1 . 

Theorem  4.2  The  space  requirement  of  the  trapezoid  tree  3  (9?)  is  O  (n  log  AO- 

Proof:  Let  the  order  of  a  V-node  p.  be  the  number  of  V-nodes  on  the  path  from  p  to  the  root  of 
3  (91),  including  p.  Because  of  the  median  line  decomposition  rule,  the  order  of  p  is  at  most 
logN.  Also,  for  any  p  £log  N,  there  are  at  most  n  nodes  of  order  p  because  each  such  node 
corresponds  to  a  trapezoid  with  at  least  one  vertex.  Hence,  the  total  number  of  V-nodes  in  3  (91) 
is  at  most  n  log  N.  Regarding  O-nodes,  we  observe  that  each  edge  is  decomposed  according  to 
the  well-known  segment-tree  scheme,  so  that  the  total  space  requirement  for  the  O-nodes  and 
their  sequences  of  spanning  edges  is  0  (n  log  IV).  Finally,  we  can  charge  at  most  two  e-nodes  to 
each  V-node,  so  that  the  total  space  for  the  trapezoid  tree  3  (91)  is  O  (n  log  N).  □ 

Theorem  4.3  The  depth  of  the  trapezoid  tree  3  (91)  is  O  (log  n  +  logN). 

Proof:  Without  significant  loss  of  generality,  assume  that  N  is  a.  power  of  2.  Let  height  (x)  be 
the  number  of  elementary  slabs  spanned  by  a  trapezoid  x.  Clearly,  height(x)  SIV.  We  show  that 
for  a  V-node  of  3(91)  with  associated  trapezoid  x,  the  depth  of  3(x)  is  at  most 
log  weight(x)  +  3  log  height(x)  + 1 .  (For  simplicity  of  notation,  we  assume  that  log  0=1.)  The 
proof  is  by  induction  on  height(x).  The  base  of  the  induction,  height {x)  =  2,  is  immediate,  since 
the  left  and  right  subtrees  of  3  (x)  consist  each  of  a  O-node  and  two  e-nodes.  For  the  inductive 
step,  define  w  =  weight  (x)  and  h=height{x).  Consider  the  deeper  of  the  two  subtrees  of  x,  say, 
the  one  associated  with  trapezoid  x&.  Let  x,  be  a  leaf  of  the  decomposition  tree  T(Xg),  and  define 
Wi  =  weight(Xi )  and  hi  -  height {Xi).  Since  T(Xg)  is  a  biased  binary  tree,  leaf  xx  has  depth  bounded 
by  log(w/w,)  +  2.  By  the  inductive  hypothesis,  the  depth  of  the  trapezoid  tree  3  (x,)  is  at  most 
log  w,-  +  3  log  hi  + 1.  Hence,  the  total  depth  of  the  trapezoid  tree  3  (x)  is  at  most 


■  (log  (w/w,  )  +  2)+ 1  +  (log  w;  +  3  log  hi  + 1). 

|  Since  hi  =  A/2,  this  is  exactly  log  w  +  3  log  A  + 1. 

i 

1 
' 

l 

[  The  auxiliary  component  of  the  data  structure,  called  dictionary ,  contains  a  record  for  each 

|  vertex,  region,  and  edge  of  91.  The  record  for  a  vertex  v  stores  the  name  of  v  and  its  coordinates. 

The  record  for  a  region  r  stores  the  name  of  r,  pointers  to  the  records  of  vertices  LOW  (r)  and 
!  HIGH  (r),  and  pointers  to  two  balanced  search  trees,  Ichain  (r)  and  rchain(r),  respectively, 

•  representing  the  edges  of  the  left  and  right  chain  of  r  in  bottom-to-top  order.  The  record  for  an 

r  edge  e  stores  the  name  of  e,  pointers  to  the  records  of  the  endpoints  of  e,  denoted  LOW  (e)  and 

;  HIGH  (e),  and  pointers  If  ace  (e)  and  rface(e)  to  the  representatives  of  e  in  the  trees  rchain(r\) 

[  and  lchain(r 2)  such  that  r*  and  r2  are  the  regions  on  the  left  and  right  sides  of  e,  respectively. 

I  Dearly,  the  dictionary  takes  O(n)  space.  The  dynamic  maintenance  of  the  dictionary  in  the 

various  update  operations  can  be  easily  performed  in  O(Iogzi)  time,  and  will  not  be  further 
described. 

I  4.5.  Query 

The  query  algorithm  is  essentially  the  same  as  in  the  static  case  [46].  It  consists  of  tracing  a 
path  in  3  (31).  The  actions  taken  at  the  current  node  p.  are  as  follows: 

I  (1)  If  p.  is  a  e-node,  we  stop. 

(2)  If  p  is  a  V-node,  we  discriminate  the  query  point  q  with  respect  to  the  median  line 
MEDIAN  (t)  of  the  trapezoid  x  associated  with  p.  We  proceed  to  the  left  or  right  child  of  p 
depending  upon  whether  q  is  below  or  above  MEDIAN  (x).  (If  q  is  on  MEDIAN  (x)  we 
proceed  to  the  left  child.) 

(3)  If  p  is  a  O-node,  we  discriminate  q  with  respect  to  the  sequence  of  spanning  edges  associ¬ 
ated  with  p,  a(p) = s  1 ,  •  •  •  ,sp.  If  q  lies  to  the  left  of  edge  jj,  we  sets*  :=sj  and  proceed  to 
the  left  child  of  p.  If  q  lies  to  the  right  of  sp,  we  set  :=sp  and  proceed  to  the  right  child 
of  p.  If  q  lies  between  two  edges  of  o(p),  say  s;  and  s,  + 1 ,  we  set  si  := s,-  and  s*  :=  s,  + 1 ,  and 
stop. 
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The  region  r  reported  is  the  one  which  lies  on  the  right  side  of  edge  Si  and  on  the  left  side  of 
edge  j/j. 

The  time  spent  is  0(1)  for  the  intermediate  nodes  on  the  search  path,  and  O  (log  n)  for  the 
last  node.  The  final  computation  of  r  from  Si  and  Sr  can  be  done  by  accessing  the  representa¬ 
tives  of  si  and  Sr  in  the  trees  lchain(r)  and  rchain(r)  pointed  by  rface(si)  and  tface(sR),  and 
walking  up  to  the  root  of  such  trees,  which  takes  0  Gog  n)  time.  Since  the  height  of  tree  3  (91)  is 
logarithmic,  we  conclude 

Theorem  4.4  The  complexity  of  a  point  location  query  is  O  Gog  n  +  log  IV). 

4.6.  Insertion  of  Edges 

First,  we  consider  operation  INSERTSEGMENT  (e,vl,v2,r;rl,r2).  The  insertion  algorithm 
is  described  by  means  of  a  recursive  procedure,  whose  actions,  depending  on  the  type  of  the 
current  node  |i,  are  as  follows. 

1.  If  p  is  a  O-node  associated  with  the  string  o(p)=s  i  •  •  sp,  we  locate  e  in  a(p). 

1.1  li  e  is  to  the  left  of  s  i ,  we  recursively  call  the  procedure  on  the  left  child  of  p. 

1.2  If  e  is  to  the  right  of  sp,  we  recursively  call  the  procedure  on  the  right  child  of 

L3  Otherwise,  e  lies  between  two  edges  of  o(p),  say  st  and  Sj  +  j.  In  this  case,  we  insert  e 
into  a(ji)  between  s,-  and  s,  + 1 ,  corresponding  to  the  string  transformation 

Sl  "  Si  Si  + 1  ’  ’  ’  Sp  S 1  '  ’  *  S,"  €  Si  + 1  *  *  *  Sp . 

Note  that  this  case  can  only  occur  when  e  spans  the  slab  of  oQi),  for,  in  the  interior  of  the 
corresponding  region  of  the  plane,  there  are  no  vertices  of  SR. 

2.  If  p.  is  a  e-node,  we  know  that  it  is  either  the  leftmost  or  rightmost  leaf  of  a  decomposi¬ 
tion  tree.  In  the  former  case,  edge  e  is  added  to  the  left  of  the  sequence  of  spanning  edges  of  the 
O-node  next([L)  immediately  following  n.  In  the  latter  case,  edge  e  is  added  to  the  right  of  the 
sequence  of  spanning  edges  of  the  O-node  prediil)  immediately  preceding  |x. 

3.  If  n  is  a  V-node  associated  with  trapezoid  t,  we  determine  first  whether  e  spans  %  or  not 


3.1  If  e  docs  not  span  x,  then  if  LOW  (e)  is  below  MEDIAN  (x)  we  recursively  call  the  pro¬ 
cedure  on  the  left  child  of  p,  and  if  HIGH  (e)  is  above  MEDIAN  (x)  we  recursively  call  the  pro¬ 
cedure  on  the  right  child  of  p.  (Note:  here  is  where  the  edge  may  be  subjected  to  the  segment- 
tree  fragmentation.) 

3.2  If  e  spans  x,  we  determine  the  leftmost  and  rightmost  regions  of  x,  i.e.,  the  region  ri  to 
the  right  of  edge  LEFT  (x)  and  the  region  rR  to  the  left  of  RIGHT  (x).  Such  regions  can  be  com¬ 
puted  by  observing  that  LEFT  (x)  is  the  last  spanning  edge  in  the  sequence  a(pred(\i)),  and 
RIGHT  (x)  is  the  first  spanning  edge  in  the  sequence  c(/iexr(p)).  We  have  three  further  sub¬ 
cases: 

3.2.1  If /-=/-£,  i.e.,  edge  e  is  to  the  left  of  each  vertex  of  x,  then  edge  e  is  added  to  the  right 
of  the  sequence  of  spanning  edges  of  the  O-node  pred(p)  immediately  preceding  p,  and  no 
further  action  is  required.  (If  pred{ p)  =  nil ,  a  new  O-node  and  a  new  e-node  arc  created,  and  the 
decomposition  tree  containing  leaf  p  is  rebalanced.)  This  case  corresponds  to  the  syntactic 
transformation  x— x. 

3.2.2  Ur=rR,  i.e.,  e  is  to  the  right  of  each  vertex  of  x,  then  edge  e  is  added  to  the  left  of  the 
sequence  of  spanning  edges  of  the  O-node  nexr(p)  immediately  following  p,  and  no  further 
action  is  required.  (If  next  (p)= nil,  a  new  O-node  and  a  new  E-node  are  created,  and  the  decom¬ 
position  tree  containing  leaf  p  is  rebalanced.)  This  case  corresponds  to  the  syntactic  transforma¬ 
tion  x— >xe,  and  is  symmetric  to  case  3.2.1. 

3.2 3  If  r*rL  and  r*rR,  we  have  that  e  lies  inside  x,  which  will  cause  a  more  substantial 
update  of  the  data  structure  (see  Fig.  4.10).  In  fact,  while  x  was  formerly  horizontally  decom¬ 
posed,  it  has  now  to  be  vertically  decomposed  by  means  of  a  vertical  cut  along  edge  e.  First,  we 
call  recursively  the  procedure  on  both  the  left  and  right  children  of  p.  At  this  point,  both  tra¬ 
pezoids  Xfl  and  Xj  have  a  vertical  cut  along  edge  e,  corresponding  to  the  string 

(Xnaneai2Xi2)/  (X21 021^022 122)- 

To  comply  with  the  rule  that  vertical  cuts  take  precedence  over  horizontal  cuts,  we  have  to  res¬ 
tructure  the  tree  according  to  the  string 
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(*11  <*11  (?l\Gl\)e  (<*12*12)1  (022*22)- 

This  is  done  as  follows: 

(1)  Let  | ii  and  4#  be  the  left  and  right  child  of  4,  respectively.  Find  node  m  of  T([iL)  associ¬ 
ated  with  the  sequence  Oi  =  Cn£Cfi2-  Similarly,  find  node  112  of  associated  with  the 
sequence  <*2  =<*2ie<*22- 

(2)  Split  T(\il)  and  T(\Lr)  at  nodes  |ii  and  42,  respectively.  This  yields  trees  3  n,  3 12.  3 21, 

and  3  22- 

(3)  Split  the  sequences  o(|ii)  and  0(1x2)  at  edge  e.  This  yields  the  sequences  On,  O12,  o2i,  and 
O22. 

(4)  If  On  5*0,  create  a  O-node  |Xn  for  it,  and  add  |Xn  and  a  e-node  to  the  right  of  3  n-  If 
0219*0,  create  a  O-node  421  for  it,  and  add  1121  and  a  e-node  to  the  right  of  3  21  •  If  Ci2*0, 
create  a  O-node  p.12  for  it,  and  add  412  and  a  e-node  to  the  left  of  3 12-  If  0229*0,  create  a 
O-node  422  f°r  it,  and  add  422  and  a  e-node  to  the  left  of  3  22- 

(5)  Create  a  V-node  Vi  (associated  with  the  portion  of  x  to  the  left  of  e),  with  left  subtree  3  n 
and  right  subtree  3  21  • 

(6)  Create  a  V-node  v2  (associated  with  the  portion  of  x  to  the  right  of  e),  with  left  subtree  3 12 
and  right  subtree  3  22- 

(7)  Replace  4  with  a  O-node  v  having  o(v)  =  e,  left  child  Vi ,  and  right  child  v2.  Finally,  rebal¬ 
ance  the  decomposition  tree  containing  node  v  as  a  leaf. 

This  completes  the  description  of  the  algorithm  for  operation  INSERTSEGMENT. 

We  now  analyze  the  time  complexity  of  operation  INSERTSEGMENT.  Let  r*  be  the  inter¬ 
section  of  region  r  and  the  horizontal  slab  spanned  by  segment  e  (see  Fig.  4.11).  Notice  that, 
since  r  is  2-supported,  so  is  r*  .  The  insertion  of  e  into  the  trapezoid  tree  can  be  viewed  as  a 
visit  of  a  subgraph  of  3  (9t)  (see  Fig.  4.12).  Specifically,  this  subgraph  consists  of  an  initial  path 
originating  at  the  root  and  consisting  of  O-nodes  such  that  e  is  always  external  to  the  sequence  of 
spanning  edges  (Case  1.1  or  1.2),  and  of  V-nodes  such  that  e  never  crosses  the  median  of  the 
corresponding  slab  (Case  3.1).  For  either  type  of  node,  0(1)  time  is  spent  at  each  node.  The 
visit  of  this  initial  path  terminates  at  a  V-node  (referred  to  as  the  fork),  where,  for  the  first  time,  e 
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crosses  the  median. 


At  the  fork,  e  is  passed  on  to  both  children,  i.e.,  the  visit  continues  with  two  distinct  paths, 
called  the  left  spine  and  right  spine.  Let  us  consider  the  left  spine  (the  right  spine  is  similarly 
analyzed).  It  is  easily  realized  that  this  path  contains  a  set  of  distinguished  V-nodes  such  that, 
for  trapezoid  t  associated  with  any  of  them,  y (TOP  (x))£y  (HIGH  (e))  and 
y(BOT (t))  <y(LOW  (e))£y (MEDIAN  (t)).  This  implies  that  e  spans  the  trapezoid  X?  associ¬ 
ated  with  the  right  subtree  of  the  distinguished  V-node  x.  The  decomposition  tree  T(xt)  is  called 
an  allocation  subtree  for  edge  e.  The  slab  spanned  by  the  roots  of  all  allocation  subtrees  of  e 
form  a  partition  of  r*  into  a  set  of  convex  polygons,  of  which  it  is  easily  realized  that  at  most 
one  is  2-supported  while  all  the  others  are  1-supported  (see  Fig.  4.1 1). 

Each  distinguished  node  of  the  left  spine  is  a  branching  point  for  the  visit,  since  the  pro¬ 
cedure  is  called  recursively  on  both  children.  The  left  child  is  on  the  left  spine,  while  the  right 
child  is  the  root  of  an  allocation  subtree  for  e.  Within  the  allocation  subtree  the  visit  continues 
along  a  path,  called  allocation  trail  until  we  reach  a  node,  called  allocation  target  for  which  one 
of  the  Cases  1.3, 2,  or  3.2  occurs.  The  time  spent  at  each  node  of  the  left  spine  and  at  each  node 


slabs  of  left  children 
of  distinguished  nodes 
on  right  path 


median  of  fork 

slabs  of  right  children 
of  distinguished  nodes 
on  left  path 


Figure  4.11  Region  r*  in  the  analysis  of  operation  INSERTSEGMENT. 


of  the  allocation  trail,  target  excluded,  is  O  (1). 

Consider  now  an  allocation  target  p  of  e.  If  p  is  an  e-node,  then  its  handling  takes  O  (log  n) 
time  (Case  2).  If  p  is  an  O-node,  then  in  O  (log  n )  time  e  is  inserted  into  the  corresponding  a(p) 
(Case  1-3).  Finally,  if  p  is  a  V-node,  we  are  in  Case  3.2.  If  the  intersection  of  r  with  the  slab  of 
p  is  a  1-supported  polygon,  then  we  have  either  Case  3.2.1  or  Case  3.2.2,  since  e  is  adjacent  to  a 
spanning  edge  of  this  slab.  For  such  cases,  the  processing  of  p  takes  time  O  (logn).  Else  (when 
the  intersection  of  r  with  the  slab  of  p  is  a  2-supported  polygon),  we  are  in  Case  3.2.3.  From 
Lemma  4.1,  we  have  that  this  case  can  occur  for  at  most  one  target  node,  called  henceforth  spe¬ 
cial  target. 

Let  p  be  the  special  target,  and  let  t  be  the  associated  trapezoid.  In  the  processing  of  p,  we 
recursively  call  the  procedure  on  both  children  of  p,  so  that  the  visit  branches  again  into  two 
paths.  Since  edge  e  spans  both  xB  and  Xj,  each  path  will  encounter  a  node  for  which  one  of  the 
Cases  1.3,  2,  or  3.2  occurs.  However,  by  Lemma  4.1,  at  most  one  of  the  polygons  obtained  by 
intersecting  region  r  with  xB  and  Xj  is  2-supported,  so  that  Case  3.2.3  can  occur  again  for  at  most 
one  such  node.  Finally,  after  the  recursive  calls  have  been  executed,  we  perform  the  restructur¬ 
ing  specified  in  the  description  of  Case  3.2.3,  which  takes  0  (log  n)  time.  Let  p  be  the  order  of 
p,  i.e.,  the  number  of  nodes  on  the  path  from  p  to  the  root  The  work  done  for  completing  the 
visit  originating  at  p  requires  in  the  worst  case  time  0(logn)  plus  the  time  for  executing  one 
recursive  call  to  the  procedure  on  a  node  of  order  p  - 1.  Since  the  order  of  a  node  is  at  most 
log  N,  this  sums  up  to  O  (log  n  log  N). 

We  summarize  the  preceding  discussion  as  follows: 

(i)  The  total  time  spent  at  the  nodes  of  the  path  from  the  root  to  the  fork  and  at  the  nodes  of  the 
left  and  right  spine  is  O  (log  n  +  log  N). 

(«)  The  total  time  spent  in  traversing  each  allocation  trail  is  OQogrt),  for  a  total  time 
O  (log  n  log  N). 

(Hi)  The  time  spent  for  processing  each  nonspecial  target  node  is  O(logn),  for  a  total  time 
0(logn  log  AO- 

(iv)  The  time  spent  for  processing  the  special  target  and  the  subsequent  visit  is  O  (log  n  log  N). 


We  therefore  conclude 


Theorem  4.5  The  complexity  of  operation  INSERTSEGMENT  (e,vi,v2tr;ri,r2)  is 
OQogn  log N). 

Remark.  The  underlying  segment-tree  structure  of  the  method  effects  the  fragmentation  into  a 
logarithmic  number  of  subsegments  of  an  edge  to  be  inserted.  Each  of  these  subsegments,  how¬ 
ever,  may  not  be  insertable  as  a  single  entity,  if  it  does  not  occur  next  to  existing  spanning  edges 
of  the  corresponding  slab.  Note,  however,  that  such  a  spanning  edge  would  always  exist  if  each 
region  of  91  were  a  1 -supported  polygon.  Relaxing  such  strong  constraint  to  2-supportedness, 
we  effectively  create  a  “channel”  (a  sort  of  “virtual  diagonal”)  between  LOW  (r)  and 
HIGH  (r)  of  each  polygon  r  of  91,  which  guarantees  that  at  most  one  segment  fragment  (the  one 
reaching  the  special  target)  is  further  subdivided  into  a  number  of  subsegments  at  most  loga¬ 
rithmic  in  N,  before  being  restored  to  its  original  length  as  specified  by  the  string  transformation 
given  for  Case  3.2.3. 

4.7.  Other  Update  Operations 

The  algorithm  for  operation  REMOVESEGMENT  (e,v  i ,v2,r t ,r2;r)  is  similar  to  the  one 
for  the  INSERTSEGMENT  operation,  and  for  brevity  we  omit  its  description. 

Theorem  4.6  The  complexity  of  operation  REMOVESEGMENT  (e,v1,V2,ri,r2;r)  is 
O  (log  n  log  AO- 

We  consider  now  operation  INSERTPOINT  (v,e;e1,e2)-  The  algorithm  consists  of  two 
phases.  The  first  phase  is  performed  by  a  recursive  procedure  consisting  of  the  following  steps: 

(1)  Using  the  query  algorithm,  we  search  for  a  node  p.  of  the  trapezoid  tree  such  that  v  lies  in 
the  trapezoid  x  of  p.,  and  e  is  a  spanning  edge  of  o(p.) = s  i ,  •  •  •  ,sp,  i.e.,  e  =r,-  for  some  i. 

(2)  We  perform  some  restructurings  at  node  p.  (to  be  specified  later),  which  correspond  to 
inserting  v  onto  line  MEDIAN  (x). 

(3)  We  recursively  call  the  procedure  if  v  is  not  on  line  MEDIAN  (x). 
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The  restructurings  performed  in  Step  2  depend  on  the  position  of  edge  e  in  the  sequence 
We  distinguish  three  cases: 

1.  If  e  =  Si  is  a  nonextreme  element  of  the  sequence  o(p),  i.e,  i*l,p,  the  restructuring  is 
dictated  by  the  string  transformation 

Si_i«sl  +  1  (c)/(e)Si+ 1. 

Omitting  the  details,  node  h  is  replaced  in  the  decomposition  tree  by  nodes  Hi ,  p/.  and  H2>  where 
Hi  and  H2  are  O-nodes  with  sequences  ct(Hi)=Si.  ■  •  •  ,r,-i  and  o(H2)=^+i.  •  •  •  ,sp,  and  h'  is  a 
V-node.  In  turn,  h*  is  the  root  of  an  elementary  decomposition  tree,  whose  left  and  right  subtree 
consist  each  of  a  O-node  with  sequence  e,  and  of  two  e-nodes. 

2.  If  e=s  1  is  the  leftmost  element  of  c(h)  and  p£2,  the  restructuring  is  dictated  by  the 
string  transformation 

(xi)l(x2)es2  -*  (Ti  e)l(x2e)s2. 

Let  \l'  =pred(\i).  Node  \i!  is  associated  with  trapezoid  xf  to  the  left  of  e.  Edge  e  is  inserted  into 
the  left  and  right  children  of  h'  as  in  the  algorithm  for  INSERTSEGMENT. 

3.  If  e  -sp  is  the  rightmost  element  of  <j(h)  and  p  £  2,  the  restructuring  is  dictated  by  the 
string  transformation 

5^  _  1 « (Xi )  /  (X2)  — »  _  i  (e  xj)/(ex2). 

This  case  is  symmetric  to  Case  2  above. 

4.  If  e  -S\  -sp  is  the  unique  element  of  a(H)>  the  restructuring  is  dictated  by  the  string 
transformation 

(Xu)  l  (x2i)e  (X12)  /  (X22)  -» (Xu  e  x12)  /  (x21  e  xn)- 

This  restructuring  merges  the  two  trapezoids  associated  with  the  V-nodes  Hi  =pred{\i)  and 
H2  -  nexr(n)  into  a  unique  trapezoid,  represented  by  a  new  node  h'.  The  tasks  to  be  performed 
are  analogous  to  the  ones  for  Case  3.2.3  of  the  edge  insertion  algorithm,  and  omitted  for  brevity. 

This  completes  the  description  of  the  restructurings  performed  in  the  three  cases.  In  Step  3,  if  v 
is  below  MEDIAN  (x)  we  recursively  call  the  procedure  on  the  left  subtree  of  \i',  while  if  v  is 


above  MEDIAN  (x)  we  recursively  call  the  procedure  on  the  right  subtree  of  p'. 

A  detailed  time  complexity  analysis  of  the  first  phase  can  be  conducted  using  techniques 
similar  to  the  ones  presented  in  the  preceding  section.  In  summary,  the  recursion  depth  is  at 
most  O  (log AO,  and  the  time  for  Step  2  is  O  (log  n)  for  all  but  (possibly)  one  node,  for  which  it  is 
O  (log  n  log IV).  Hence  the  first  phase  takes  time  O  (log  n  log  N). 

The  second  phase  rebalances  the  decomposition  trees  embedded  in  3  (31)  for  which  the 
weight  of  some  leaf  has  increased  by  one.  There  are  at  most  0  (log  AO  such  trees,  so  that  the 
second  phase  also  takes  time  O  (log  n  log  N). 

We  therefore  conclude 

Theorem  4.7  The  complexity  of  operation  INSERTPOINT  (y,e;e i,e2)  is  O  (log/i  log  N). 

The  algorithm  for  operation  REMOVEPOINT  (v,e  j  ,e2;e)  is  similar  to  the  one  for  operation 
INSERTPOINT,  and  we  have 

Theorem  4.8  The  complexity  of  operation  REMOVEPOINT  {y,e i,e2;«)  is  O  Gog  n  log N). 

Finally,  as  regards  operation  MOVEPOINT  (v\x),  we  observe  that  it  does  not  affect  the  tra¬ 
pezoid  tree.  Hence,  it  can  be  simply  carried  out  by  updating  the  dictionary  in  time  O  (log  n). 

In  conclusion,  we  have  proved  the  following  theorem: 

Theorem  4.9  Let  31  be  a  convex  subdivision  whose  n  vertices  lie  on  a  fixed  set  L  of  N  horizon¬ 
tal  lines.  There  exists  a  dynamic  point  location  data  structure  for  91  with  space  O  ( n  log  AO  and 
query  time  0(log/i  +  logA0,  which  allows  for  insertion/deletion  of  vertices  and  edges  in  time 
O  (log/i  log  AO,  all  time  bounds  being  worst-case. 

4.8.  Extensions 

The  analysis  of  the  query  and  update  operations  shows  that  the  essential  requirement  for  the 
method  is  the  monotonicity  and  2-supportedness  of  the  regions  of  the  subdivision,  while  the  con¬ 
vexity  requirement  can  be  relaxed.  Hence,  we  obtain  the  following  corollary: 


Corollary  4.1  Let  9?  be  a  monotone  2-supported  subdivision  whose  n  vertices  lie  on  a  fixed  set 
L  of  N  horizontal  lines.  There  exists  a  dynamic  point  location  data  structure  for  SR  with  space 
O  ( n  log  N)  and  query  time  O  (log  n  +  logiV),  which  allows  for  insertion/deletion  of  vertices  and 
edges  in  time  O  (log  n  log  N),  all  time  bounds  being  worst-case. 
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CHAPTER  5 


DYNAMIC  POINT  LOCATION:  RECONSTRUCTION  METHOD 


5.1.  Introduction 

In  this  section,  we  present  a  fully  dynamic  technique  for  point  location  in  a  triangulation 
with  n  vertices.  This  method  features  a  simple  and  versatile  O(n)- space  data  structure,  which 
can  be  used  in  conjunction  with  any  of  the  known  static  point  location  data  structures,  and 
allows  a  tradeoff  between  query  and  insertion  time.  Namely,  we  show  that  for  any  2  £b  <V/T, 
there  exists  a  dynamic  point  location  data  structure  with  0(Iogi,  n  logn)  query  time  and 
0((logfe  n)2  b  log  b)  update  time. 

By  setting  b- 2,  we  obtain  (9 (log2/!)  query  and  update  times,  which  match  the  result  of 
Chapter  3.  By  setting  b  =  log  n,  we  obtain  <9(log2»/loglogn)  query  time  and 
O  (log3  n  /  log  log  n)  update  time.  Finally,  by  setting  b  =  Vn ,  we  obtain  <2 (log  n)  query  time  and 
0(yfn  logn)  update  time.  This  method  is  based  on  the  dynamization  techniques  for  decompos¬ 
able  search  problems  described  in  [3, 42], 

The  rest  of  the  chapter  is  organized  as  follows.  Section  5.2  contains  preliminary 
definitions.  The  dynamic  technique  is  described  in  Section  5.3. 

5.2.  Preliminaries 

Let  91  be  a  triangulation.  We  consider  the  following  update  operations  on  91: 

INSERTSTAR  (v,r;r1,r2,r3):  Add  vertex  v  inside  region  r  and  edges  between  v  and  the 
vertices  of  r,  which  is  decomposed  into  new  regions  r  i ,  r2,  and  r3  (see  Fig.  5.1(a)). 

REMOVESTAR  (v,r1,r2,r3;r):  Remove  degree-3  vertex  v  and  its  incident  edges,  which 
merges  the  three  regions  formerly  containing  v  into  a  new  region  r  (see  Fig.  5.1(a)). 

SWAlPDIAGONAL  (c;rltr2):  Let  e  be  an  edge  such  that  the  union  of  the  two  regions  on 
the  left  and  right  of  e  is  a  convex  quadrilateral.  Remove  e  and  reinsert  it  as  the  other  diagonal  of 
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the  quadrilateral,  thus  creating  two  new  regions  rx  and  r2.  (see  Fig.  5.1(b)). 

The  problem  of  performing  point  location  in  SR  under  the  above  set  of  dynamic  updates  will 
be  called  problem  PL-1.  The  following  theorem  demonstrates  the  adequacy  of  the  above  reper¬ 
tory. 


Figure  5.1  (a)  Example  of  operation  INSERTSTAR/ 
REMOVESTAR-,  (b)  Example  of  operation 
SWAPDIAGONAL . 


f« 


Theorem  5.1  An  arbitrary  triangulation  SR  with  n  vertices  can  be  assembled  starting  from  the 
empty  subdivision,  and  disassembled  to  obtain  the  empty  subdivision,  by  a  sequence  of  O(n) 
operations  from  the  above  repertory. 

Proof:  We  constructively  demonstrate  only  the  disassembly  of  SR,  since  the  assembly  is  obtain¬ 
able  by  replacing  each  operation  with  its  inverse  and  by  reversing  their  order.  If  91  has  a  vertex 
u  of  degree  3,  we  remove  u  by  means  of  a  REMOVESTAR  operation.  Otherwise,  we  know  that 
SR  has  at  least  one  vertex  v  of  degree  de  {4,5 }  (an  easy  consequence  of  Euler’s  formula).  Con¬ 
sider  the  neighbors  v<j,v  i,  •  •  •  ,Vd-\  of  vertex  v,  in  clockwise  order  around  v.  Each  pair  of  two 
consecutive  neighbors  forms  with  v  a  triangle  /,•  =  (v,  v(-,v; + 2 )  (sum  over  the  indices  is  modulo  d). 
Since  the  sum  of  the  internal  angles  of  the  polygon  (vo.vj,  •  •  •  ,vj_i)  is  equal  to  (d - 2)  jc,  it  is 
easy  to  verify  that  for  at  least  one  index  i,  the  union  of  triangles  r,  and  r;  +  1  is  a  convex  quadrila¬ 
teral.  Hence,  by  a  sequence  of  at  most  two  SWAPDIAGONAL  and  one  REMOVESTAR  opera¬ 
tion,  we  can  remove  vertex  v.  A  simple  inductive  argument  completes  the  proof.  □ 

First,  we  consider  another  variant  of  the  dynamic  point  location  problem,  denoted  as  prob¬ 
lem  PL-2,  which  is  described  as  follows:  Let  T  be  a  set  of  triangles  such  that  the  interiors  of  any 
two  triangles  of  T  are  disjoint.  We  want  to  perform  the  following  operations  on  T: 

FIND  ( q,T ):  Return  the  triangle  t  of  T  that  contains  the  query  point  q;  if  q  is  on  some  edge 
(or  vertex)  of  t,  the  edge  (or  vertex)  itself  is  returned;  if  no  triangle  contains  q,  a  null  value  is 
returned. 

ADD  ( t,T ):  Add  triangle  f  to  T;  the  operation  is  allowed  only  if  t  does  not  intersect  the 
interior  of  any  other  other  triangle  of  T. 

DELETE  ( t,T ):  Remove  triangle  t  from  T. 

We  denote  with  Find(m),  Add{m),  and  Delete(m)  the  time  to  perform  operations  FIND,  ADD, 
and  DELETE,  respectively,  where  m  is  number  of  triangles  of  T. 

Since  a  triangulation  31  with  n  vertices  has  m  =  2  n  -  4  regions,  and  each  of  the  update 
operations  for  problem  PL-1  corresponds  to  at  most  four  ADD/DELETE  operations  in  problem 
PL-2,  we  have 
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Lemma  5.1  Given  a  data  structure  for  problem  PL-2,  there  exists  a  data  structure  for  problem 
PL-1  with  the  same  (order  of)  space  requirement  and  query  time,  and  such  that  each  update 
operation  takes  time  O  (Add(n)+ Delete  (n)). 

S3.  The  Incremental  Reconstruction  Method 

In  this  section,  we  show  an  efficient  solution  to  problem  PL-2  based  on  the  following 
straightforward  decomposition  property: 

Lemma  5.2  Denoting  by  (J\,T2)  a  partition  of  T,  we  have 

FIND  (p,T)=F!ND  ( p.Tj)  u  FIND  (p,T2)-  □ 

The  main  idea  is  to  maintain  a  partition  of  the  set  T  into  O  (log  m)  subsets,  where  each  sub¬ 
set  has  a  separate  static  data  structure  for  point  location.  Using  the  decomposition  property,  the 
answer  to  FIND  (p,T)  can  be  obtained  by  performing  a  separate  FIND  in  each  subset,  and  then 
combining  the  answers.  Since  each  subset  has  no  more  than  m  elements,  each  subquery  takes 
O(logm)  time,  so  that  the  global  time  for  operation  FIND  is  O  (log2m).  Also,  since  the  space 
requirement  for  each  subset  is  linear,  the  overall  space  requirement  is  O  (m).  In  order  to  show 
that  operation  ADD  can  also  be  efficiently  performed,  we  describe  the  data  structure  in  more 
detail. 

Let  av,av- 1* '  ’  *  be  the  binary  representation  of  m  (the  number  of  triangles  of  T ), 
where  v=  [log/nj  .  The  dynamic  data  structure  D  for  T  consists  of  v  + 1  static  point  location 
data  structures  So,  *  •  *  ,Sy,  where  S;  represents  a  subdivision  induced  by  exactly  a,-2‘  triangles. 
After  the  addition  of  a  new  triangle,  m  increases  to  m  + 1,  so  that  in  order  to  maintain  the  data 
structure  D,  we  have  to  discard  the  static  structures  5,  corresponding  to  the  bits  that  change  from 
one  to  zero,  and  create  new  static  structures  for  the  bits  that  change  from  zero  to  one.  In  the 
worst  case,  namely  when  m  + 1  is  a  power  of  two,  this  takes  O  (m)  time,  since  we  have  to  discard 
all  the  existing  static  structures  and  build  from  scratch  a  new  data  structure  for  m  + 1  triangles. 
However,  the  average  time  complexity  over  a  sequence  of  m  ADD  operations  to  an  initially 
empty  data  structure,  denoted  Add(m),  is  Add(m)  =  0 (log^m).  To  see  this,  assume  for  simpli¬ 
city  that  m  =  2V  - 1 .  From  elementary  properties  of  binary  counting,  the  static  structure  S,  is 


built  exactly  2V~‘ " 1  times,  so  that  the  total  insertion  cost  over  a  sequence  of  m  insertions  is 

mA55(m)=Vj;12v_i'1P(2‘) 

i=0 


where  P(k)  =  0(klogk)  is  the  (preprocessing)  cost  for  building  an  optimal  static  point  location 
structure  fori  triangles  [13].  With  simple  manipulations,  we  obtain 

=  0(mlog2m). 


fv-l  .  .  .1 

•  * 
,v-l 

2 V_1  £i 

m  Add(m)=0 

M 

7 

=  0 

l«=o  j 

i  =  o. 

At  this  point,  we  have  an  0(m)  space  semidynamic  data  structure  with  0{m)  space 
requirement  that  supports  operation  FIND  in  time  0(log2m)  (worst-case)  and  operation  ADD  in 
time  Oilog^m)  (amortized).  The  amortized  time  bound  for  ADD  can  be  turned  into  a  “worst- 
case”  bound  by  spreading  the  work  for  building  each  static  structure  Si  over  several  consecutive 
steps,  so  that  after  each  ADD  operation  0  (log 2m)  work  is  performed  [3, 42].  This  data  structure 
has  a  special  type  of  worst-case  time  bound,  since  the  restructuring  operations  are  distributed 
over  time.  We  will  refer  to  such  time  bound  as  a  distributed  worst-case  time. 

The  above  discussion  can  be  summarized  as  follows: 


Lemma  S3  There  is  a  data  structure  for  problem  PL-2  with  0{m)  space  requirement  that  sup¬ 
ports  operations  FIND  and  ADD  in  time  O  (log2m)  (distributed  worst-case). 

In  order  to  perform  efficiently  operation  DELETE  also,  we  use  a  “lazy”  deletion  scheme. 
First,  we  modify  the  aforementioned  data  structure  by  adding  a  mark  field  to  every  record  that 
describes  a  triangle  of  T,  so  that  DELETE  ( t,T )  can  be  simply  performed  by  setting  the  mark  of  t. 
Also,  we  modify  the  procedure  for  constructing  the  static  point  location  structures,  so  that  only 
the  unmarked  triangles  are  considered  in  search  operations,  while  the  marked  ones  serve  the  only 
purpose  of  maintaining  the  nominal  size  of  the  structure. 

This  modified  structure  D  allows  for  efficient  queries  and  updates  as  long  as  the  number  k 
of  marked  triangles  does  not  exceed  a  fixed  fraction  of  the  total  number  m  of  triangles,  say  one 
half.  If  after  a  deletion  we  have  k>m/ 2,  we  start  building  a  new  structure  D'  of  the  same  type 
for  the  unmarked  triangles  only.  We  spread  the  construction  of  D'  over  m/4  steps,  where  a  step 
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corresponds  to  an  update  operation.  At  each  step  we  add  two  “old”  triangles  to  D'  and  we  per¬ 
form  the  current  update  operation  ( ADD  or  DELETE)  on  both  D  and  D'.  After  the  completion  of 
all  the  steps,  D'  is  ready  and  replaces  D.  Notice  that  the  total  number  m  of  triangles  in  D’  and 
the  number  of  marked  triangles  X  satisfy  the  relations 

m  .  ,.m  ,  m  3  n  sis* m 
2  2  4  4  4 

so  that  we  have  Elm  £  1/2,  which  means  that  also  D'  can  support  queries  and  updates  efficiently. 
During  the  construction  of  D',  we  perform  the  FIND  operations  on  D,  for  which  the  number  of 
marked  triangles  does  not  exceed  one  fourth  of  the  total  number  of  triangles. 

Therefore,  by  using  at  most  two  data  structures  that  allow  only  for  insertions  and  lazy  dele¬ 
tions,  and  maintaining  the  invariant  that  the  number  of  marked  triangles  does  not  exceed  a  fixed 
fraction  of  the  total  number  of  triangles,  we  obtain  a  complete  dynamic  solution  for  problem 
PL-2  with  O  (log2m)  query  and  update  time: 

Theorem  5.2  There  is  a  data  structure  for  problem  PL-2  with  O  (m)  space  requirement  that  sup¬ 
ports  operations  FIND,  ADD,  and  DELETE  in  time  O  (log2m)  (distributed  worst-case). 

Recalling  Lemma  5.1,  we  obtain  the  following  result  for  the  original  problem  PL-1. 

Corollary  5.1  There  is  a  data  structure  for  the  dynamic  point  location  problem  PL-1  with  O  (n) 
space  requirement  and  O  (log2n)  (distributed  worst-case)  query  and  update  times. 

The  aforementioned  technique  for  problem  PL-2  makes  use  of  the  properties  of  counting  in 
base  2.  By  using  a  larger  base,  we  can  reduce  the  number  of  static  data  structures  used.  In  fact, 
using  a  base  b  St  2  the  above  data  structure  can  be  modified  to  work  with  a  partition  of  T  into  sub¬ 
sets  Sij,  i  =0,  •  •  •  ,v;  j  =  1,  ■  •  •  ,(b  - 1);  where  v  =  log*,  m,  and  each  subset  has  size  b‘. 

Using  arguments  similar  to  the  ones  already  given  for  the  case  b  =  2,  we  can  show  that  this 
data  structure  has  O  (m)  space  requirement  and 

Find(m) =0  (log*  m  log  m);  Add(m)  =  Delete (m)  =  0((log*  m)2  b  log  b). 

To  obtain  an  improvement  in  the  search  time,  we  can  let  the  base  b  grow  with  m.  Applying  the 


results  described  in  [38],  we  obtain 


Theorem  5.3  Let  91  be  a  triangulation  with  n  vertices,  and  let  b -b{n)  be  a  smooth  nondecreas¬ 
ing  integer  function  with  2 There  exists  a  data  structure  for  point  location  problem 
PL-1  with  O(n)  space  requirement,  O  (log*  n  logn)  query  time  and  0^(log*,  n)2  b  log&j  update 
time. 
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CHAPTER  6 


DYNAMIC  TRANSITIVE  CLOSURE 


6.1.  Introduction 

The  notion  of  a  planar  sf-graph  -  i.e.,  a  planar  acyclic  digraph  embedded  in  the  plane  with 
exactly  one  source,  s,  and  one  sink,  t,  both  on  the  external  face  -  was  first  introduced  in  the 
planarity  testing  algorithm  of  Lempel  et  al.  [36],  and  was  fruitfully  used  in  a  number  of  applica¬ 
tions,  which  include  planar  graph  embedding  [9,30,54],  graph  planarization  [29,44],  graph 
drawing  algorithms  [10,53,59],  floor  planning  [4,57]  planar  point  location  [13,35],  visibility 
representations  [41,49,51,58],  motion  planning  [21,48],  and  VLSI  layout  compaction  [25,57]. 
Also,  planar  sr-graphs  are  important  in  the  theory  of  partially  ordered  sets  since  they  are  associ¬ 
ated  with  planar  lattices  [32]. 

In  this  chapter  we  further  the  investigation  of  these  structures,  and  show  that  any  planar  st- 
graph  G  admits  two  total  orders  (referred  to  as  leftist  and  rightist  orders)  on  the  set  VkjE  uf, 
where  V,  E,  and  F  are  respectively  the  set  of  vertices,  edges,  and  faces  of  G.  Each  of  these  two 
orders  yields  a  unique  representation  of  G  as  a  string  whose  terms  are  symbols  representing  all 
the  topological  constituents  of  G.  Graph  G  can  be  dynamically  modified  by  means  of  insertion 
of  edges  and  expansions  of  vertices,  and  of  their  inverses.  These  operations  form  a  complete  set, 
since  any  n-vertex  planar  sr-graph  can  be  assembled  or  disassembled  by  an  appropriate  sequence 
of  O(n)  such  operations. 

The  central  result  of  this  chapter  is  that  the  string  representation  of  the  graph  resulting  from 
one  of  the  postulated  updating  operations  is  obtained  as  a  syntactic  transformation  of  the  original 
string  representation.  This  transformation  consists  of  the  execution  of  O  (1)  primitives,  such  as 
insertions,  deletions,  and  swaps  of  substrings. 

This  general  framework  provides  the  theoretical  underpinning  and  unifying  viewpoint  for 
three  significant  applications:  point  location  in  a  planar  monotone  subdivision,  transitive-closure 
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query  in  planar  sr-graphs,  and  contact-chain  query  in  convex  subdivisions.  In  this  chapter  we 
shall  only  briefly  illustrate  (in  Section  6.4)  the  connection  between  planar  st- graphs  and  mono¬ 
tone  subdivisions,  since  the  point  location  problem  in  the  latter  has  been  treated  earlier  in 
Chapter  3  in  purely  geometric  terms. 

A  transitive-closure  query  for  a  planar  sr- graph  G  consists  of  testing  for  the  existence  of 
(and/or  reporting)  a  directed  path  between  two  vertices  u  and  v  of  G.  We  are  interested  in  a 
graph  G  that  can  be  dynamically  modified. 

The  previous  best  results  concern  semi-dynamic  versions  of  the  transitive-closure  query 
problem  in  digraphs  (namely,  edge  insertions  in  general  digraphs  and  edge  deletions  in  acyclic 
digraphs),  and  have  0(1)  query  time,  O  (n)  amortized  update  time,  and  0(n2)  storage 
[27, 28, 45].  In  this  chapter  we  establish  the  following  result: 

Theorem  6.1  Let  G  be  a  planar  sr-graph  with  n  vertices.  There  exists  an  O  (n)-space  dynamic 
data  structure  for  the  transitive-closure  query  problem  on  G,  which  supports  queries  and  updates 
in  time  O  (log  n)  (worst-case). 

Finally,  we  consider  the  problem  of  contact-chain  query  in  convex  subdivisions,  which 
arises  in  motion  planning  and  computer  graphics,  and  is  described  as  follows  [8, 21, 48].  Given  a 
convex  subdivision  91  of  the  plane  and  an  (oriented)  direction  9,  we  say  that  region  rj  pushes 
region  r2  if  rx  and  r2  are  adjacent  and  there  exists  a  line  in  direction  0  which  intersects  rj  and 
r2  in  that  order.  A  contact  chain  in  91  is  a  sequence  of  regions  r1(r2,  such  that  r,  pushes 

ri  +  i  for  i  =  1,  •  •  •  ,k  - 1  (see  Fig.  6.1).  Assume  that  the  regions  of  91  are  rigid  objects,  and  we 
want  to  translate  them  one  at  a  time  in  direction  0  avoiding  collisions.  Then  the  existence  of  a 
contact  chain  from  r  j  to  r2  implies  that  r2  obstructs  rj,  i.e.,  r2  must  be  translated  before  r1. 

A  contact-chain  query  consists  of  testing  the  existence  of  (and/or  reporting)  a  contact  chain 
between  two  regions  of  91.  We  are  interested  in  answering  contact-chain  queries  in  a  very  gen¬ 
eral  dynamic  environment,  where  91  can  be  updated  by  means  of  insertions/deletions  of  vertices 
and  edges,  and  the  direction  0  can  be  changed  by  elementary  increments/decrements.  (An  ele¬ 
mentary  increment/decrement  of  direction  is  such  that  the  push  relation  is  inverted  for  exactly 
one  pair  of  adjacent  regions.)  Casting  this  problem  in  the  planar  jr-graph  framework,  we 
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Figure  6.1  Example  of  contact  chain. 


establish  the  following  result: 


Theorem  6.2  Let  31  be  a  convex  subdivision  with  n  vertices.  There  exists  an  <3(n)-space 
dynamic  data  structure  for  the  contact-chain  query  problem  in  91,  which  supports  queries  and 
updates  in  time  O (log n)  ( worst-case). 

The  rest  of  this  chapter  is  organized  as  follows.  Section  6.2  provides  preliminary 
definitions  and  properties  of  planar  sr-graphs.  Section  6.3  presents  the  technique  for  the  dynamic 
maintenance  of  planar  sr-graphs.  Applications  to  planar  point  location,  transitive  closure,  and 
contact  chains  are  described  in  Section  6.4. 

6.2.  Planar  s/*Graphs 

Basic  definitions  on  graphs  and  posets  can  be  found  in  textbooks  such  as  [5, 15], 


Let  G  be  a  directed  graph,  for  brevity  digraph,  and  v  a  vertex  of  G.  As  usual,  we  denote  by 
deg-(v)  the  indegree  of  v,  i.e.,  the  number  of  incoming  edges  of  v,  and  by  deg+(v)  the  outdegree 
of  v,  i.e.,  the  number  of  outgoing  edges  of  v.  A  source  of  G  is  vertex  s  with  deg~(r) =0.  A  sink 
of  G  is  vertex  t  with  deg+(r) = 0.  A  transitive  edge  of  G  is  an  edge  e  =  (u, v)  such  that  there  exists 
another  directed  path  from  u  to  v  consisting  of  at  least  two  edges. 

A  planar  st- graph  is  a  planar  acyclic  digraph  G  with  exactly  one  source,  s,  and  exactly  one 
sink,  r,  which  is  embedded  in  the  plane  so  that  s  and  t  are  on  the  boundary  of  the  external  face 
(see  Fig.  6.2). 

These  graphs  were  first  introduced  in  the  planarity  testing  algorithm  of  Lempel  et  al.  [36]. 
Several  important  properties  of  planar  rr-graphs  are  expressed  by  the  following  lemmas: 


v0  =  s 

Figure  6.2  Example  of  planar  j/- graph. 
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Lemma  6.1  [36]  Every  vertex  of  G  is  on  some  directed  path  from  s  to  t. 

Lemma  6.2  [51]  For  every  vertex  v  of  G ,  the  incoming  (outgoing)  edges  appear  consecutively 
around  v  (See  Fig.  6.3(a)). 

Lemma  6.3  [51]  For  every  face  /  of  G ,  the  boundary  of  /consists  of  two  directed  paths  with 
common  origin  and  destination.  (See  Fig.  6.3(b)). 

Lemma  6.4  [10, 33]  G  admits  a  planar  upward  drawing,  i.e.,  a  planar  drawing  such  that  every 
edge  (u,v)  is  a  curve  monotonically  increasing  in  the  vertical  direction. 

Let  P  be  a  poset  (partially  ordered  set),  where  c  denotes  the  partial  order  on  the  elements 
of  P.  The  Hasse  diagram  (also  called  covering  digraph)  of  P  is  a  digraph  G  whose  vertices  are 
the  elements  of  P,  and  such  that  (u,v)  is  an  edge  of  G  if  and  only  if  u«v  and  there  is  no  other 
element  x  of  P  such  that  «<tct<cv.  G  is  acyclic  and  has  no  transitive  edges,  i.e.,  G  is  the 


HIGH  (/) 


LOW  {f) 

(a)  (b) 

Figure  6.3  (a)  Example  for  Lemma  6.2.  (b)  Example  for  Lemma  6.3. 
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transitive  reduction  of  the  graph  of  <c  on  vertex  set  P.  Hasse  diagrams  are  usually  represented 
by  straight-line  drawings  such  that  for  each  edge  (u,v)  the  ordinate  of  vertex  u  is  smaller  that 
that  of  vertex  v. 

A  planar  lattice  is  a  poset  whose  Hasse  diagram  is  a  planar  sr-graph.  Also,  every  plane  st- 
graph  without  transitive  arcs  is  the  Hasse  diagram  of  some  planar  lattice.  Several  properties  of 
planar  lattices  are  described  in  [32]. 

A  linear  extension  of  a  poset  P  is  a  total  order  <  on  the  elements  of  P  such  that,  for  any 
two  elements  u  and  v  of  P,  u<v  implies  u  <  v.  A  linear  extension  corresponds  to  a  topological 
sorting  of  the  vertices  of  the  Hasse  diagram  of  P.  We  say  that  P  has  dimension  k  if  G  admits  k 
linear  extensions  <1 ,  <2  ,*•*><*,  such  that  ucv  if  and  only  if  u  <\  v  ,  u  <2  v  ,  •  •  • ,  u  <*  v, 
and  k  is  minimum. 

It  is  known  that  planar  lattices  have  dimension  2  [5,  p.  32,  ex.  7(c)]  [31, 32],  which  implies 
the  following  lemma: 

Lemma  6.5  [5,31,32]  Let  G  be  a  planar  st- graph  with  n  vertices.  There  exist  two  total  orders 
on  the  vertices  of  G,  denoted  <1  and  <r  ,  such  that  there  is  a  directed  path  from  u  to  v  if  and 
only  if  u  <1  v  and  u  <R  v.  Furthermore,  orders  <1  and  <r  can  be  computed  in  O  (n)  time. 

Lemma  6.5  is  based  on  the  fact  that  the  underlying  partial  order  of  a  planar  lattice  admits  a 
“complementary”  partial  order  (see  [32]).  Figure  6.4(a)  shows  a  planar  sr-graph  where  each 
vertex  is  labeled  by  its  ranks  in  the  orders  <1  and  <r  . 

In  the  following  definitions,  the  concepts  of  left  and  right  refer  to  the  orientation  of  the 
edges.  For  example,  the  face  to  the  left  of  an  edge  («,v)  is  the  face  containing  edge  e  which 
appears  on  the  left  side  when  traversing  edge  («,v)  from  vertex  u  to  vertex  v.  Also,  the  reader 
will  find  it  convenient  to  visualize  the  planar  sr-graph  G  as  being  drawn  in  the  plane  with  edges 
monotonically  increasing  in  the  vertical  direction  (see  Lemma  6.4). 

Given  vertices  u  and  v  of  G  such  that  there  exists  a  path  from  u  to  v,  the  set  of  paths  from  u 
to  v  defines  a  planar  sr-graph  with  source  u  and  sink  v  which  is  an  induced  subgraph  of  G.  The 
two  paths  that  form  the  external  boundary  of  this  subgraph  will  be  called  the  leftmost  path  and 
rightmost  path  from  u  to  v,  respectively.  For  example,  the  external  boundary  of  G  consists  of  the 
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(a)  (b) 


|  Figure  6.4  (a)  Orders  <i  and  <r  on  the  vertices  of  a 

planar  w-graph;  (b)  A  planar  st- graph  G  and  its  dual  G*. 

leftmost  and  rightmost  paths  from  s  to  t. 

Let  G*  be  the  digraph  obtained  from  the  dual  graph  of  G  as  follows  (see  Fig.  6.4(b)): 
(1)  the  dual  edge  e*  of  an  edge  e  is  directed  from  the  face  to  the  left  of  e  to  the  face  to  the  right 
of  e;  (2)  the  external  face  of  G  is  dualized  to  two  vertices  of  G*,  denoted  s *  and  t* ,  which  are 
_  incident  with  the  duals  of  the  edges  on  the  leftmost  and  rightmost  paths  from  s  to  r,  respectively. 

Vertices  s*  and  t*  can  be  thought  of  as  being  the  “left”  and  “right  external  face”  of  G,  respec¬ 
tively.  It  is  simple  to  verify  that  G*  is  a  planar  jr-graph  with  source  s*  and  sink  t *  [41,51]. 
.  Notice  that  G  *  might  have  multiple  arcs. 

Let  V ,  E,  and  F  denote  the  set  of  vertices,  edges,  and  faces  of  G,  respectively,  where  F  has 
elements  s*  and  t*  representing  the  external  face.  We  will  show  that  the  orcers  <i  and  <r  can 
be  extended  to  the  set  V^jEkjF,  thereby  giving  a  unique  total  order  of  all  topological  consti¬ 
tuents  of  G. 
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First,  for  each  element  xofVuEuF,  we  define  vertices  LOW (x)  and  HIGH (x),  and  faces 
LEFTix)  and  RIGHT (x),  as  follows: 

(1)  If  x=ve  V,  we  define  LOW (v)= HIGH (v)=v.  Also,  with  reference  to  Lemma  6.2  and 
Fig.  6.3(a),  we  denote  by  LEFT(v )  and  RIGHT (v)  the  two  faces  that  separate  the  incoming 
and  outgoing  edges  of  a  vertex  v*s,r,  where  LEFT  (v)  is  the  face  to  the  left  of  the  leftmost 
incoming  and  outgoing  edges,  and  RIGH1  (v)  is  the  face  to  the  right  of  the  rightmost 
incoming  and  outgoing  edges.  For  v  =  5  or  v  =r,  we  conventionally  define  LEFT  ly)-s* 
and  RIGHT (v)  =  t\ 

(2)  If  x=e  e  £,  we  define  LOW (e)  and  HIGH ( e )  as  the  tail  and  head  vertices  of  e,  respec¬ 
tively.  Also,  we  denote  by  LEFT  (e)  and  RIGHT  {e)  the  faces  oc  the  left  and  right  side  of  e, 
respectively. 

(3)  If  x  =/  e  F  and  f*s  * , ;  * ,  we  denote  by  LOW  (f )  and  HIGH  (f )  the  two  vertices  that  are  the 
common  origin  and  destination  of  the  two  paths  forming  the  boundary  off  (see  Lemma  6.3 
and  Fig.  6.3(b)).  Vertices  LOW  ( / )  and  HIGH  if)  are  called  the  extreme  vertices  of  face  /. 
For  f=s*  or  f=t*,  LOW(f)  and  HIGHif)  are  undefined.  Also,  we  define 
LEFT(f)=RIGHT(f)=f.  Finally,  the  two  directed  paths  forming  the  boundary  of  /arc 
called  the  left  path  and  right  path  of  /,  respectively. 

We  say  that  x  is  below  y,  denoted  xty,  if  there  is  a  path  in  G  from  HIGH (x)  to  LOW(y). 
Also,  we  say  that  x  is  to  the  left  ofy,  denoted  x-*y,  if  there  is  a  path  in  G*  from  RIGHT (x)  to 
LEFTiy). 

For  example,  in  the  planar  tt-graph  shown  in  Fig.  6.2,  we  have  £2Tv4,/4  tv4,  v5-»/4,  and 

Lemma  6.6  Relations  T  and  — >  are  partial  orders  on  VkjEkjF. 

Proof:  A  consequence  of  the  fact  the  graphs  G  and  G*  are  acyclic.  □ 

The  following  lemma  shows  that  T  and  — »  are  complementary  partial  orders. 


76 


Lemma  6.7  Let  x  and  y  be  any  two  elements  of  V  uE  uF.  Then  one  and  only  one  of  the  fol¬ 
lowing  holds: 

xty,  ylx,  x-*y,  y-*x. 

Proof:  We  prove  the  theorem  for  the  case  when  y  =  v  is  a  vertex  of  G.  The  other  two  cases  can 
be  proved  using  similar  arguments. 

Let  7ti  and  JC2  be  the  leftmost  and  rightmost  paths  from  s  to  v,  respectively.  Also,  let  713 
and  7C4  be  the  leftmost  and  rightmost  paths  from  v  to  t,  respectively.  These  paths  partition 
VuEuF  into  five  subsets,  one  of  which  is  v,  and  the  others  are  defined  as  follows  (see  Fig. 
6.5): 

(1)  A  contains  the  vertices,  edges,  and  faces  enclosed  by  paths  7ti  and  7t2,  including  the  vertices 
and  edges  of  these  paths,  but  excluding  v; 


Figure  6  J  Partiton  ofVu£uF  with  respect  to  vertex  v. 
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(2)  B  contains  the  vertices,  edges,  and  faces  enclosed  by  paths  713  and  including  the  vertices 
and  edges  of  these  paths,  but  excluding  v; 

(3)  C  contains  the  vertices,  edges,  and  faces  to  the  left  of  paths  7Cj  and  713 ,  excluding  the  ver¬ 
tices  and  edges  of  these  paths; 

(4)  D  contains  the  vertices,  edges,  and  faces  to  the  right  of  paths  ard  714,  excluding  the  ver¬ 
tices  and  edges  of  these  paths. 

It  is  easy  to  verify  that  the  edges  of  A  are  those  of  a  planar  sr-graph  with  source  s  and  sink 
v,  which  is  an  induced  subgraph  of  G,  and,  similarly,  the  edges  of  B  are  those  of  a  planar  st- 
graph  with  source  v  and  sink  t.  Notice  that  the  vertices  w  of  A  are  exactly  those  such  that  there  is 
a  directed  path  in  G  from  w  to  v,  and  analogously  for  the  vertices  of  B. 

Using  simple  duality  arguments,  we  can  show  that  the  duals  of  the  edges  of  C  are  those  of  a 
planar  rr-graph  with  source  s*  and  sink  LEFT (v),  which  is  an  induced  subgraph  of  G*.  Simi¬ 
larly,  the  duals  of  the  edges  of  D  are  those  of  a  planar  sf-graph  with  source  RIGHT  (v)  and  sink 
t*.  Notice  that  the  faces  /of  C  are  exactly  those  such  that  there  is  a  directed  path  in  G*  from / 
to  LEFT  (v),  and  analogously  for  the  faces  of  D. 

By  Lemma  6.1,  there  are  directed  paths  from  every  vertex  of  A  to  v,  and  from  v  to  every 
vertex  of  B.  Since  for  every  edge  or  face  x  of  A  ( B ),  both  LOW (x)  and  HIGH(x)  are  in  A  (B), 
we  conclude  that  x  e  A  implies  xtv  and  x  e  B  implies  vtx.  With  similar  arguments,  we  con¬ 
clude  that  x  e  C  implies  x— >v  and  x  e  D  implies  v—*x. 

It  remains  to  be  shown  that  relations  T  and  ->  are  mutually  exclusive.  Let  xe  A  uB,  i.e., 
either  xTv  or  vTjc.  Suppose  xtv;  if  x— >v,  then  there  is  a  path  in  G*  from  RIGHT (x)  to 
LEFT{v).  This  implies  that  RIGHT  (x)  e  C,  a  contradiction.  An  analogous  contradiction  is 
reached  if  we  assume  that  xtv  and  v— »x  jointly  hold.  Finally,  let  x  e  C  uD,  i.e.,  either  x— »v  or 
v— >x.  Suppose  x— >v;  if  x Tv, then  there  is  a  path  in  G  from  HIGH(x)  to  v.  This  implies  that 
HIGH  (x)g  A,  a  contradiction.  An  analogous  contradiction  is  reached  if  we  assume  that  x-»v 
and  vtx  jointly  hold.  □ 

We  define  relations  <1  and  <r  on  VkjE  uF,  as  follows: 
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x  <Ly  «*xTy  or x->y\  x  <*y  <=>xty  ory-wc. 


r 


As  a  consequence  of  Lemma  6.7,  we  obtain 
Theorem  6.3  The  relations  <l  and  <r  on  V  u£  uF  are  total  orders. 

We  also  note  that  there  is  a  path  in  G  from  vertex  u  to  vertex  v  if  and  only  if  u  <l  v  and 
u  <r  v,  since  such  path  exists  if  and  only  if  «Tv. 

We  define  the  left-sequence  of  G  as  the  sequence  of  elements  of  VuE  uF,  sorted  accord¬ 
ing  to  <i  (leftist  order).  The  right-sequence  of  G  is  defined  similarly  with  respect  to  <r 
(rightist  order). 

For  example,  the  right-sequence  of  the  graph  of  Fig.  6.2  is: 

f5Voe3Ue2Vie4V2eiif2esv3e1V4elof2e6fle1v5e9v6f0. 

We  will  use  a  convenient  string  notation  for  such  sequences.  Namely,  we  use  terminal 
symbols  (lower-case  letters)  for  the  elements  of  VuE  uF,  and  variables  (upper-case  letters)  for 
substrings  of  the  left-  or  right-sequence.  For  example,  the  left-sequence  of  the  graph  of  Fig.  6.2 
can  be  represented  by  the  string 

foVoeiAv3e6v5e9f2B 

where  A  =f\e2v\e5  w&B  =e-tf3e *f  \e3v2e%v le^v (J$. 

6.3.  On-Line  Maintenance  of  a  Planar  sf-Graph 

In  this  section  we  define  a  complete  set  of  update  operations  on  a  planar  sr-graph,  and  show 
that  the  restructuring  of  the  orders  <i  and  <#  resulting  from  any  such  update  operation  can  be 
expressed  by  means  of  a  simple  string  transformation.  From  this  result,  we  derive  an  efficient 
data  structure  for  the  on-line  maintenance  of  the  two  orders  of  a  planar  tf-graph. 

The  update  operations  on  a  planar  st-  graph  are  defined  as  follows: 

INSERTEDGE  (e,u,v,f;fl,f2):  Add  edge  e  =(n,v)  inside  face  /,  which  is  decomposed 
into  faces  f\  and  f2,  with  f\  to  the  left  of  e  and  f2  to  the  right  (see  Fig.  6.6(a)). 
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REMOVEEDGE  (e,u,v,fi  ,/2;/) :  Delete  edge  e  =  (u,v)  and  merge  the  two  faces  f\  and  /2 
formerly  on  the  two  sides  of  e  into  a  new  face  /  (see  Fig.  6.6(a)). 

EXPANDVERTEX  (e,/,g,v;v i,V2):  Expand  vertex  v  into  vertices  vj  and  v2,  which  are 
connected  by  a  new  edge  e  with  face  /  to  its  left  and  face  g  to  its  right  (see  Fig.  6.6(b)). 

CONTRACTVERTEX  (e,f,g,v i,v2;v):  Contract  edge  e  =(v1,v2),  and  merge  its  endpoints 
into  a  new  vertex  v.  Faces  /  and  g  are  to  the  left  and  right  of  e,  respectively  (see  Fig.  6.6(b)). 
Parallel  edges  resulting  from  the  contraction  are  merged  into  a  simple  edge. 

Each  operation  is  allowed  if  the  resulting  graph  is  itself  a  planar  rr-graph.  It  is  interesting 
to  observe  that  operations  EXPANDVERTEX  and  CONTRACTVERTEX  are  dual  of  INSERT- 
EDGE  and  REMOVEEDGE ,  respectively,  since  performing  one  on  G  corresponds  to  performing 
the  other  on  G*. 

We  say  that  an  edge  e  of  G  is  removable,  if  operation  REMOVEEDGE  (e,u,v,fi  ,/2/)  on  G 
yields  a  planar  rr-graph.  We  say  that  e  is  contractible  if  operation 
CONTRACTVERTEX  ( e,f,g,v  i ,v2;v)  on  G  yields  a  planar  sr-graph. 

Lemma  6.8  Let  e  be  an  edge  of  G. 

(1)  If  e  is  not  removable  then  it  is  contractible. 

(2)  If  e  is  not  contractible  then  it  is  removable. 

Proof:  From  the  definition  of  planar  sr-graph,  it  is  easy  to  see  that  an  edge  e-{u,v)  is  remov¬ 
able  if  and  only  if  deg+(«)  £  2  and  deg_(v)  £  2,  and  it  is  contractible  if  and  only  if  it  is  not  a  tran¬ 
sitive  edge.  (1)  Assume  that  edge  e  =  (m,v)  is  not  removable.  Then  we  have  deg+(«)  =  1  and/or 
deg“(v)  =  1.  This  implies  that  there  is  no  other  path  in  G  from  u  to  v,  so  that  e  cannot  be  a  transi¬ 
tive  edge.  Hence,  edge  e  is  contractible.  (2)  Conversely,  assume  that  edge  e  =(u,v)  is  not  con¬ 
tractible.  Then  e  is  a  transitive  edge,  which  implies  deg+(w)>2  and  deg_(v)£2,  so  that  e  is 
removable.  □ 

A  simple  induction  based  on  Lemma  6.8  yields 
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(a) 


Figure  6.6  (a)  Operations  INSERTEDGE  and 
REMOVEEDGE.  (b)  Operations  EXPANDVERTEX  and 
CONTRA  CTVERTEX. 
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Lemma  6.9  Let  Go  be  the  trivial  planar  st- graph  consisting  of  a  single  vertex.  Any  planar  st- 
graph  with  n  vertices  can  be  assembled  starting  form  Go  by  means  of  O(n)  INSERTEDGE  and 
EXPANDVERTEX  operations,  and  can  be  disassembled  to  yield  Gq  by  means  of  O(n)  REMO¬ 
VEEDGE  and  CONTRA  CTVERTEX  operations. 

Now,  we  describe  the  transformation  of  the  leftist  order  <i  as  a  consequence  of  operations 
INSERTEDGE  2).  Similar  arguments  hold  for  the  order  <*  and  for  operation 

EXPANDVERTEX  (e,/,g,v;vltv2). 


I 
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Theorem  6.4  Let  G  be  a  planar  rt-graph,  and  G'  be  the  graph  obtained  from  G  after  the  execu¬ 
tion  of  operation  INSERTEDGE  Depending  on  the  relative  orders  of  u,  v,  and  / 

we  have  the  following  transformations  ( left-sequence  of  G)  =>  ( left-sequence  of  Gy. 

(1)  u<iv<Lf:  AuB  v  C  fD  =>  AuB  f\ev  C  fiD\ 

(2)  f<LU<Lv  '  AfBuCvD  =>  A  f\  B  ue  /2C  v  D; 

(3)  «<l/<lv:  AuBfCvD  =*  AuBfiefiCvD; 

(4)  v<Lf<Lu :  Av  B  f  C  uD  =>  Af\C  uev  B  f2D. 

Proof:  The  four  cases  are  illustrated  in  Fig.  6.7.  First,  we  observe  that  the  union  of  the  elements 
of  VkjE  uF  associated  with  any  one  of  the  substrings  A,  B,  C,  and  D,  is  a  topologically  con¬ 
nected  region  of  the  plane.  The  above  regions,  together  with  u,  v,  and  /,  form  a  partition  of  the 
entire  plane,  which  is  determined  by  four  paths,  and  specifically 

(i)  the  leftmost  path  from  HIGH  ( f)  to  r, 

07)  the  rightmost  path  from  s  to  LOW  (f );  and 

Oh)  depending  respectively  on  each  of  the  four  cases,  the  following  two  paths: 

(1)  the  leftmost  paths  from  u  to  t  and  from  v  to  t  (see  Fig.  6.7(a-b)); 

(2)  the  rightmost  paths  from  s  to  u  and  from  s  to  v  (see  Fig.  6.7 (c-d)); 

(3)  the  leftmost  path  from  u  to  t  and  the  rightmost  path  from  s  to  v  (see  Fig.  6.7(e-f)); 

(4)  the  leftmost  path  from  v  to  t  and  the  rightmost  path  from  s  to  u  (>ee  Fig.  6.7(g-h)). 

We  discuss  in  detail  Case  4  (see  Fig.  6.7  (g-h)).  The  proof  for  the  other  cases  can  be  derived  with 
similar  arguments.  The  insertion  of  edge  e  causes  every  vertex  in  C  to  be  connected  with  a 
directed  path  to  every  vertex  of  B.  At  the  same  time,  the  insertion  of  e  breaks  all  the  paths  of  G* 
from  the  faces  of  B  to  the  faces  of  C.  Hence,  we  have  the  following  relations: 

/i-»C,  Ctw,  wte,  etv,  vt B,  B->f2,  fi<iD, 

where  a  substring  represents  compactly  all  of  its  elements.  These  relations  yield  immediately 
the  updated  left- sequence.  □ 
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Theorem  6.4  shows  that  the  update  of  the  order  <i  is  a  simple  syntactic  transformation  of 
the  left-sequence,  consisting  of  at  most  four  insertions/deletions  of  elements,  and  at  most  one 
swap  of  substrings.  S;nce  operation  REMOVEEDGE  is  the  inverse  of  operation  INSERTEDGE, 
the  order  before  and  after  the  deletion  can  be  obtained  by  reversing  the  transformations  given  in 
Theorem  6.4.  The  same  situation  arises  with  respect  to  operations  EXPANDVERTEX  and  CON¬ 
TRACT VERTEX.  We  can  summarize  these  results  as  follows: 

Theorem  6.5  Let  G  be  a  planar  .rr-graph,  and  G'  be  the  graph  obtained  from  G  after  update  n, 
where  II  is  one  of  INSERTEDGE,  REMOVEEDGE,  EXPANDVERTEX,  or  CONTRA CTVERTEX 
operations.  Then  the  left-sequence  of  G'  can  be  obtained  from  the  left-sequence  of  G  by  means 
of  at  most  four  insertions/deletions  of  elements,  and  at  most  one  swap  of  substrings. 

Theorem  6.5  allows  us  to  design  a  simple  and  yet  efficient  data  structure  for  maintaining 
on-line  the  orders  of  a  planar  sr-graph  G.  We  represent  orders  <i  and  <R  by  means  of  two  bal¬ 
anced  binary  trees  (such  as  red-black  trees  [55,  pp.  52-53]),  denoted  7/,  and  TR,  where  the  left- 
to-right  order  of  the  leaves  of  Tl  gives  the  left-sequence  of  G,  and  the  left-to-right  order  of  the 
leaves  of  TR  gives  the  right-sequence  of  G.  From  Euler’s  formula,  trees  TL  and  TR  have  O(n) 
nodes,  so  that  their  depth  is  O  (log  n). 

An  order-query  on  a  planar  sr-graph  G  consists  of  determining,  given  elements  x  and  y  of 
VuEuF,  whether x  <^y  oxy  <i x,  and  similarly  with  respect  to  order  < R . 

Lemma  6.10  Let  T  be  a  binary  tree,  Xt  and  two  leaves  of  T,  and  |i  the  lowest  common  ances¬ 
tor  of  and  X2.  Then  Xi  precedes  X2  in  the  left-to-right  order  if  and  only  if  Xj  is  in  the  left 
subtree  of  p.  (and  X\  is  in  the  right  subtree  of  n). 

Proof:  The  theorem  is  an  immediate  consequence  of  the  following  recursive  definition  of  the 
left-to-right  order  of  the  leaves  of  T: 

(1)  If  leaf  Xi  is  in  the  left  subtree  of  T  and  X2  is  the  right  subtree,  then  precedes  X2. 

(2)  If  Xj  and  X2  are  in  the  same  subtree  T'  of  T,  say  the  left  subtree,  then  X\  precedes  X2  in  T  if 

and  only  if  Xi  precedes  X2  in  T'.  □ 
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Theorem  6.6  An  order-query  can  be  executed  in  O  (log  n)  time. 

Proof:  The  order-query  algorithm  is  as  follows.  We  access  the  leaves  of  tree  Ti  associated  with 
elements  x  and  y,  and  we  trace  the  paths  px  and  py  from  these  leaves  to  the  root  of  T £,.  This 
allows  us  to  find  the  lowest  common  ancestor  |i  of  leaves  x  and  y.  From  Lemma  6.10,  we  have 
that  x<iy  if  and  only  if  the  node  of  px  immediately  preceding  (j.  is  the  left  child  of  p..  Since 
paths  px  and  py  have  length  O  (log  n),  we  obtain  the  stated  time  bound.  □ 

Let  T  be  a  balanced  binary  tree.  The  left-to-right  sequence  of  the  leaves  of  T  will  be 
denoted  by  A(T).  Two  basic  operations  on  balanced  binary  trees  are  defined  as  follows: 

SPLIT  (T,X;T\,T2):  Construct  from  tree  T  two  balanced  binary  trees  T\  and  r2,  such  that 
A(7i)  is  the  portion  of  A(T)  from  its  leftmost  leaf  to  X,  and  A(T2)  is  the  remaining  portion  of 
A(T).  Tree  T  is  destroyed  by  the  operation. 

SPLICE  (Tx,T2;T):  Construct  from  the  balanced  binary  trees  Tx  and  T2  a  new  balanced 
binary  tree  T  such  that  A(T)  is  the  concatenation  of  A (Tj)  and  A (T2),  with  A (Tx)  occurring  to  the 
left  of  A(T2).  Trees  Tx  and  T2  are  destroyed  by  the  operation. 

Let  m  be  the  number  of  leaves  of  tree  T.  Standard  techniques  allow  the  performance  of 
each  of  the  above  operations  in  O  (log m)  time  [55,  pp.  52-53]. 

Regarding  the  update  operations  on  the  planar  sr- graph  G,  the  syntactic  transformations  on 
the  left-  and  right-sequence  of  G  correspond  to  performing  0(1)  insertions/deletions  and 
SPUTISPUCE  operations  on  the  trees  Ti  and  Tr.  Notice  that  the  elements  of  V^jEkjF 
involved  in  the  update  identify  the  elements  of  the  left- sequence  that  are  inserted,  deleted,  or  are 
at  the  boundary  of  substrings  to  be  swapped.  For  example,  the  algorithm  for  operation  INSERT- 
EDGE  is  as  follows: 

Algorithm  INSERTEDGE  (e,  u,  v,f  ;/j  ,/2 ) 

(1)  Determine  the  relative  order  of  u,  v,  and  /  in  the  left-sequence  of  G  by  applying  the  order- 
query  algorithm  outlined  in  the  proof  of  Theorem  6.6.  This  determines  which  of  the  four 
cases  of  Theorem  6.4  applies. 
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(2)  Access  leaves  u,  v,  and  /  in  tree  Ti  and  remove  them.  Also,  by  means  of  at  most  three 
SPLIT  operations,  construct  from  Ti  four  trees  associated  with  substrings  A,  B,  C,  and  D. 

(3)  Destroy  leaf  / and  create  new  leaves  f\  and  fi- 

(4)  Assemble  the  updated  tree  Ti  from  the  leaves  u,  v,  f\ ,  and  fi,  and  from  the  trees  associated 
with  A,  B ,  C,  and  D  by  a  sequence  of  SPLICE  operations  and  insertions.  The  correct  left- 
to-right  order  of  these  constituents  is  selected  according  to  the  specifications  of  Theorem 
6.4. 

(5)  Perform  the  above  Steps  1-4  on  the  right-sequence  and  tree  TR. 

Analogous  algorithms  can  be  formulated  for  the  other  update  operations,  and  we  have 

Theorem  6.7  The  restructuring  of  trees  TL  and  T#  after  any  one  of  the  update  operations 
INSERTEDGE,  REMOVEEDGE,  EXPANDVERTEX,  and  CONTRACTVERTEX  can  be  per¬ 
formed  in  0(log  ri)  time. 

6.4.  Applications 

The  general  framework  for  the  maintenance  of  orders  <l  and  <r  in  a  planar  sr-graph  can 
be  profitably  used  in  three  interesting  applications:  (i)  dynamic  point  location  in  monotone  sub¬ 
divisions,  («)  dynamic  transitive-closure  query  in  planar  st- graphs,  and  (iii)  dynamic  contact- 
chain  query  in  convex  subdivisions. 

6.4.1.  Dynamic  planar  point  location  revisited  A  monotone  subdivision  31  is  associated  with 
a  planar  rr-graph  G  such  that  (see  Fig.  6.8): 

(1)  the  vertices  of  G  are  the  vertices  of  91,  plus  two  special  vertices  s  and  t,  associated  with  ver¬ 
tices  at  infinity  in  the  vertical  direction; 

(2)  the  arcs  of  G  are  associated  with  the  edges  of  91,  and  oriented  from  the  lower  to  the  upper 
endpoint  (horizontal  edges  are  oriented  from  left  to  right);  also  G  contains  arcs  connecting 
consecutive  vertices  of  91  at  infinity. 

Note  that  the  vertices  on  the  external  boundary  of  G  are  the  vertices  of  91  at  infinitity,  plus  s  and 
r.  Since  monotone  chains  in  91  correspond  to  directed  paths  in  G,  Theorem  3.2  of  Chapter  3  can 
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(a)  (b) 

Figure  6.8  (a)  Monotone  subdivision;  (b)  The  planar  si- 
graph  associated  with  the  monotone  subdivision  of  part  (a). 

be  viewed  as  the  geometric  counterpart  of  Lemma  6.9. 

Let  ri  'r 2>  '  ‘  *  >rf  be  the  regions  of  'K,  sorted  according  to  some  total  order  <  compatible 
with  relation  — i.e.,  rt— implies  /  <  j.  The  common  boundary  of  the  regions  with  index  less 
than  or  equal  to  i  and  of  the  ones  with  index  greater  than  i  is  a  separator  of  9?  denoted  or,. 
Clearly,  G;  is  to  the  left  of  a;  for  i  <  j.  Hence,  the  order  <  defines  a  complete  family  of  separa¬ 
tors  £={cJi,  •  •  •  Notice  that  region  is  the  portion  of  the  plane  between  separators 

and  Gj.  Conversely,  a  complete  family  Z=  {Gj,  •  •  •  .<5f  \ }  of  /  - 1  separators  defines  a 
total  order  rltr2,  •  •  •  ,rf  on  the  regions  compatible  with  relation  where  r,  is  the  unique 
region  contained  between  separators  g,  _i  and  g,-. 

In  the  point  location  technique  of  Lee-Preparata  [35]  the  leaves  of  the  separator-tree  3  are 
the  regions  of  91,  sorted  from  left  to  right  according  to  < ,  and  the  sequence  of  the  nodes  of  3  in 
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symmetric  order  is  •  •  •  <Sf-\Tf.  A  REMOVECHAIN  operation  might  cause  the  order  of 

the  regions  given  by  3  to  become  inconsistent  with  relation  (see  Fig.  6.9),  and  the  consequent 
rearrangement  of  this  order  might  require  the  reconstruction  of  large  portions  of  3.  Our 
dynamic  point  location  technique  overcomes  the  above  obstacle  using  the  leftist  order  <i  on  the 
regions  of  SR  to  define  the  family  of  separators  I,  so  that  the  rearrangement  of  the  separator  tree 
after  an  update  can  be  performed  by  a  simple  transformation. 

With  regard  to  channels,  we  observe  that  the  channel  between  two  vertically  consecutive 
regions  rx  and  r*i  consists  exactly  of  the  vertices  and  edges  between  HIGH  (r j )  and  LOW (ri)  in 
the  left-sequence  of  G.  For  example,  the  left-sequence  of  the  subdivision  of  Fig.  6.10  (omitting 


Figure  6.9  Example  of  REMOVECHAIN  operation  which 
causes  the  order  of  the  regions  given  by  the  separator-tree 
to  become  inconsistent  with  relation  — ». 
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vertices  and  edges  at  infinity)  is  given  by  the  string  rx  ex  r-i  e 2  A  e$  r$  r4,  where  the  sub¬ 
string  A  =v\  «3  V2  e4  V3  is  associated  with  the  channel  between  rx  and  rj.  With  regard  to  the 
INSERTCHAIN  operation,  we  observe  that  Theorem  3.8  of  Chapter  3  is  the  geometric  counter¬ 
part  of  Theorem  6.4. 

6.4.2.  Transitive-closure  query  Recall  that  a  transitive-closure  query  on  a  planar  sr-graph  G 
consists  of  determining  the  existence  of  a  directed  path  between  vertices  u  and  v  of  G.  Such 
query  is  equivalent  to  testing  whether  both  u<^v  and  u  <r  v  so  that,  by  Theorem  6.6,  it  takes 
O  (log  n)  time.  This  establishes  Theorem  6.1  of  Section  6.1. 

A  variant  of  query  reports  a  path  between  u  and  v,  and  can  be  executed  in  time 
0(logn+k),  where  k  is  the  number  of  path  edges.  First,  we  query  (in  O(logn)  time)  the 
existence  of  a  path  between  u  and  v.  Suppose  that  such  path  exists  and,  say,  uTv.  We  know  that 
the  leftmost  path  from  u  to  t  and  the  leftmost  path  from  s  to  v  have  at  least  one  vertex  in  com¬ 
mon.  Resorting  to  a  standard  DCEL  representation  of  the  planar  sr-graph  (see  [47,  pp.  15-17]), 
we  can  trace  each  of  these  two  paths.  Alternating  between  them  one  edge  at  a  time,  we  trace  the 


Figure  6.10  Example  of  channel. 
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path  between  u  and  t  forward  from  u,  and  the  path  between  s  and  v  backward  from  v.  In  this 
traversal  we  mark  each  visited  vertex.  The  process  terminates  when  we  reach  a  vertex  for  the 
second  time.  If  k  is  the  length  of  the  path  to  be  reported,  clearly  at  most  2 k  vertices  have  been 
visited  by  the  process.  This  establishes  that  the  report-type  query  is  executed  in  time 
0(logrt+£). 

6.4.3.  Contact-chain  query  We  can  reformulate  the  problem  of  contact  chains  by  assuming 
without  any  sacrifice  of  generality  that  the  reference  direction  6  is  always  the  x-axis.  In  this  set¬ 
ting,  we  have  that  region  rx  pushes  region  r 2  if  and  only  if  rj  is  to  the  left  of  r2.  Hence,  the 
transitive  closure  of  the  “push”  relation  is  the  same  as  relation  and  variations  of  0 
correspond  to  rotations  of  the  subdivision. 

We  assume,  with  negligible  loss  of  generality,  that  the  slopes  of  the  edges  are  all  distinct 
(In  the  case  of  parallel  edges,  a  virtual  perturbation  of  their  slopes  achieves  this  simplifying  con¬ 
dition.)  Thus,  if  we  continuously  rotate  the  subdivision,  only  one  edge  at  a  time  becomes  hor¬ 
izontal.  An  elementary  clockwise  rotation  from  a  given  position  of  91  is  the  minimal  nonzero 
clockwise  rotation  such  that  an  edge  becomes  horizontal.  An  elementary  counterclockwise  rota¬ 
tion  is  correspondingly  defined.  Thus,  a  full  2n-rotation  of  91  is  a  sequence  of  elementary  rota¬ 
tions. 

Since  a  convex  subdivision  91  is  also  a  monotone  subdivision,  we  consider  the  planar  st- 
graph  G  associated  with  91,  and  its  dual  G*.  It  is  easy  to  see  that  contact  chains  of  91  are  in 
one-to-one  correspondence  with  paths  in  the  graph  G*. 

We  consider  the  following  update  operations  on  91: 

INSERTPOINT (v,e;ei,e2)-  Split  the  edge  e=(u,w)  into  two  edges  e\  =(n,v)  and 
e2  =  (v,w),  by  inserting  vertex  v. 

REMOVEPOINT  (v,e1,e2;e):  Let  v  be  a  vertex  of  degree  2  whose  incident  edges, 
e\  =(u,v)  and  e2  =(v,w),  are  on  the  same  straight  line.  Remove  v  and  replace  e\  and  e2  with 
edge  e=(u,w). 

INSERTSEGMENT (e,u,v,r\ri,r2):  Add  edge  e=(u,v)  inside  region  r,  which  is  decom¬ 
posed  into  regions  rt  and  r2,  with  r  j  to  the  left  of  e  and  r2  to  the  right 


REMOVESEGMENT  (e,u,v,r i, r^r):  Remove  edge  e  =  (n,v)  and  merge  the  regions  rx 
and  r-i  formerly  on  the  two  sides  of  e  into  region  r .  [  The  operation  is  allowed  only  if  the  subdi¬ 
vision  9T  so  obtained  is  convex.  ] 

ROTATE  (8):  Perform  an  elementary  rotation  of  the  subdivision  SR.  The  binary  parameter 
5  indicates  whether  the  rotation  is  clockwise  or  counterclockwise. 

To  maintain  information  on  the  paths  of  G*,  we  use  the  theoretical  framework  developed  in 
Sections  6.2  and  6.3,  and  exchange  the  roles  of  G  and  G*.  Operations  INSERTPOINT  and 
REMOVEPOINT  on  SR  correspond  to  performing  operations  INSERTEDGE  and  REMOVEEDGE 
on  G  *.  Operations  INSERTSEGMENT  and  REMOVESEGMENT  on  91  correspond  to  performing 
operations  EXPAND  VERTEX  and  CONTRACTVERTEX  on  G*.  These  correspondences  yield 
algorithms  for  performing  contact-chain  queries  and  insertions/deletions  of  vertices  and  edges  in 
time  O  (log  /»). 

With  regard  to  the  operation  ROTATE,  let  e  be  the  edge  of  91  that  becomes  horizontal  at 
some  time  during  the  rotation.  The  effect  of  such  rotation  on  G*  is  to  invert  the  direction  of  the 
dual  edge  e*  of  e  (see  Fig.  6.11).  Hence,  operation  ROTATE  on  91  corresponds  to  performing  a 
REMOVEEDGE  operation  on  G*,  followed  by  an  INSERTEDGE  operation  of  the  same  edge  in 
the  reverse  orientation. 

Let  the  azimuth  of  a  directed  edge  be  defined  counterclockwise  with  respect  to  the  x-axis, 
so  that  it  lies  in  the  range  [0,7t],  The  edge  e  involved  in  the  rotation  can  be  identified  by  main¬ 
taining  a  list  of  the  edges  of  91  sorted  by  increasing  azimuth.  Specifically,  the  edge  involved  in  a 
clockwise  (counterclockwise)  elementary  rotation  is  the  first  Oast)  edge  of  this  list,  and  is  moved 
to  the  end  (front)  of  the  list  after  the  rotation.  The  list  is  implemented  as  a  balanced  binary  tree, 
so  that  edges  can  be  efficiently  inserted/deleted  as  specified  by  the  operations  INSERTPOINT, 
REMOVEPOINT,  INSERTSEGMENT,  and  REMOVESEGMENT. 

In  conclusion,  all  the  update  operations  have  G(Iogn)  time  complexity,  which  establishes 


Theorem  6.2  of  Section  6.1. 


Figure  6.11  (a)  Convex  subdivision  91.  (b)  Subdivision  91 
after  an  elementary  clockwise  rotation  (edge  e  becomes 
horizontal  at  some  time  during  the  rotation),  (c)  Graph  G* 
before  the  rotation,  (d)  Graph  G*  after  the  rotation  (the 
orientation  of  edge  e*  is  reversed). 


CHAPTER  7 


DYNAMIC  PLANAR  GRAPH  EMBEDDING 


7.1.  Introduction 

Embedding  a  graph  in  the  plane  is  a  fundamental  problem  in  several  areas  of  computer  sci¬ 
ence,  including  circuit  layout,  graphics,  and  computer-aided  design.  The  problem  of  testing  the 
planarity  and  of  constructing  a  planar  embedding  of  a  graph  has  been  extensively  studied  in  the 
past  years,  and  the  development  of  linear  time  algorithms  for  it  has  brought  significant  advances 
in  algorithm  design  and  analysis  [7,24].  Nevertheless,  as  confirmed  by  recent  results  [9, 16, 17], 
graph  planarity  is  still  a  vital  area  of  research,  rich  in  interesting  issues  to  be  explored. 

In  this  chapter  we  consider  the  problem  of  incrementally  constructing  a  planar  embedding 
of  a  graph.  We  investigate  a  dynamic  data  structure  that  allows  us  to  perform  efficiently  the  fol¬ 
lowing  operations: 

(1)  queries:  given  two  vertices  u  and  v,  determine  whether  there  is  a  face  of  the  current 
embedding  whose  boundary  contains  both  u  and  v; 

(2)  updates:  modify  on-line  the  current  embedding  by  adding  and/or  removing  vertices  and 
edges. 

The  performance  of  such  a  data  structure  will  be  measured  in  terms  of  (1)  the  space  requirement, 
(2)  the  query  and  update  times,  and  (3)  the  preprocessing  time. 

Formally,  our  problem  can  be  defined  as  follows:  Let  G  be  a  planar  graph  embedded  in  the 
plane,  referred  to  henceforth  as  a  plane  graph .  For  generality,  we  allow  G  to  have  multiple 
edges,  and  we  denote  with  n  and  m  the  number  of  vertices  and  edges  of  G,  respectively.  We 
consider  the  dynamic  embedding  problem ,  which  consists  of  performing  the  following  opera¬ 
tions  on  G: 

TEST  (u,v):  Test  whether  there  is  a  face  /  that  has  both  vertices  u  and  v  on  its  boundary. 
If  such  a  face  exists,  output  its  name. 
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UST  (u,v):  List  all  the  faces  that  have  both  vertices  u  and  v  on  their  boundary. 

INSERTEDGE  {e,u,v,f-jx  ,fiY  Add  the  edge  e=(u,v)  to  G  inside  face  /,  which  is  decom¬ 
posed  into  faces  f\  and  f^.  Vertices  u  and  v  must  both  be  on  the  boundary  of  face  /. 

INSERTVERTEX  {e,v\e\,e-iY  Split  the  edge  e=(u,w)  into  two  edges  e\  =  (u,v)  and 
=  (v,w),  by  adding  vertex  v. 

REMOVEEDGE  (e,u,v,fi  Remove  the  edge  e=(u,v),  and  merge  faces  /i  and  fi 

formerly  on  the  two  sides  of  e  into  face  /. 

REMOVEVERTEX  (e\,e2,v,e):  Let  v  be  a  vertex  of  degree  two.  Remove  v  and  replace 
its  incident  edges  e\  =(u,v)  and  =(v,w)  with  edge  e  -  ( u,w ). 

It  is  a  relatively  simple  exercise  to  show  that  the  above  repertory  of  operations  is  complete  for 
the  class  of  sf-2-connectible  planar  graphs  defined  in  Section  7.2. 

The  dynamic  embedding  problem  naturally  arises  in  interactive  CAD  layout  environments. 
Applications  include  the  design  of  integrated  circuits,  modon  planning  in  robodcs,  architectural 
floor  planning,  and  graphic  editing  of  block  diagrams. 

We  present  a  data  structure  that  uses  0(m)  space,  supports  all  of  the  above  operations  in 
O  (logm)  time,  and  can  be  constructed  in  O  ( m )  time.  If  G  is  simple,  i.e.  it  has  no  multiple  edges 
and  no  self-loops,  then  m=0(n),  and  the  above  bounds  become  O(n)  space  and  preprocessing 
time,  and  O(logn)  query  and  update  times.  In  addition  to  the  good  space/time  performance 
from  a  theoretical  viewpoint,  our  data  structure  is  also  practical  and  easy  to  implement,  and 
therefore  suited  for  real-world  applications. 

These  results  are  obtained  by  maintaining  on-line  an  orientation  of  the  graph,  called  spheri¬ 
cal  st-orientation  and  exploiting  the  partial  order  among  the  vertices,  edges,  and  faces  induced 
by  this  orientation.  Besides  its  relevance  to  this  problem,  the  concept  of  spherical  sr-orientation 
is  of  theoretical  interest  in  its  own  right,  and  extends  the  results  on  bipolar  orientations  and 
cylindric  orientations  of  planar  graphs  presented  in  [49, 51,52]. 

This  work  constitutes  also  a  first  step  toward  the  development  of  an  efficient  data  structure 
for  the  dynamic  planarity  testing  problem,  which  consists  of  performing  the  following  operations 
on  a  planar  graph  G:  (1)  testing  if  a  new  edge  can  be  added  to  G  so  that  the  resulting  graph  is 
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itself  planar;  (2)  adding  and  removing  vertices  and  edges. 

The  rest  of  this  chapter  is  organized  as  follows:  Section  7.2  contains  definitions  and  prel¬ 
iminary  results.  Section  7.3  deals  with  orientations  of  planar  graphs.  In  Section  7.4,  we  study 
the  topological  location  problem,  which  consists  of  performing  operations  TEST  and  LIST.  Sec¬ 
tion  7.5  describes  the  full-fledged  data  structure  for  the  dynamic  embedding  problem.  Finally, 
further  applications  are  discussed  in  Section  7.6. 

7 2.  Preliminaries 

We  consider  only  finite  connected  graphs  without  self-loops.  We  allow  multiple  edges 
between  two  vertices.  For  the  basic  terminology  about  graphs  and  planarity,  see  [6, 15].  Unless 
otherwise  specified,  paths  and  cycles  of  a  directed  graph  are  assumed  to  be  directed. 

First,  we  recall  some  definitions  on  graph  connectivity.  A  cutvertex  of  a  graph  G  is  a  vertex 
whose  removal  disconnects  G.  Graph  G  is  said  to  be  2 -connected  if  it  has  no  cutvertices,  and  1- 
connected  otherwise.  A  block  of  a  1 -connected  graph  G  is  a  maximal  2-connected  subgraph  of 
G.  The  proper  vertices  of  a  block  B  of  G  are  the  vertices  of  B  that  are  not  cutvertices.  The 
block-cutvertex  tree  of  G  is  a  tree  whose  nodes  represent  the  blocks  and  cutvertices  of  G,  and 
whose  edges  connect  each  cutvertex  v  to  the  blocks  that  contain  v  (see  Fig.  7.1).  A  graph  G  is 
st-2 -connectible  if  adding  the  edge  (s,t)  to  G  makes  G  2-connected  [36],  Clearly,  a  2-connected 
graph  is  also  sf-2-connectible  for  every  pair  of  vertices  s  and  t.  An  rt-numbering  of  a  graph  G 
with  vertex  set  V  is  a  bijection  V— >{  1,2,  •  •  • ,  I V  I }  such  that  every  vertex  v*s,f  has  neighbors 
u  and  v  with  £(w)  <  2;(v)  <  £(w).  A  graph  admits  an  st-numbering  if  and  only  if  it  is  st-2- 
connectible  [36]. 

Regarding  the  dynamic  embedding  problem,  we  assume  that  the  vertices,  edges,  and  faces 
of  the  graph  are  identified  by  names,  which  are  elements  of  an  ordered  set.  For  example,  names 
can  be  integers,  alphanumeric  strings,  or  pairs  of  coordinates.  The  total  order  among  names  will 
be  referred  to  as  alphabetic  order.  Regarding  the  complexity  analysis,  we  assume  that  a  name 
uses  0(1)  space,  and  that  the  alphabetic  comparison  between  two  names  can  be  done  in  0(1) 
time.  Also,  for  generality,  we  assume  that  the  query  and  update  operations  use  the  names  as 
input  parameters. 
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r  (a)  (b) 

Figure  7.1  (a)  A  1 -connected  graph,  and  (b)  its  block-cutvertex  tree. 

Let  m  be  the  current  number  of  edges  of  the  graph.  We  will  denote  the  storage  by 
Space  (m)  and  the  preprocessing  time  by  Preprocess  (m).  Also,  the  time  complexity  of  the  vari- 
0  ous  operations  will  be  denoted  by  Test  (m),  List  ( m,k ),  InsertEdge  (m),  InsertVertex  (m), 

Remove  Edge  (m),  and  RemoveVertex  (m),  where  k  is  the  number  of  faces  retrieved  by  the  LIST 
operation.  Finally,  throughout  this  chapter,  log*  means  max  { 1 ,  log2* }. 

_  73.  Orientations  of  Planar  Graphs 

A  spherical  st-graph  is  a  plane  digraph  G  such  that: 

Property  1:  G  has  exactly  one  source  (vertex  without  incoming  edges),  s,  and  exactly  one  sink 
(vertex  without  outgoing  edges),  t. 

Property  2:  Every  vertex  v  of  G  is  on  some  directed  simple  path  from  s  to  t. 

Property  3:  Every  directed  cycle  y  separates  s  from  t,  i.e.,  one  of  s  and  t  is  inside  the  region  of 
l  the  plane  bounded  by  y,  and  the  other  is  outside. 
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Wc  can  visualize  a  spherical  sr- graph  as  embedded  on  the  surface  of  a  sphere,  with  s  and  t 
at  the  South  and  North  poles,  respectively  (see  Fig.  7.2). 

The  concept  of  spherical  sr-graph  extends  the  one  of  planar  st-graph  introduced  in  [36], 
which  has  important  applications  in  the  test  of  graph  planarity  [36]  and  the  construction  of 
planar  drawings  [10,49, 51]  (see  Chapter  6  for  the  definition  and  properties  of  planar  sr-graphs). 
A  spherical  sr-graph  differs  from  a  planar  st-graph  because  it  admits  (directed)  cycles.  However, 
such  cycles  must  verify  the  aforementioned  Property  3. 

The  following  two  lemmas  show  that  spherical  sr-graphs  have  the  same  properties  as  planar 
sf-graphs  with  regard  to  the  circular  sequences  of  edges  that  are  incident  upon  a  vertex  (Lemma 
6.2  of  Chapter  6),  and  that  form  the  boundary  of  a  face  (Lemma  6.3  of  Chapter  6). 

Lemma  7.1  For  every  vertex  v  of  G,  the  incoming  (outgoing)  edges  appear  consecutively 
around  v  (See  Fig.  7.2(b)). 

Proof:  Assume,  for  a  contradiction,  that  there  is  a  vertex  v,  v*s,r,  for  which  the  lemma  is  not 
true  (see  Fig.  7.3).  Then  there  must  be  four  edges  incident  upon  v,  denoted  e1=(w1,v), 
e2  =  (v,w2),  «3  =(w3,v),  and  e4  =  (v,w4),  which  appear  in  this  order  counterclockwise  around  v. 
By  Property  2,  there  are  (directed)  paths  from  s  to  w  j  and  w3.  Let  s'  he  the  vertex  farthest  from 
s  that  is  on  both  these  paths.  We  denote  with  Kj  and  7t3  the  portions  of  such  paths  from  s'  to  W] 
and  w3,  respectively.  The  union  of  n j,  Jt3,  e\,  and  e3  forms  an  undirected  cycle  y,  which 
separates  w2  from  w4.  The  two  regions  of  the  plane  delimited  by  cycle  y  will  be  denoted  by  A 
and  B,  where  A  is  the  region  that  contains  vertex  w2.  We  assume  that  both  A  and  B  contain 
cycle  y.  By  Property  2,  there  must  be  paths  Jt2  and  rc4  from  w2  and  w4  to  r,  respectively.  Now, 
we  have  four  cases  for  the  relative  placement  of  s  and  r  with  respect  to  cycle  y.  If  both  s  and  t 
are  in  A ,  then  n4  intersects  y  at  some  vertex  (see  Fig.  7.3(a)).  This  creates  a  cycle  that  does  not 
separate  s  from  r.  If  s  is  in  A  and  t  is  in  B,  then  rc2  intersects  y,  and  again  we  have  a  cycle  that 
does  not  separate  s  from  t  (see  Figs.  7.3(b-c)).  The  cases  when  both  s  and  t  are  in  B,  or  s  is  in  B 
and  t  is  in  A,  are  treated  similarly.  We  conclude  the  proof  by  observing  that  in  all  cases  we  have 
a  contradiction  to  Property  3.  □ 
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Lemma  7.2  For  every  face  /  of  G,  the  boundary  of  /  consists  of  two  directed  paths  with  com¬ 
mon  origin  and  destination  (See  Fig.  7.2(c)). 

Proof:  Assume,  for  a  contradiction,  that  there  is  a  face  /  for  which  the  lemma  is  not  true  (see 
Fig.  7.4).  Then  there  are  distinct  vertices  u  and  v  on  the  boundary  of  /  such  that  the  edges  of  the 
boundary  of  /incident  upon  them  are  all  outgoing.  We  denote  these  edges  with  e\=(u,w\), 
e2-(u,W2),  e-$  —  (v, w 3 ),  and  e4  =  (v,w4),  in  counterclockwise  order  on  the  boundary  of /.  From 
Property  2,  there  are  directed  paths  from  s  to  u  and  v.  Let  s'  be  the  vertex  farthest  from  s  that  is 
on  both  these  paths.  We  denote  with  nu  and  jcv  the  portions  of  such  paths  from  s'  to  u  and  v, 
respectively.  Also,  we  denote  by  n'  the  portion  of  nu  from  s  to  s'.  The  union  of  7CU,  7tv,  and  the 
portion  of  the  boundary  of  /  clockwise  from  v  to  u  forms  an  undirected  cycle  y,  which  contains 
vertices  w»2  and  w3-  The  two  regions  of  the  plane  delimited  by  cycle  y  will  be  denoted  by  A  and 
B,  where  A  is  the  region  that  does  not  contain  face  /.  We  assume  that  both  A  and  B  contain  cycle 
y.  From  Property  2,  there  must  be  paths  Jtj,  Jt2,  1x3,  and  714  from  wj,  W2,  w3,  and  w4  to  t, 
respectively.  Now,  we  have  four  cases  for  the  relative  placement  of  s  and  t  with  respect  to  cycle 
Y 

(0  s  and  t  are  both  in  A. 

Path  Jtj  must  intersect  at  least  one  of  nu  and  Kv.  If  it  intersects  first  JtM,  then  we  have 
immediately  a  cycle  that  does  not  separate  s  from  t  (see  Fig.  7.4(a)).  Otherwise,  let  r  be  the 
intersection  vertex  of  7ti  with  7tv.  Path  jc4  must  intersect  either  nv,  or  the  portion  of  Tti 
from  wi  tor  and  then  ku.  In  both  cases,  we  have  again  a  cycle  that  does  not  separate  s 
from  t  (see  Fig.  7.4(b-c)). 

(»)  s  is  in  A  and  t  is  in  B. 

Let  tC  be  a  path  from  s  to  s'.  Path  Jt2  must  intersect  at  least  one  of  Jt„  and  Ky.  If  it  inter¬ 
sects  first  nu  at  vertex  r,  then  we  have  a  cycle  formed  by  edge  e2{u,w>2),  the  subpath  of  JI2 
from  W2  to  r,  and  the  subpath  of  ku  from  r  to  u.  If  this  cycle  does  not  separate  s  from  t,  we 
are  done  (see  Fig.  7.4(d)).  Otherwise,  path  712  must  intersect  7t'  at  some  vertex  q  and  path 
713  must  intersect  the  directed  path  consisting  of  the  portion  of  K2  from  W2  to  q,  the  portion 
of  n'  from  q  to  s',  and  path  nv.  Hence,  we  have  a  cycle  that  does  not  separate  s  from  t  (see 
Fig.  7.4(e)).  The  last  subcase  to  examine  is  when  K2  intersects  first  Tty.  Let  r  be  the  first 
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proof  of  Lemma  7.2  (a-c) 


intersection  of  JC2  with  nv.  Path  n 3  must  intersect  the  directed  path  consisting  of  the  portion 
of  Jt2  from  h>2  to  r  and  the  portion  of  path  jiv  from  r  to  v.  Hence,  713  forms  a  cycle  with  the 
above  path.  If  this  cycle  separates  s  from  t,  then  rtj  must  intersect  7t'  at  some  vertex  q  and 
form  a  cycle  that  does  not  separate  s  from  t  (see  Fig.  7.4(f)). 

(iii)  and  (iv)  s  and  t  are  both  in  B,  or  s  is  in  B  and  t  is  is  in  A. 

These  cases  are  analogous  to  the  ones  above,  and  their  treatment  is  omitted  for  brevity.  In 
all  cases  we  have  a  contradiction  to  Property  3,  and  the  proof  is  completed.  □ 

Motivated  by  the  previous  lemmas,  we  introduce  for  spherical  sr-graphs  the  same  terminol¬ 
ogy  as  for  planar  sr-graphs,  given  in  Section  6.2  of  Chapter  6,  which  is  repeated  here  for  the 
reader’s  convenience.  Let  V,  £,  and  F  denote  the  set  of  vertices,  edges,  and  faces  of  G,  respec¬ 
tively.  For  each  element  x  of  VuE  uF,  we  define  vertices  LOW{x)  and  HIGH  lx),  and  faces 
LEFTQc)  and  RIGHT  (x),  as  follows: 

(1)  If  x  =  v  e  V,  we  define  LOW(v )  =  HIGHiy ) = v.  Also,  with  reference  to  Lemma  7. 1  and  Fig. 
7.2(b),  we  denote  by  LEFT(v)  and  RIGHTly)  the  two  faces  that  separate  the  incoming  and 
outgoing  edges  of  a  vertex  v*y,r,  where  LEFTly)  is  the  face  to  the  left  of  the  leftmost 
incoming  and  outgoing  edges,  and  RIGHT (v)  is  the  face  to  the  right  of  the  rightmost 
incoming  and  outgoing  edges. 

(2)  If  x=ee  E,  we  define  LOW(e)  and  HIGH(e)  as  the  tail  and  head  vertices  of  e,  respec¬ 
tively.  Also,  we  denote  by  LEFT(e)  and  RIGHT (e)  the  faces  on  the  left  and  right  side  of  e, 
respectively. 

(3)  If  x  =/€  F,  we  denote  by  LOW(f)  and  HIGH(f)  the  two  vertices  that  are  the  common  ori¬ 
gin  and  destination  of  the  two  paths  forming  the  boundary  of  /  (see  Lemma  7.2  and  Fig. 
7.2(c)).  Vertices  LOW (J )  and  HIGH(f )  are  called  the  extreme  vertices  of  face  /.  Also,  we 
define  LEFT(f ) = RIGHT (f)  =  /.  Finally,  the  two  directed  paths  forming  the  boundary  of  / 
are  called  the  left  path  and  right  path  of/,  respectively. 

We  assume  that  the  left  and  right  paths  of  /  do  not  include  their  endpoints,  so  that  these 
paths  and  the  extreme  vertices  of  /  form  a  partition  of  the  boundary  of  /.  Notice  that  a  vertex  v  is 
on  the  left  (respectively,  right)  path  of  face  /  if  and  only  if  RIGHT(v)=f  (respectively. 


LEFT(v)  -f).  Hence,  from  LEFT(v),  RIGHT(v),  LOW(f),  and  HIGHif)  we  can  decide  in  O  (1) 
time  whether  vertex  v  is  on  the  boundary  of  face  /. 

Let  G  be  a  plane  graph,  and  /  and  /  two  distinct  vertices  of  G.  A  spherical  st-orientation  of 
G  is  a  spherical  //-graph  whose  undirected  version  is  isomorphic  to  G.  G  is  said  to  be  st- 
orierttable  if  it  admits  a  spherical  //-orientation.  The  following  theorem  provides  a  characteriza¬ 
tion  of  /z-orientable  graphs,  and  is  similar  to  the  characterization  of  //-numerable  graphs  given  in 
[36]. 

Theorem  7.1  Let  G  be  a  plane  graph.  The  following  statements  are  equivalent: 

(1)  G  is  //-orientable; 

(2)  G  admits  an  acyclic  spherical  //-orientation; 

(3)  G  admits  an  //-numbering; 

(4)  G  is  //-2-connectible. 

Also,  there  are  O  (m)  time  algorithms  for  testing  if  G  is  //-orientable  and  constructing  a  spherical 
//-orientation  for  G. 

Proof:  It  is  proved  in  [36]  that  (4)  =>  (3).  Given  an  //-numbering  for  G,  we  can  construct  an 
acyclic  spherical  .//-orientation  by  orienting  each  edge  from  the  lowest  to  the  highest  numbered 
vertex.  We  have  thus  (3)  =*(2).  Clearly,  (2)  =>  (1).  To  complete  the  proof  of  the  characteriza¬ 
tion,  we  show  that  (1)  =>  (4).  Assume,  for  a  contradiction,  that  G  is  not  //-2-connectible.  Then 
there  is  a  cutvertex  v  of  G  such  that  one  of  the  components  generated  by  the  removal  of  v, 
denoted  by  C,  does  not  contain  either  s  or  /.  Let  u  be  a  vertex  of  C.  Any  path  from  s  tot  through 
v  is  not  simple,  which  is  a  contradiction. 

The  algorithm  for  testing  if  G  is  //-orientable  consists  of  verifying  that  each  cutvertex  of  G 
belongs  to  exactly  two  blocks  (connected  components)  of  G  and  that  each  block  of  G  contains  no 
more  than  two  cutvertices.  This  takes  O  (m)  time.  Finally,  since  computing  the  //-numbering  of 
a  planar  graph  can  be  done  in  0  (m)  time  [14],  we  also  have  the  result  that  constructing  a  spheri¬ 
cal  //-orientation  takes  O  ( m )  time.  □ 
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Now,  wc  turn  our  attention  to  operations  that  update  a  spherical  sr- graph  by  additions  and 
deletions  of  vertices  and  edges.  The  operations  INSERTEDGE,  INSERTVERTEX,  REMO- 
VEEDGE,  and  REMOVEVERTEX,  defined  in  the  introduction,  are  suitable  for  this  purpose. 
However,  further  restrictions  must  be  imposed  on  their  applicability  in  order  to  ensure  that  the 
resulting  graph  is  itself  a  spherical  sr- graph. 

Lemma  7.3  Let  G  be  a  spherical  rr-graph,  and  G'  be  the  directed  plane  graph  obtained  by  per¬ 
forming  operation  IT  on  G.  Depending  on  n,  G'  is  a  spherical  rr- graph  if  and  only  if 

(1)  for  H=INSERTEDGE  ie,u,vff\ff),  edge  e  must  not  create  a  (directed)  cycle  with  the 
edges  of  face/, 

(2)  for  n = INSERTVERTEX  (e,v;e  i,«2)»  there  is  no  restriction; 

(3)  for  n=REMOVEEDGE e  =  (u,v)  must  be  an  edge  such  that  deg+(«)£2 
and  deg~(v)22  (where,  as  usual,  deg+(w)  and  deg"(w)  denote  the  outdegree  and  indegree 
of  a  vertex  w,  respectively); 

(4)  for  II = REMOVEVERTEX (e  l.ej.v.e),  v  must  be  a  (degree-2)  vertex  distinct  from  s  and  r. 

Proof:  The  proof  of  (2),  (3),  and  (4)  is  straightforward.  For  operation  INSERTEDGE,  we  con¬ 
sider  two  cases.  First,  assume  that  there  is  a  path  7t  on  the  boundary  of  /  from  u  to  v.  If  edge 
e=(u,v)  creates  a  cycle  that  does  not  separate  s  from  t,  then  by  replacing  e  with  tc  we  have  that 
G  already  had  a  nonseparating  cycle,  a  contradiction.  Now,  if  there  is  no  path  on  the  boundary 
of  /  from  u  to  v,  we  are  in  the  situation  shown  in  Fig.  7.5.  Let  y  be  the  cycle  that  does  not 
separate  s  from  t,  and  iti  and  *2  be  paths  from  s  to  LOWif)  and  from  HIGHif)  to  t,  respec¬ 
tively.  One  of  these  two  paths,  say  Jti ,  must  intersect  y  at  some  vertex  x.  This  implies  that  G 
already  had  a  nonseparating  cycle,  formed  by  the  subpath  of  y  from  v  to  x,  the  subpath  of  Jti 
from  x  to  LOW if),  and  the  subpath  of  the  right  path  of  /from  LOW if)  to  v.  Again,  we  reach  a 
contradiction.  □ 

Corollary  7.1  Let  G  be  a  spherical  rr-graph,  and  G  j  and  G  2  be  the  graphs  obtained  by  perform¬ 
ing  operations  INSERTEDGE  ie,u,v,f  ,f\\fi)  and  INSERTEDGE  ie,v,u,f  ,f  iff)  on  G,  respec¬ 
tively,  where  both  vertices  u  and  v  are  on  the  boundary  of  face  /.  Then  at  least  one  of  G !  and 
G2  is  a  spherical  rr- graph. 
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Figure  7.5  Example  for  the  proof  of  Lemma  7.3 


Proof:  By  Lemma  7.3,  we  only  have  to  ensure  that  the  edge  to  be  added  does  not  form  directed 
cycles  with  the  boundary  off.  Hence,  if  there  is  a  (directed)  path  in  the  boundary  of  /from  u  to 
v,  we  perform  operation  INSERTEDGE  (e,u,v,/,/i \fi),  while  if  there  is  a  (directed)  path  from  v 
to  u,  we  perform  operation  INSERTEDGE  ;/2).  Both  operations  are  allowed  when  no 

directed  path  exists  between  u  and  v.  □ 


7.4.  Topological  Location 

In  this  section  we  consider  the  topological  location  problem ,  which  consists  of  performing 
efficiently  the  TEST  and  LIST  operation  on  a  plane  graph. 
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7.4.1.  rf*orientable  graphs  Let  G  be  a  spherical  st-graph.  From  Lemma  7.2,  vertices  u  and  v 
of  G  are  on  the  boundary  of  face  /if  and  only  if  one  of  the  following  cases  occurs  (see  Fig.  7.6): 

Case  1:  Both  u  and  v  are  nonextreme  vertices  of  /,  i.e., 

(f=LEFT(u)  or  f= RIGHT iu))  and  (f=LEFT(v )  or  f= RIGHT (v)); 

Case  2:  u  is  an  extreme  vertex  of  /  and  v  is  not,  i.e., 

C u  =LOW(f)  or  u  = HIGH  if))  and  (f=LEFT(v)  or  f= RIGHT  (v)); 

Case  3:  v  is  an  extreme  vertex  of  /  and  u  is  not,  i.e., 

(v  =  LOW  if)  or  v  = HIGH  if))  and  if = LEFT  {u)  or  f= RIGHT  iu)); 

Case  4:  both  u  and  v  are  extreme  vertices  of  /,  i.e., 

iu  =LOW(f)  and  v=HIGHif))  or  iu  =  HIGH  if)  and  v  -LOWif)). 

In  order  to  check  Cases  1-3,  we  store  for  each  vertex  v  the  faces  L£FT(v)  and  RIGHTiv), 
and  for  each  face  /  the  vertices  HIGH  if)  and  LOW  if).  For  each  of  these  cases,  the  test  is  car¬ 
ried  out  in  O  (1)  time.  To  check  the  remaining  Case  4,  we  store  with  each  vertex  v  a  search  table 
BELOWiy)  whose  elements  are  the  pairs  if,LQW(J))  such  that  v -HIGHif),  sorted  according 
to  the  alphabetic  order  of  the  name  of  LOWif).  Hence,  the  test  for  Case  4  consists  of  searching 
for  vertex  u  in  BELOW (v)  and  for  vertex  v  in  BELOW iu),  which  takes  G(logm)  time.  This 


Figure  7.6  The  four  cases  for  two  vertices  on  the  same  face 


proves  the  following  theorem: 


Theorem  7.2  There  exists  a  data  structure  for  the  topological  location  problem  in  spherical  st- 

graphs  with  the  following  performance: 

Space  (m)  =  Preprocess  ( m)  =  0 (m); 

Test  (m)  =  0  (log  m );  List  ( m,k )  =  0  (log  m  +  k). 


Corollary  7.2  There  exists  a  data  structure  for  the  topological  location  problem  in  w-orientable 

graphs  with  the  following  performance: 

Space  (m)= Preprocess  (m)=0(m)\ 

Test  (m)  =  0(logm);  List  (m,k)  =  0  (log  m  +k). 

Proof:  Construct  a  spherical  sr-orientation  for  G  and  apply  Theorem  7.2.  □ 


7.4.2.  General  plane  graphs  For  plane  graphs  that  are  not  sr-orientable,  the  algorithm  of  the 
preceding  subsection  cannot  be  directly  applied.  Instead,  we  will  combine  the  above  technique 
with  a  data  structure  that  takes  into  account  the  plane  arrangement  of  the  blocks  of  the  graph. 

In  the  following,  we  will  be  interested  in  preprocessing  a  graph  G  in  order  to  determine 
quickly  whether  there  is  a  block  that  contains  two  given  vertices  u  and  v.  This  can  be  done 
efficiently  by  orienting  the  block-cutvertex  tree  T  of  G  so  that  it  becomes  a  rooted  source  tree 
with  an  arbitrarily  selected  block  at  the  root.  In  the  example  of  Fig.  7.1,  the  edges  of  the  block- 
cutvertex  tree  are  oriented  from  bottom  to  top.  Also,  we  store  in  every  vertex  v  a  pointer 
node(v),  where,  if  v  is  a  cutvertex,  node(v )  points  to  the  representative  node  of  cutvertex  v  in  7, 
while  if  v  is  a  proper  vertex  of  a  block  B ,  node(v)  points  to  the  representative  node  of  block  B  in 
7.  It  is  simple  to  verify  that  there  is  a  block  containing  vertices  u  and  v  if  and  only  if  one  of  the 
following  cases  is  verified: 

(1)  node («)  =  node (v); 

(2)  node(u)  is  the  father  of  node(v)  in  7; 

(3)  node (v )  is  the  father  of  node(u)  in  7 ; 

(4)  the  fathers  of  node(u)  and  node(v)  are  the  same  node  of  7,  associated  with  a  block. 


We  conclude 


Lemma  7.4  There  exists  a  data  structure  with  O  (m)  space  and  preprocessing  time  that  allows 
for  testing  if  two  vertices  belong  to  the  same  block  in  O  (1)  time. 

Now,  let  G*  be  the  dual  graph  of  a  plane  graph  G.  A  face  of  G  that  is  a  cutvertex  of  G*  is 
called  a  outface  (of  G).  The  following  simple  lemma  shows  that  there  is  a  bijection  between  the 
blocks  of  G  and  the  ones  of  G  * . 

Lemma  7.5  [22,  p.  124,  ex.  1 1.4]  A  subset  of  edges  of  G  forms  a  block  of  G  if  and  only  if  their 
duals  form  a  block  of  G  * . 

We  define  now  a  partial  order  on  the  set  of  blocks,  vertices,  and  faces  of  G.  In  the  follow¬ 
ing  definitions  we  adopt  the  convention  that  the  faces  of  a  block  B  of  G  have  the  same  names  as 
the  corresponding  faces  in  G. 

Let  B  be  a  block  of  G.  The  outer  face  of  B,  denoted  OUTER  (B),  is  the  face  /  whose  boun¬ 
dary  includes  the  external  boundary  of  B.  In  turn,  B  is  called  an  inner  block  of  face  OUTER  (B). 
Notice  that  OUTER  (B)  is  either  a  cutface  or  the  external  face.  Now,  let  /  be  a  face  distinct  from 
the  external  face.  We  denote  with  OUTER(f)  the  block  B  such  that  /  is  an  internal  face  of  B. 
Finally,  let  v  be  a  vertex.  If  there  is  a  block  B  containing  v  such  that  v  is  not  on  the  external 
boundary  of  B,  then  we  define  OUTER(v)=B.  Otherwise,  there  is  a  (unique)  face /containing  v 
such  that,  for  no  block  B  containing  v,  B  =  OUTER  (/),  and  we  define  OUTER  (v)=/.  As  a 
straightforward  consequence  of  the  above  definitions,  we  have 

Lemma  7.6  The  graph  of  relation  OUTER  is  a  directed  source  tree,  whose  root  is  the  external 
face. 


Figure  7.7  shows  a  1-connected  graph  and  the  corresponding  OUTER  relation  for  the 
blocks,  cutvertices,  and  cutfaces. 


Theorem  7.3  Let  u  and  v  be  vertices  of  G  that  are  on  the  boundary  of  the  same  face  /.  If  u  and  v 
belong  to  the  same  block  B,  then  they  are  also  on  the  boundary  of  face  /  in  B.  Otherwise,  one  of 
the  following  cases  arises: 

(1)  /  =  OUTER (u)  and /=  OUTER (v); 

(2)  /= OUTER  (u)  and  v  is  on  face  /of  block  OUTER  if ); 

(3)  f= OUTER (v)  and  u  is  on  face/of  block  OUTERif). 

Proof:  A  straightforward  consequence  of  the  definition  of  the  OUTER  relation.  □ 


The  data  structure  for  the  topological  location  problem  in  general  plane  graphs  consists  of: 

(1)  A  data  structure  to  test  whether  two  vertices  belong  to  the  same  block,  see  Lemma  7.4. 

(2)  A  separate  data  structure  for  performing  operations  TEST  and  LIST  in  a  2-connected  graph 
(see  the  previous  subsection),  for  each  block  B  of  G. 

(3)  OUTER  pointers  for  the  blocks,  vertices,  and  faces  of  G. 

From  Theorem  7.3,  we  conclude 

Theorem  7.4  There  exists  a  data  structure  that  solves  the  topological  location  problem  for  gen¬ 
eral  plane  graphs  with  the  following  performance: 

Space  (m)= Preprocess  ( m)-0 (m); 

Test  (m)~0 (log m);  List  ( m,k)=0 (log  m  +  k). 


7.4.3.  Average  query  time  Let  Quv(n,m)  be  the  time  to  perform  the  query  operation 
TEST(u,v),  We  have  shown  in  the  previous  subsection  that  in  the  worst  case 
QUv(n,m)  =  O  (log  n).  Here,  we  consider  the  average  query  time  over  all  possible  queries. 


defined  by:  Q(n,m)  = 


u,v 

u*v 


Theorem  7.5  In  the  previously  described  data  structure  for  the  topological  location  problem,  the 


average  time  for  the  TEST  operation  over  all  possible  queries  is  Q(n,m)  =  0 


log  — 
n 


.  In 
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particular,  if  the  graph  is  simple,  then  Q(n,m)=0  (1). 

Proof:  From  the  description  of  the  algorithm  for  the  TEST  operation  we  have  that 
Quv(n,m)  =  O (log deg" (u)  +  log deg“(v)).  Hence,  fljl  Q(n,m)  =  G(n£logdeg_(v)).  The  proof 

^  J  V 

is  completed  observing  that  Y  deg~(v)  -m. 

V  □ 


IS.  Dynamic  Planar  Graph  Embedding 

In  this  section,  we  present  the  complete  data  structure  for  the  dynamic  embedding  problem. 
First,  we  consider  the  problem  in  spherical  sr-graphs,  and  then  extend  the  results  to  undirected 
graphs  using  spherical  sr-orientations. 

7.5.1.  Spherical  s/-graphs  In  this  subsection  we  describe  a  data  structure  for  efficiently  solv¬ 
ing  the  dynamic  embedding  problem  for  spherical  st-  graphs.  Let  G  be  a  spherical  rr-graph. 

The  data  structure  has  a  record  for  each  vertex,  edge,  and  face  of  G.  The  records  for  the 
vertices  are  arranged  in  a  vertex-tree  Tv,  which  is  a  balanced  search  tree  whose  nodes  are 
ordered  according  to  the  alphabetic  order  of  the  names  of  the  vertices.  Similarly,  the  records  of 
the  edges  and  faces  are  arranged  in  alphabetic  order  in  a  face-tree  Tp,  and  in  an  edge-tree  Tp, 
respectively.  The  above  trees  allow  us  to  access  in  G(logm)  time  the  records  associated  with 
the  vertices,  edges,  and  faces  involved  in  the  current  operation. 

The  record  for  a  face  /  stores  the  following  information: 

(1)  /:  name  of  the  face; 

(2)  HIGH (J):  pointer  to  the  record  of  the  topmost  vertex  of  /, 

(3)  LOW(f ):  pointer  to  the  record  of  the  bottommost  vertex  of  /, 

(4)  pointers  to  two  balanced  search  trees,  LPATH(f)  and  RPATH(f),  associated  with  the  left 
and  right  paths  of  face  /,  respectively.  The  nodes  of  each  such  tree  represent  the  vertices 
and  edges  of  the  corresponding  path,  and  are  sorted  according  to  the  direction  of  the  path. 
The  roots  of  LPATHif)  and  RPATH(f)  point  back  to  the  record  of  /. 
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The  record  for  a  vertex  v  stores  the  following  information: 

(1)  v:  name  of  the  vertex; 

(2)  deg+(v),  deg"(v):  outdegree  and  indegree  of  v. 

(3)  PLEFT(y ):  pointer  to  the  representative  of  v  in  the  tree  RPATH(f),  where  / = LEFT(v ); 

(4)  PRIGHTiy):  pointer  to  the  representative  of  v  in  the  tree  LPATH(g),  where  g  -  RIGHT {y)\ 

(5)  BELOW (v):  pointer  to  a  balanced  search  tree  whose  nodes  store  pointers  to  the  records  of 
the  faces  /  such  that  v  =HIGH(f).  The  nodes  of  BELOW (v)  are  sorted  according  to  the 
alphabetic  order  of  the  names  of  the  vertices  LOW(f). 

The  record  for  an  edge  e  stores  the  following  information: 

(1)  e:  name  of  the  edge; 

(2)  pointers  to  the  records  of  the  vertices  LOW(e)  and  HIGH(e); 

(3)  pointers  PLEFT(e)  and  P RIGHT  (e)  to  the  representatives  of  e  in  the  trees  RPATH(f)  and 
LPATH{g),  where  f=LEFT(e)  and  g  =  RIGHT  (e). 

We  show  in  Fig.  7.8  a  spherical  si-graph  and  a  fragment  of  the  data  structure  for  it. 

Using  the  above  data  structure,  the  TEST  operation  can  be  performed  with  the  same  stra¬ 
tegy  as  in  the  static  case.  The  only  difference  is  that  now  the  faces  to  the  left  and  right  of  u  and  v 
are  not  immediately  available,  and  must  be  retrieved  by  walking  up  to  the  roots  of  the  trees  that 
contain  the  representatives  of  u  and  v  pointed  to  by  PLEFT(u),  P RIGHT (u),  PLEFTiy),  and 
PRIGHTiy ). 

With  regard  to  the  INSERTEDGE  operation,  testing  for  its  applicability  can  be  done  by  a 
simple  modification  of  the  TEST  algorithm. 

Now,  assume  that  f\  is  to  the  left  of  e  and  fi  is  to  the  right  of  e.  We  partition  the  left  path 
of  /into  subpaths  L\,  Li,  and  £3,  where  Li  is  the  left  path  of  f\.  Notice  that  L\  and/or  L^ 
might  be  empty.  Analogously,  we  partition  the  right  path  of  / into  subpaths  R\,  R 2,  and  R 3, 
where  R  2  is  the  right  path  of  fi-  After  the  insertion  of  edge  e,  the  new  boundaries  of  f\  and  fi 
are  as  follows  (see  Fig.  7.9): 
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vertex  record 


Figure  7.8  A  spherical  tf-graph  and  a  fragment  of  the 
dynamic  data  structure  for  it.  The  edges  are  oriented 
upward. 


LPATH(fO=L2-,  RPATH[f\)=R\eRy,  LPAT//(/2)=L1eL3;  RPATH(f2)=R2. 

Hence,  the  LPATH  and  RPATH  trees  of  the  new  faces  are  obtained  by  splitting  LPATHif)  and 
RPATHif)  at  vertices  u  and  v,  and  splicing  appropriately  the  resulting  trees.  Standard  tech¬ 
niques  allow  us  to  perform  the  above  split  and  splice  operations  in  logarithmic  time  (see  for 
example  (39,  pp.  213-216]). 

The  remaining  updates  of  the  data  structure  are  as  follows: 

(1)  Insert  into  Te  a  new  record  for  edge  e,  and  increment  counters  deg+(u)  and  deg~(v). 

(2)  Delete  from  7>  the  record  of  face  f,  and  insert  new  records  for  faces  f\  and  f2. 

(3)  Delete  from  BELOW  {HIGH  if))  the  node  pointing  to  /,  and  insert  into  BELOW (HIGH  (fx )) 
and  BELOW (HIGH{f2))  nodes  pointing  to/j  and  /2,  respectively. 

The  INS ERTVERTEX  (e,v,e  1,62)  operation  is  performed  as  follows: 

(1)  Delete  from  Te  the  record  for  edge  e,  and  insert  new  records  for  edges  e  1  and  e2. 


Figure  7.9  Restructuring  of  the  face  boundaries  after  an 
INSERTEDGE  operation,  (a)  u  and  v  on  left  path;  (b)  u  on 
left  path  and  v  on  right  path;  (c)  u  bottommost  and  v  on 
right  path;  (d)  u  bottommost  and  v  topmost. 
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(2)  Insert  into  Ty  a  new  record  for  vertex  v. 

(3)  Find  the  faces  /  and  g  respectively  to  the  left  and  right  of  edge  e  by  walking  up  to  the  roots 
of  the  trees  that  contain  PLEFT(e)  and  P RIGHT (e). 

(4)  Replace  the  leaf  representative  of  e  in  RPATH(f )  by  a  subtree  with  root  v  and  children  e  i 
and  «2.  and  rebalance  RPATH(f). 

(5)  Replace  the  leaf  representative  of  e  in  LPATH{g)  by  a  subtree  with  root  v  and  children  e  \ 
and  &  2.1  and  rebalance  LPATH{g). 

(6)  Set  BELOW (v)  :=  0. 

The  above  data  structure  also  supports  operations  REMOVEEDGE  and  REMOVEVERTEX. 
In  fact,  the  REMOVEEDGE  operation  can  be  performed  by  reversing  the  transformations  on  the 
data  structure  realized  by  the  algorithm  for  the  INSERTEDGE  operation.  Similarly,  the  REMO¬ 
VEVERTEX  operation  is  the  reverse  of  the  INSERTVERTEX  operation.  Notice  also  that,  by 
Lemma  7.3,  the  counters  deg+(v)  and  deg“(v)  allow  us  to  test  the  feasibility  of  each  such  opera¬ 
tion  in  0(1)  time.  The  time  complexity  analysis  of  the  various  operations  is  straightforward, 
and  we  conclude 

Theorem  7.6  There  is  a  data  structure  for  the  dynamic  embedding  problem  in  spherical  st- 
graphs  with  the  following  performance: 

Space  (m ) = Preprocess  {m)~0{m)\ 

Test  (m)=0(logm);  List  (m,k)  =  0(\ogm  +k); 

InsertEdge  (m )  =  InsertVertex  ( m)=RemoveEdge  (m ) = RemoveVertex  (m)  =  0(\ogm). 

In  the  execution  of  an  update  operation  we  can  distinguish  the  search  time  spent  in  finding 
the  nodes  of  the  various  trees  involved  in  the  operation,  and  the  restructuring  time  that  takes  into 
account  the  update  and  rebalancing  of  the  trees.  The  next  theorem  shows  that  in  our  data  struc¬ 
ture  the  amortized  restructuring  time  for  a  sequence  of  INSERTEDGE  and  INSERTVERTEX 
operations  is  optimal.  For  the  definition  of  amortized  time  complexity,  see  [56]. 

Theorem  7.7  There  exists  a  data  structure  for  the  dynamic  embedding  problem  in  spherical  st- 


(1)  The  space  occupation  and  time  complexity  of  the  various  operations  are  the  same  as  in 
Theorem  7.6. 

(2)  In  a  sequence  of  INSERTEDGE  and  INSERTVERTEX  operations,  the  amortized  restructur¬ 
ing  time  complexity  of  each  such  operation  is  0(1). 

Proof:  Use  2-4  trees  (which  are  equivalent  to  red-black  trees  [20])  to  realize  trees  Tp,  Ty,  Tp, 
and  BELOWiy).  Such  trees  have  0(1)  amortized  rebalancing  time  for  insertions  and  deletions 
[26].  With  regard  to  the  LPATH  and  RPATH  trees,  their  manipulation  in  a  sequence  of  INSER¬ 
TEDGE  and  INSERTVERTEX  operations  involves  insertions  and  the  kind  of  generalized  split¬ 
tings  considered  in  [23].  It  is  shown  there  that  circular  level-linked  2-4  trees  support  efficiently 
a  sequence  of  insertions  and  generalized  splittings.  With  arguments  similar  to  the  ones 
developed  in  [23]  we  can  show  that  by  realizing  the  LPATH  and  RPATH  trees  by  circular  level- 
linked  2-4  trees,  the  amortized  rebalancing  time  for  LPATH  and  RPATH  is  O  (1).  □ 

7.5.2.  s/-orientable  graphs  For  sr-orientable  plane  graphs,  we  maintain  on-line  a  spherical  st- 
orientation  and  use  the  data  structure  previously  described.  Operations  TEST,  LIST,  and 
INSERTVERTEX  do  not  require  any  modifications.  In  connection  with  operation 
INSERTEDGE  (e,u,v,f-Ji  ,/2),  we  have  to  select  the  direction  of  edge  e  so  that  it  does  not  intro¬ 
duce  cycles  internal  to  face  /.  By  Corollary  7.1,  this  can  be  done  easily  by  reversing  the  direc¬ 
tion  whenever  the  INSERTEDGE  algorithm  rejects  the  operation.  This  proves 

Theorem  7.8  There  is  a  data  structure  that  supports  operations  TEST,  LIST,  INSERTEDGE,  and 
INSERTVERTEX  in  a  sf-orientable  plane  graph  with  the  following  performance: 

Space  (m)  =  Preprocess  (m)  =  0(m); 

Test  (m)  =  0(log  m);  List  ( m,k )  =  0(logm  +  k)\ 

InsertEdge  {m)~InsertVertex  (m)  =  0(logm). 

In  regard  to  the  REMOVEEDGE  and  REMOVEVERTEX  operations,  we  are  faced  with  the 
difficulty  that  the  data  structure  acts  on  a  spherical  jr-orientation  of  the  graph,  therefore  permit¬ 
ting  only  deletions  that  preserve  the  sr-structure  of  the  orientation.  We  say  that  a  vertex  (edge)  is 
free  if  operation  REMOVEVERTEX  ( REMOVEEDGE )  can  be  performed  on  it  in  the  spherical 


sf-orientation;  we  say  that  it  is  locked  otherwise.  At  any  time  the  vertices  and  edges  of  the  graph 
are  partitioned  into  free  and  locked,  and  we  are  allowed  to  delete  only  vertices  and  edges  that  are 
free.  We  thus  have 

Theorem  7.9  The  data  structure  of  Theorem  7.8  supports  also  operations  REMOVEEDGE  and 
REMO V EVERT EX  on  free  edges  and  vertices  in  0  (logm)  time. 

In  several  layout  applications,  design  methodologies  limit  the  freedom  of  the  designer  in 
making  arbitrary  updates  to  the  layout.  For  example,  well  known  hierarchical  design  strategies 
for  VLSI  circuits  build  a  layout  in  a  top-down  fashion  by  means  of  successive  refinements. 

In  the  following,  we  show  that  the  class  of  free  edges  and  vertices  is  sufficiently  large  to 
support  a  hierarchical  deletion  scheme  that  allows  us  to  “undo”  any  INSERTEDGE  and 
INSERTVERTEX  operation  performed  in  the  past  (not  just  the  last  operation).  For  instance,  we 
define  the  hierarchical  embedding  problem  as  a  variation  of  the  previously  discussed  dynamic 
embedding  problem,  where  the  following  restrictions  are  placed  on  the  REMOVEEDGE  and 
REMOVEVERTEX  operations: 

(1)  An  edge  can  be  deleted  by  a  REMOVEEDGE  operation  only  if  it  was  created  (at  any  time 
in  the  past)  by  means  of  an  INSERTEDGE  operation. 

(2)  A  vertex  can  be  deleted  by  a  REMOVEVERTEX  operation  only  if  it  was  created  (at  any 
time  in  the  past)  by  an  INSERTVERTEX  operation. 

In  Fig.  7.10  we  show  a  sequence  of  update  operations  in  an  instance  of  the  hierarchical  embed¬ 
ding  problem. 

The  aforementioned  restrictions  on  the  REMOVEEDGE  and  REMOVEVERTEX  operations 
can  be  enforced  by  storing  with  each  edge  e  two  flags,  denoted  FREE -TAIL(e)  and 
FREE -HEAD  (e),  which  are  associated  with  the  head  and  tail  of  edge  e  in  the  spherical  st¬ 
onemason,  respectively.  We  use  these  flags  to  maintain  the  invariant  that  an  edge  e  can  be 
removed  if  and  only  if  both  FREE-TAIL(e )  and  FREE-HEAD(e)  are  set.  The  manipulation  of 
the  flags  in  the  various  operations  is  straightforward.  In  the  example  of  Fig.  7.10,  a  flag  in  the 
“set”  condition  is  indicated  by  showing  the  corresponding  endpoint  unconnected. 
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It  is  not  difficult  to  show  that  the  free  edges  and  vertices  include  the  ones  that  can  be 
deleted  in  an  instance  of  the  hierarchical  embedding  problem.  We  have  thus  the  following 
theorem: 

Theorem  7.10  There  exists  a  data  structure  that  allows  us  to  solve  the  hierarchical  embedding 
problem  for  st-orientable  plane  graphs  with  the  following  performance: 

Space  (m)= Preprocess  (m)=0(m); 

Test  (m)  =  0(log  m);  List  ( m,k )  =  0(logm  +  k); 

InsertEdge  (m)  =InsertVertex  ( m)=RemoveEdge  (m)  =RemoveVertex  (m)  =  0(log  m). 

7.5.3.  General  plane  graphs  For  general  plane  graphs  we  add  to  the  repertory  the  following 
operations,  which  allow  respectively  to  add  and  remove  an  elementary  block  to  the  graph: 

ATTACH  (e,v,u,f):  Add  vertex  v  and  edge  (u,v)  inside  face  /. 

DETACH  (e,v,u,f):  Remove  the  degree-1  vertex  v  and  its  incident  edge  e=(u,v),  which 
lie  in  face  /. 

In  the  following,  we  will  provide  an  amortized  O(logm)  time  bound  for  operation  INSER¬ 
TEDGE,  and  a  worst-case  Oflogm)  time  bound  for  the  remaining  operations.  According  to 
standard  conventions  in  amortized  complexity  analysis  [56],  this  means  that  a  sequence  of  m 
operations  starting  from  an  initial  graph  consisting  of  a  single  edge  takes  O  (m  logm)  time. 

First,  we  dynamize  the  data  structure  that  tests  whether  two  vertices  belong  to  the  same 
block  of  a  graph  G.  We  represent  the  oriented  block-cutvertex  tree  T  of  G  by  means  of  a  bal¬ 
anced  search  tree  for  the  children  of  each  node.  Also,  we  store  with  each  block  B  a  pointer 
PROPER  (B )  to  a  balanced  search  tree  whose  nodes  represent  the  proper  vertices  of  B.  In  turn, 
each  proper  vertex  v  of  B  has  a  pointer  to  its  representative  in  the  tree  PROPER  (B).  We  have 

Lemma  7.7  There  exists  a  data  structure  with  0(m )  space  and  preprocessing  time  that  allows 
for  testing  whether  two  vertices  belong  to  the  same  block  in  O(logm)  time  (worst-case),  and 
supports  operations  INSERTVERTEX  and  ATTACH  in  O(logm)  worst-case  time,  and  operation 
INSERTEDGE  in  O(logm)  amortized  time. 
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We  now  show  how  to  dynamically  maintain  the  OUTER  relation.  For  each  block  B,  we 
define  INT(B)  as  the  list  of  the  cutvertices  of  B  such  that  OUTER(y)=B.  Also,  we  define  a  list 
EXT(B)  for  each  block  B  and  a  circular  list  INNER  (f)  for  each  cutface  /  by  means  of  the  follow¬ 
ing  constructive  procedure. 

Consider  a  cutface  /  whose  outer  block  is  Bq.  Initially,  all  cutvertices  of  /  in  B  0  are 
marked,  while  the  remaining  cutvertices  of  /  are  unmarked.  Also,  all  the  inner  blocks  of  /  are 
initially  unmarked.  Suppose  now  that  we  traverse  the  boundary  of  /  counterclockwise,  starting 
at  a  cutvertex.  Whenever  we  traverse  an  edge  of  an  unmarked  block  B,  we  add  (a  representative 
of)  B  into  INNER  (f)  and  mark  B.  Also,  whenever  we  encounter  an  unmarked  cutvertex  v,  we 
mark  v  and  add  (a  representative  of)  v  into  the  list  EXT(B)  of  the  first  block  B  in  the  current 
INNER  if)  that  contains  v.  At  the  end  of  the  traversal,  INNER  (v)  will  contain  all  the  inner 
blocks  of  /,  and  the  lists  EXT(B)  will  contain  all  the  cutvertices  v  of  /such  that  f -OUTER  (v). 
Notice  that  each  block  B  is  contained  in  exactly  one  INNER  if)  list,  and  each  cutvertex  v  is  con¬ 
tained  in  exactly  one  EXT(B)  or  INT(B)  list 

In  the  example  of  Fig.  7. 1 1(a)  we  have 

INNER(f)*(B2,Bl,B3,B5,B4,B6.B9,B7,Bsy, 

£XT(£2)=(vi,v2);  £XT(fl3)  =  (v4);  £XT(B5)  =  (v5,v6);  £XT(fl7)  =  (v7); 

EXT(Bx)=EXT(BA)  =  EXT(B6)=EXTiBs)=EXT(B9)  =  0. 

Implementing  the  above  lists  by  means  of  balanced  search  trees,  OUTER  i  )  can  be  easily 
retrieved  in  O(logm)  time,  so  that  the  complexity  of  the  TEST  and  LIST  operations  is  the  same 
as  in  the  static  case. 

Now,  we  describe  the  transformations  of  the  data  structure  due  to  the  execution  of  opera¬ 
tion  INSERTEDGE  (e,n,v,/;/i,/2).  If  face  /  is  not  a  cutface  and  is  part  of  block  B,  we  simply 
perform  the  operation  in  the  data  structure  of  5.  Otherwise,  edge  e  and  a  sequence  of  blocks  of 
G,  called  old-blocks  are  merged  into  a  new  block,  called  new-block.  In  the  example  of  Fig. 
7.11(b),  the  old-blocks  are  B\,B2,  and  £3,  and  the  new-block  is  the  union  of  B\,  B2,  B3,  and 
edgee=(u,v). 

Specifically,  let  p(w)  be  the  node  of  the  block-cutvertex  tree  T  of  G  associated  with  vertex 
w  (i.e„  if  w  is  a  cutvertex  then  |i(w)  is  the  node  representing  w,  and  otherwise  p(w)  is  the  node 
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representing  the  block  of  w).  The  old-blocks  are  exactly  the  ones  on  the  path  of  T  from  |i(u)  to 
p(v).  Also,  in  general  both  f\  and  f2  will  be  cutfaces  in  the  updated  embedding. 

In  conclusion,  the  modification  of  the  block  structure  caused  by  the  INSERTEDGE  opera¬ 
tion  requires  that  we  set  up  a  new  data  structure  for  the  new-block  B,  and  new  INNER  lists  for/i 
and/2. 

Now,  we  show  that  the  data  structure  for  B  can  be  set  up  by  reorienting  all  but  one  old- 
block,  called  base-block,  which  is  chosen  as  an  old-block  with  maximum  number  of  edges. 

Indeed,  let  the  old-blocks  be  B 1,  •  *  *  ,B[,  with  u  e  B\  and  ve  5/,  and  let  v,-,  1  =  1,  •  •  •  ,/-l, 
be  the  cutvertex  between  Bt  and  2?l  +  1.  Also,  let  v0  =  u  and  v/=v.  (In  the  example  of  Fig. 
7.11(b)  1  =  3.)  Denoting  by  Bj  the  base-block,  at  least  one  of  the  directed  edges  (v;_ltv;)  and 
(Vj,Vj- 1),  say  can  be  added  to  Bj  while  preserving  the  spherical  sr-orientation.  In  this 

case,  we  construct  a  new  acyclic  spherical  jr-orientation  for  each  £,,  i*j,  using  v,„i  as  the 
source  and  vf-  as  the  sink,  and  we  direct  e  from  v/  to  v0.  (The  case  when  (vj-\,Vj)  can  be  added 
to  Bj  is  analogous.)  The  resulting  orientation  of  the  new-block  B  is  a  spherical  tf-orientation, 
from  which  the  data  structure  for  B  can  be  constructed  in  time: 

/ 

O  (log  m  +  mi  -  max  mi), 

where  m,  is  the  number  of  edges  of  fl;. 

Lists  INNER (f\)  and  INNER (f  2)  are  obtained  from  INNER (f)  by 

(1)  removing  the  (representatives  of  the)  old-blocks; 

(2)  splitting  the  resulting  list  into  INNER  ( f  \ )  and  INNER  (fi)  with  an  appropriate  cut;  and 

(3)  possibly  adding  to  one  of  these  lists  (a  representative  of)  the  new-block  B. 

With  regard  to  the  EXT  lists,  we  form  INT(B)  and  EXT{B)  by  a  sequence  of  O  (/)  split  and 
splice  operations  on  the  EXT  lists  of  the  old-blocks. 

In  the  example  of  Fig.  7.1 1,  we  have 

INNER (f \)  =  (B,B 9, By, B%);  INNER(f2)=(Bs,BA,B6y,  £XT(fi)  =  (v2);/Nr(B)  =  (v1,v4). 
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Theorem  7.11  There  is  a  data  structure  that  supports  operations  TEST,  LIST,  INSERTEDGE, 

INSERTVERTEX ,  and  ATTACH  with  the  following  performance: 

Test  (m)  =  O (log m);  List  (m,k)  =  0 (log m  +k); 

InsertEdge  (m)  =  0(logm)  (amortized);  InsertVertex  (m)=Attach(m)=0  dog  m). 

Proof:  The  amortized  complexity  analysis  makes  use  of  the  following  potential  function : 

b  1 

0  =  £  m,  log  — 

i  =  1  m» 

b 

where  B  x ,  •  •  •  ,2?*  are  the  current  blocks  of  G,  and  mx  is  the  number  of  edges  of  B,  ( £  m,  =/n). 

t  =  1 

Notice  that  -mlogm^O^O.  It  is  interesting  to  observe  that  the  function  O  is  similar  to  the 
entropy  function  used  in  information  theory.  Informally,  we  can  say  that  when  the  old-blocks 
are  merged  together  into  the  new-block,  O  decreases  to  compensate  for  the  work  spent  in  the 
merge  process.  For  simplicity,  consider  the  case  when  two  old-blocks,  B  i  and  B 2,  are  merged 
into  a  new-block  B* ,  where  mx  is  the  number  of  edges  of  B,  (i  =  1,2)  and  m*  is  the  number  of 
blocks  of  B*.  By  an  appropriate  choice  of  the  time  unit,  we  have  that  operation  INSERTEDGE 
is  executed  in  time 


r  =  logm  +  2min  [m\ ,  m2). 
Now,  the  variation  AO  of  potential  is  given  by 


AO=m*log— 

m 


m  1  log  +m  2  log  -^—1 . 
rn  1  m  2] 


m  1  m  2 

Define x \  = — -  and *2  =  — 7-  Clearly,  X\  +x2  =  1.  We  can  write 
m  m 


AO  .  ^  1 

—  =  (x1+x2)\og  —  - 
m  m 


X\  log  — + x 2  log  — -  I 
m  1  m2 


Xi  log  — +x2log  — 
X\  X2 


-H(x\pt2). 


where  H{xx,x2)  is  the  binary  entropy  function  of  information  theory.  We  show  in  Fig.  7.12  a 
plot  of //(*!,  1-xj)  in  the  interval  xx  s  [0,1].  It  is  easy  to  see  that 
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Figure  7.12  Plot  of  the  binary  entropy  function  H(xi,l-xi). 


Therefore,  we  obtain 


//(*!,! -*!>£  2  min  (xi,!-*!). 


-2  m*  min  =  -2 min  [mltm 2). 


Hence,  the  amortized  time  complexity  of  operation  INSERTEDGE  is 

7=r+A0^1ogm+2min  {mi,/n2)  ~2min  [mi,m2)  =log m.  q 

With  regard  to  operations  REMOVEEDGE ,  REMOVEVERTEX,  and  DETACH,  again  we 
can  partition  the  edges  of  the  graph  into  free  and  locked,  and  we  have  the  following: 

Corollary  7.3  The  data  structure  of  Theorem  7.11  also  supports  operations  REMOVEEDGE, 
REMOVEVERTEX,  and  DETACH  on  free  edges  and  vertices  in  0  (logm)  time  (worst-case). 

7.6.  Applications 

7.6.1.  A  dual  embedding  problem  The  dual  dynamic  embedding  problem  consists  of  perform¬ 
ing  the  following  operations  on  a  plane  graph  G : 
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TEST*  (f,g ):  Test  if  there  is  a  vertex  v  that  is  on  the  boundaries  of  both  faces  /  and  g.  In 
case  such  a  vertex  v  exists,  output  its  name. 

LIST*  (f,g ):  List  all  the  vertices  that  are  on  the  boundaries  of  both  faces  / and  g. 

EXPANDVERTEX  (e,/,g,v;v  \ ,V2):  Expand  vertex  v  into  vertices  v  i  and  v2  connected  by 
an  edge  e  on  the  boundary  of  faces  /and  g. 

CONTRACTVERTEX  (e,f,g,vx,v 2;v):  Contract  the  edge  e  =  (v1,v2)  on  the  boundary  of 
faces  /  and  g,  and  call  v  the  vertex  resulting  from  the  contraction  of  v  x  and  v  2. 

DUPL1CATEEDGE  (,e,f\e i,c2):  Replace  the  edge  e  with  two  multiple  edges,  ex  and  e2, 
with  the  same  endpoints,  and  call  /  the  resulting  face  between  them. 

MERGEEDGE  (ex,e2,f;e):  Let  /  be  a  face  whose  boundary  consists  of  two  multiple 
edges,  ex  and  e2.  Remove /by  merging  ex  and  e2  into  a  new  edge  e. 

We  can  show  that  this  problem  is  the  dual  of  the  dynamic  embedding  problem  by  extending 
the  notion  of  duality  of  undirected  plane  graphs  to  spherical  sr-graphs. 

The  dual  graph  G*  of  a  plane  digraph  G  is  the  plane  digraph  defined  as  follows:  the  ver¬ 
tices  of  G  *  are  the  faces  of  G;  for  each  edge  e  of  G  there  is  an  edge  e  *  of  G  *  from  face  LEFT(e) 
to  face  RIGHT (e).  A  cylindrical  s*t* -graph  is  a  plane  digraph  G  whose  dual  G*  is  a  spherical 
sr-graph.  An  (undirected)  plane  graph  is  said  to  be  s* t* -orientable  if  it  can  be  oriented  to  become 

a  cylindrical  sY-graph.  Notice  that  every  2-connected  graph  is  sY-orientable.  Using  the 
results  of  Section  7.4  and  duality  arguments,  we  obtain 

Theorem  7.12  There  is  a  data  structure  that  allows  us  to  solve  the  dual  dynamic  embedding 

problem  for  cylindrical  sY-graphs  with  the  following  performance: 

Space  (m )  =  Preprocess 

Test*  ( m )  ~  G(log  m );  List*  ( m,k ) = G(log  m  +  k); 

ExpandVertex  (m )  =DuplicateEdge  (m )  =  ContractVertex  (m ) = MergeEdge  ( m)  =  0  (log  m). 

Further  results  can  be  obtained  by  dualizing  the  remaining  theorems  of  the  previous  sec¬ 
tion. 
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The  hierarchical  dual  embedding  problem  is  defined  as  a  variation  of  the  dynamic  embed¬ 
ding  problem  where  the  C ONTRA CTVERTEX  and  MERGEEDGE  operations  can  performed 
only  on  edges  (respectively,  faces)  previously  inserted  by  the  EXPANDVERTEX  (respectively, 
DUPUCATEEDGE)  operation.  We  have 

Theorem  7.13  There  is  a  data  structure  that  allows  us  to  solve  the  hierarchical  dual  embedding 

problem  for  s Y-orientable  plane  graphs  with  the  following  performance: 

Space  (m) = Preprocess  ( m)=0 (m); 

Test*  ( m)=0  (log  m);  List*  ( m,k )  -  O  (log  m + Jfc); 

ExpandVertex  (m ) =DuplicateEdge(m)  -  ContractVertex(m)=MergeEdge(m)  =  0(logm). 

7.6.2.  Computing  separating  pairs  Let  G  be  a  2-connected  plane  graph.  A  separating  pair  of 
graph  G  is  a  pair  of  distinct  vertices  of  G  whose  removal  disconnects  G.  The  identification  of 
separating  pairs  is  important  in  problems  of  fault-tolerance  of  networks. 

Lemma  7.8  Let  G  be  a  2-connected  plane  graph.  Vertices  u  and  v  form  a  separating  pair  of  G  if 
and  only  if 

(1)  u  and  v  are  adjacent  and  there  are  at  least  three  faces  whose  boundary  contains  both  u  and 
v  ;  or 

(2)  u  and  v  are  not  adjacent  and  there  are  at  least  two  faces  whose  boundary  contains  both  u 
and  v. 

Hence,  we  can  test  whether  two  vertices  form  a  separating  pair  by  a  simple  modification  of 
the  algorithm  for  operation  LIST  (u,v),  which  halts  as  soon  as  two  or  three  faces  have  been 
reported.  We  obtain 

Corollary  7.4  The  data  structure  of  Theorem  7.10  allows  us  to  test  whether  two  given  vertices 
form  a  separating  pair  in  0(1 og  m)  time. 
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